
Citation: | LIU H Y,CHEN W T,LI Y Y,et al.,2024. Application of integrated model based on EEMD-CNN-LSTM for landslide-displacement prediction[J]. Journal of Geomechanics,30(4):633−646 doi: 10.12090/j.issn.1006-6616.2023145 |
人工地震勘探是探测深部构造的常用方法, 但这种方法只能探测到两盘具有不同阻抗值的断层。而通过监测断层活动产生的微地震, 不但可以了解断层的空间位置, 还能够查明断层的活动规律。近年来在油田开发中广泛应用的微地震监测技术, 就是借助天然地震学方法和思路来监测微地震的空间和时间分布规律, 找出微地震与深部构造活动的关系, 并据此推断构造的活动情况。
当应力增大或屈服强度降低时, 构造薄弱面上一些局部部位会首先因材料失稳而产生裂缝及其微弱活动。水的注入等因素可以导致断层本身屈服强度降低。在裂缝活动过程中, 部分能量通过地震弹性波的方式释放, 并在之后达到新的平衡。当失稳范围和规模从小到大, 就可能发生地质灾害或较大的地震。利用观测到的裂缝微地震活动发出的高频波动信号可以获得发生破裂的位置和特征[1]。长期地震观测显示, 很多地震是由断裂活动引发的[2]。所以通过对微地震活动的监测, 可以得到断裂活动的有关信息。本文就是通过对微地震活动的分析来获得断裂活动信息的。
与人工地震勘探方法相比, 天然地震监测中的震源位置、发震时刻和震源强度都是未知的。与油田常用的微地震监测技术相比, 本文讨论的微地震主要是由布设在地表的地震台网来记录的, 而非像油田常用微震监测那样将探头置于井下[3~4]。
微地震监测研究主要包括两部分:①通过层析成像反演得到地下深部的波速结构; ②探知微地震发生的位置及性质。两者是相辅相成的。因为只有在可靠的深部波速结构的基础之上, 才能得到可靠的地震发生位置及性质。反之亦然。
地震层析成像是一种利用地震波对地球内部三维结构进行成像以发现地下波速异常的方法。自从30多年前Aki等[5]首次把天然地震层析成像技术应用于研究三维地球结构之后, 该技术在大尺度的地球深部结构研究中得到了广泛的应用。天然地震层析成像技术在矿产资源勘探领域的应用是近些年才兴起的。
按照地震波射线传播理论, 地震波沿着波的传播路径到达地震台站所需要的时间可以表示为对波的传播路径上的走时积分。为了计算上的简便, 通常需要把研究区划分为离散的网格, 并假设每个网格内的波速是均匀的。这样台站到震源的走时积分就可以简化为走时叠加。对于很多的震源和台站, 可以建立很多路径的走时叠加方程。而地震层析成像的目的正是要通过这些方程组中已知的其中两个量求第三个量。比如已知地震波走时和地震震源位置求速度分布, 或已知地震波走时和速度分布求地震震源位置。所以反演震源位置和波速结构这两个过程是相辅相成的, 实际上就是一个迭代反演解大型方程组的问题。
微地震层析成像[6~7]是利用密集地震台阵来监测岩层中微地震活动所产生的地震波, 然后运用地震学方法来反演微地震的活动特征和研究区三维横波、纵波波速分布的一种勘探方法。它主要利用被动震源, 也可以利用人工地震, 并适用于各种复杂的地区和环境。
本文所用的地震检波器的主要监测频率介于5~100Hz之间, 而油田监测的地震信号频率通常高达200~2000Hz。
本文研究区(见图 1)是柴达木盆地西部狮子沟油田100km2的范围, 在该区布设的16个台站自2004年10月开始记录了大量连续的高质量的地震资料。地震波形实例见图 2, 地震台站位置见图 3。本文分析处理了2004年10月到2005年12月间共14个月的数据, 最终得到了研究区内952个微地震数据。
在数据处理过程中, 首先从连续地震记录中挑出研究区内信噪比高的地震事件波形, 然后对挑出的事件进行P波和S波震相到时拾取。图 2显示了8个地震台站记录到的同一个微地震的地震波, 数据信噪比高且震相到时清楚。图中标识为虚竖线的位置为首先到达的纵波P波震相, 标识为点竖线的位置为后到达的横波S波震相。
研究区西南部分为地势平坦的盆地, 东北部地区海拔较高, 中部地形变化较大(图 3)。西南部存在一个走向近北西的明显的盆山边界, 被认为是花土沟断裂的地表出露位置(图 3中虚线标识其大概位置)。
图 3显示了14个月时间内记录到的大于-2级地震的二维分布。震中分布的一个明显特点是, 地震主要集中在研究区中部海拔变化较大(由海拔3100m迅速增加到3300m或更高)、地形复杂的部位。另外, 研究区西北角的地震也较多, 这里正是花土沟断裂和一条北东东走向大沟的交会区。由于盆山交界区域附近没有地震台站, 在该区域观测到的地震数量较少。多数地震呈带状集中在花土沟断裂带北东侧, 并且主体走向与花土沟断裂平行(图 3)。由于微地震通常沿着断层分布, 因此推断花土沟断裂为北东倾向。
根据对较大地震(震级大于1级)的统计, 一般地, 一个地区一段时间内地震发生频次(N)的对数与地震震级(M)成近似线性关系(即log (N) =a-bM, 其中a和b为常数, b在地震学中称为b值)[8], 反映了该地区在一定时间段内的应力场特征。某个时期斜率常数b值发生的变化, 意味着区域应力场等性质发生变化, 因此, b值常常被用于地震预报中[9]。
图 4为观测到的微地震发生频次与地震震级关系分布图。本文-0.25到2级的微地震也表现出地震频次对数随震级增大而减小的线性关系, 但小于-0.25级的地震表现为频次随地震震级减小而减小的关系(图 4)。这说明, 现有设备可以记录所有震级大于-0.25的地震。但是, 可能由于设备及其部署等原因, 记录到的更小地震事件有所遗漏, 并随震级减小而遗漏增多。
均匀应力场下统计的b值通常小于1, 但构造应力场非常不均匀时可以产生很多地震群, 这时b值可以高达2.5[8]。研究显示[10~11], 钻孔注入流体或其他原因导致的孔隙压力变化同样可以触发局部微震。图 4数据显示的b值约为1.35, 略大于均匀应力场下的b值(< 1.0)。因此总体看来, 研究区应力场分布可能是均匀的。但是, 油田开采或钻孔注水等改变了孔隙压力, 从而触发了一些微震, 最终可能导致b值偏高。
油田生产或注水可以导致孔隙压力变化, 从而可能触发微震。如图 3左上角钻孔密度较大的地区(即钻孔SZ10, S11和S1等附近)的地震比周围偏多; 中部的S29和S28等钻孔较多的地区地震事件相对也较多。但是, 其他一些钻孔密集区域的地震事件并未明显增加。
图 5为穿过研究区的一条北东向垂直剖面(剖面位置见图 3中黑实线)。由于盆山交界(花土沟断裂出露位置)附近没有布设地震台站, 所以该区基本没有检测到地震。另外, 图 5中的地震分布相对离散。虽然利用相对定位方法(如近年常用的双差分方法[12])可以使定位后的地震分布相对集中, 但由于没有证据表明微地震事件之间具有关联性, 因此这里没有采用。
图 5显示, 地震多集中在2km深度左右, 而且总体(图 5中虚线标出的范围)具有北东倾向。按照图中地震总体分布向左上方延伸, 就可以与花土沟断裂出露位置基本重合。考虑到图 3中地震带主体走向与花土沟断裂平行, 因此可以认为剖面显示的地震主要与以花土沟断裂为主的断裂带有关。
本文通过对柴达木盆地西部狮子沟地区密集高频地震台网一年多的微地震观测数据的分析, 检测到大量微地震, 可能与花土沟断裂的深部地震活动有关。该微地震活动带为北西走向, 北东倾向, 与前期地震观测结果[13]相吻合。
地震事件发生频次和震级之间关系的分析显示, 现有设备在现有布置格局下可以监测到研究区内所有震级大于-0.25的地震事件。因此, 虽然地震震级与频次对数的线性关系是在研究较大震级地震的过程中得到的, 但本文证实它同样适用于震级大于-0.25的微震。本区微震分布的b值约为1.35, 略大于均匀应力场下b值的范围。这可能与油田开采和钻孔注水等改变了孔隙压力而触发了一些微震有关。
[1] |
AN D, SONG K, YI Z, et al., 2021. A prediction model for reservoir landslide step-like displacements using combined EEMD and RFR method[J]. Mountain Research, 39(1): 143-150. (in Chinese with English abstract
|
[2] |
BINIYAZ A, AZMOON B, SUN Y, et al., 2022. Long short-term memory based subsurface drainage control for rainfall-induced landslide prevention[J]. Geosciences, 12(2): 64. doi: 10.3390/geosciences12020064
|
[3] |
FENG Y, ZENG H N, TU P F, 2023. A time series decomposition method for landslide displacement based on sliding detection algorithm[J]. Journal of Changjiang River Scientific Research Institute, 41(3): 126-133, 147. (in Chinese with English abstract
|
[4] |
FU W T, HU D M, WANG S H, et al. , 2020. Landslide displacement prediction based on EEMD-GM-ELM model[C]//Proceedings of the 11th China satellite navigation annual conference: S01 satellite navigation industry application. Beijing: Academic Exchange Center of China Satellite Navigation System Management Office: 89-93. (in Chinese with English abstract
|
[5] |
HAN S, LIU M J, WU J B, et al., 2022. Risk assessment of slope disasters induced by typhoon-rainfall in the southeast coastal area, China: A case study of the Shiyang north slope[J]. Journal of Geomechanics, 28(4): 583-595. (in Chinese with English abstract
|
[6] |
HARIRI S, KIND M C, BRUNNER R J, 2021. Extended isolation forest[J]. IEEE Transactions on Knowledge and Data Engineering, 33(4): 1479-1489. doi: 10.1109/TKDE.2019.2947676
|
[7] |
JIANG H W, LI Y Y, ZHOU C, et al., 2020. Landslide displacement prediction combining LSTM and SVR algorithms: a case study of Shengjibao landslide from the three gorges reservoir area[J]. Applied Sciences, 10(21): 7830. doi: 10.3390/app10217830
|
[8] |
KANG E S, ZHAO Z X, MENG H D, 2022. Displacement prediction of dump slope based on EEMD-HW-PSO-ELM coupling model[J]. Gold Science and Technology, 30(4): 594-602. (in Chinese with English abstract
|
[9] |
KATTENBORN T, LEITLOFF J, SCHIEFER F, et al., 2021. Review on Convolutional Neural Networks (CNN) in vegetation remote sensing[J]. ISPRS Journal of Photogrammetry and Remote Sensing, 173: 24-49. doi: 10.1016/j.isprsjprs.2020.12.010
|
[10] |
KRKAČ M, ŠPOLJARIĆ D, BERNAT S, et al., 2017. Method for prediction of landslide movements based on random forests[J]. Landslides, 14(3): 947-960. doi: 10.1007/s10346-016-0761-z
|
[11] |
Aggarwal A, Alshehri M, Kumar M, et al., 2020. Landslide data analysis using various time-series forecasting models[J]. Computers & Electrical Engineering, 88: 106858.
|
[12] |
LEI Y G, LIN J, HE Z J, et al., 2013. A review on empirical mode decomposition in fault diagnosis of rotating machinery[J]. Mechanical Systems and Signal Processing, 35(1-2): 108-126. doi: 10.1016/j.ymssp.2012.09.015
|
[13] |
LI H J, XU Q, HE Y S, et al., 2018. Prediction of landslide displacement with an ensemble-based extreme learning machine and copula models[J]. Landslides, 15(10): 2047-2059. doi: 10.1007/s10346-018-1020-2
|
[14] |
LI L J, YAO X, ZHOU Z K, et al., 2022. The applicability assessment of Sentinel-1 data in InSAR monitoring of the deformed slopes of reservoir in the mountains of southwest China: A case study in the Xiluodu Reservoir[J]. Journal of Geomechanics, 28(2): 281-293. (in Chinese with English abstract
|
[15] |
LIAN C, ZENG Z G, YAO W, et al., 2013. Displacement prediction model of landslide based on a modified ensemble empirical mode decomposition and extreme learning machine[J]. Natural Hazards, 66(2): 759-771. doi: 10.1007/s11069-012-0517-6
|
[16] |
LIU F T, TING K M, ZHOU Z H, 2008. Isolation forest[C]//2008 Eighth IEEE international conference on data mining. Pisa: IEEE: 413-422.
|
[17] |
MA Y, YU B, HE Y X, et al., 2023. Rainfall Threshold and Development Characteristics of Shallow Landslides Induced by Rainfall: A Case Study of the "June 10th, 2019" Disaster in the Dajishan Area, Quannan County, Jiangxi Province[J]. Geology and Exploration, 59(5): 1065-1073. (in Chinese with English abstract
|
[18] |
NAN X C, LIU J F, ZHANG Y X, et al., 2023. Dynamic prediction of landslide displacement based on multi-source time series[J]. Pearl River, 44(4): 54-62. (in Chinese with English abstract
|
[19] |
PAK U, KIM C, RYU U, et al. , 2018. A hybrid model based on convolutional neural networks and long short-term memory for ozone concentration prediction[J]. Air Quality, Atmosphere & Health, 11(8): 883-895.
|
[20] |
PENG L, NIU R Q, 2011. Analysis on deformation characteristics and influential factors of Baijiabao landslide in the Three Gorges Reservoir area[J]. The Chinese Journal of Geological Hazard and Control, 22(4): 1-7. (in Chinese with English abstract
|
[21] |
QU J L, WANG X F, GAO F, et al., 2014. Noise assisted signal decomposition method based on complex empirical mode decomposition[J]. Acta Physica Sinica, 63(11): 110201. (in Chinese with English abstract doi: 10.7498/aps.63.110201
|
[22] |
SHANG M, XIONG D B, ZHANG H Q, et al., 2022. Landslide displacement prediction model based on time series and mixed kernel function SA-SVR[J]. Journal of Engineering Geology, 30(2): 575-588. (in Chinese with English abstract
|
[23] |
SHI X L, HAN X D, YANG X Y, et al., 2023. Factors inducing the Xigouwan landslide in the Three Gorges Reservoir area and the influence of antecedent precipitation[J]. Journal of Geomechanics, 29(2): 253-263. (in Chinese with English abstract
|
[24] |
SINGH P, KEYVANLOU M, SADHU A, 2021. An improved time-varying empirical mode decomposition for structural condition assessment using limited sensors[J]. Engineering Structures, 232: 111882. doi: 10.1016/j.engstruct.2021.111882
|
[25] |
SUN D L, CHEN D L, MI C L, et al., 2023. Evaluation of landslide susceptibility in the gentle hill-valley areas based on the interpretable random forest-recursive feature elimination model[J]. Journal of Geomechanics, 29(2): 202-219. (in Chinese with English abstract
|
[26] |
SUN N K, WANG Y, JI Z C, 2022. Anomaly detection method of electrical power consumption based on deep autoencoder[J]. Journal of System Simulation, 34(12): 2557-2565. (in Chinese with English abstract
|
[27] |
WANG C H, ZHAO Y J, GUO W, et al., 2022. Displacement prediction model of landslide based on ensemble empirical mode decomposition and support vector regression[J]. Acta Geodaetica et Cartographica Sinica, 51(10): 2196-2204. (in Chinese with English abstract
|
[28] |
WANG Y K, TANG H M, HUANG J S, et al., 2022. A comparative study of different machine learning methods for reservoir landslide displacement prediction[J]. Engineering Geology, 298: 106544. doi: 10.1016/j.enggeo.2022.106544
|
[29] |
WU J N, 2021. Research on landslide geological disaster prediction method based on deep learning[D]. Chengdu: University of Electronic Science and Technology of China. (in Chinese with English abstract
|
[30] |
WU Z H, HUANG N E, 2009. Ensemble empirical mode decomposition: a noise-assisted data analysis method[J]. Advances in Adaptive Data Analysis, 1(1): 1-41. doi: 10.1142/S1793536909000047
|
[31] |
XU S L, NIU R Q, 2018. Displacement prediction of Baijiabao landslide based on empirical mode decomposition and long short-term memory neural network in Three Gorges Area, China[J]. Computers & Geosciences, 111: 87-96.
|
[32] |
XU Z H, YANG X, SUN Q C, et al. , 2022. Step-type landslide deformation prediction based on multivariable self-optimizing dynamic neural network[J]. Metal Mine(3): 74-82. (in Chinese with English abstract
|
[33] |
YAN S J, TANG H M, XIANG W, 2007. Effect of rainfall on the stability of landslides[J]. Hydrogeology & Engineering Geology, 34(2): 33-36. (in Chinese with English abstract
|
[34] |
YANG B B, YIN K L, DU J, 2018. A model for predicting landslide displacement based on time series and long and short term memory neural network[J]. Chinese Journal of Rock Mechanics and Engineering, 37(10): 2334-2343. (in Chinese with English abstract
|
[35] |
YANG B B, YIN K L, LACASSE S, et al., 2019. Time series analysis and long short-term memory neural network to predict landslide displacement[J]. Landslides, 16(4): 677-694. doi: 10.1007/s10346-018-01127-x
|
[36] |
YAO W M, LI C D, ZUO Q J, et al., 2019. Spatiotemporal deformation characteristics and triggering factors of Baijiabao landslide in three gorges reservoir region, China[J]. Geomorphology, 343: 34-47. doi: 10.1016/j.geomorph.2019.06.024
|
[37] |
YUAN W, SUN R F, ZHONG H Y, et al., 2023. Research on comprehensive deformation prediction and monitoring and early warning methods for step-type like landslide[J]. Journal of Hydraulic Engineering, 54(4): 461-473. (in Chinese with English abstract
|
[38] |
ZHOU C, YIN K L, CAO Y, et al., 2018. A novel method for landslide displacement prediction by integrating advanced computational intelligence algorithms[J]. Scientific Reports, 8(1): 7287. doi: 10.1038/s41598-018-25567-6
|
[39] |
安冬,宋琨,仪政,等,2021. 一种基于EEMD-RFR的水库滑坡台阶状位移预测模型[J]. 山地学报,39(1):143-150.
|
[40] |
冯谕,曾怀恩,涂鹏飞,2023. 一种滑动检测算法下的滑坡位移时序分解方法[J]. 长江科学院院报,41(3):126-133,147.
|
[41] |
付文涛,胡丁梅,王守华,等,2020. 基于EEMD-GM-ELM模型的滑坡位移预测[C]//第十一届中国卫星导航年会论文集:S01卫星导航行业应用. 北京:中国卫星导航系统管理办公室学术交流中心:89-93.
|
[42] |
韩帅,刘明军,伍剑波,等,2022. 东南沿海台风暴雨型单体斜坡灾害风险评价:以泰顺仕阳北坡为例[J]. 地质力学学报,28(4):583-595. doi: 10.12090/j.issn.1006-6616.2021168
|
[43] |
康恩胜,赵泽熙,孟海东,2022. 基于EEMD-HW-PSO-ELM耦合模型的排土场边坡位移预测模型[J]. 黄金科学技术,30(4):594-602.
|
[44] |
李凌婧,姚鑫,周振凯,等,2022. Sentinel-1数据在西南山区水库变形斜坡InSAR监测中的适用性评价:以溪洛渡水库为例[J]. 地质力学学报,28(2):281-293.
|
[45] |
马 煜,余 斌,何元勋,等,2023. 降雨激发浅层滑坡发育特征与阈值研究:以江西省全南县大吉山“2019.6. 10”灾害为例[J]. 地质与勘探,59(5):1065-1073. doi: 10.12134/j.dzykt.2023.05.012
|
[46] |
南骁聪,刘俊峰,张永选,等,2023. 基于多源时间序列的滑坡位移动态预测[J]. 人民珠江,44(4):54-62.
|
[47] |
彭令,牛瑞卿,2011. 三峡库区白家包滑坡变形特征与影响因素分析[J]. 中国地质灾害与防治学报,22(4):1-7. doi: 10.3969/j.issn.1003-8035.2011.04.001
|
[48] |
曲建岭,王小飞,高峰,等,2014. 基于复数据经验模态分解的噪声辅助信号分解方法[J]. 物理学报,63(11):110201. doi: 10.7498/aps.63.110201
|
[49] |
尚敏,熊德兵,张惠强,等,2022. 基于时间序列与混合核函数SA-SVR的滑坡位移预测模型研究[J]. 工程地质学报,30(2):575-588.
|
[50] |
史学磊,韩旭东,杨秀元,等,2023. 三峡库区溪沟湾滑坡的诱发因素及前期降雨影响[J]. 地质力学学报,29(2):253-263. doi: 10.12090/j.issn.1006-6616.2022049
|
[51] |
孙德亮,陈丹璐,密长林,等,2023. 基于随机森林-特征递归消除模型的可解释性缓丘岭谷地貌滑坡易发性评价[J]. 地质力学学报,29(2):202-219. doi: 10.12090/j.issn.1006-6616.2022128
|
[52] |
孙宁可,王艳,纪志成,2022. 基于深度自编码器的电力能耗异常检测方法[J]. 系统仿真学报,34(12):2557-2565.
|
[53] |
王晨辉,赵贻玖,郭伟,等,2022. 滑坡位移EEMD-SVR预测模型[J]. 测绘学报,51(10):2196-2204. doi: 10.11947/j.AGCS.2022.20220291
|
[54] |
吴俊男,2021. 基于深度学习的滑坡地质灾害预测方法研究[D]. 成都:电子科技大学.
|
[55] |
徐志华,杨旭,孙钱程,等,2022. 基于多变量自优化动态神经网络的“阶跃型”滑坡变形预测[J]. 金属矿山(3):74-82.
|
[56] |
严绍军,唐辉明,项伟,2007. 降雨对滑坡稳定性影响过程分析[J]. 水文地质工程地质,34(2):33-36. doi: 10.3969/j.issn.1000-3665.2007.02.008
|
[57] |
杨背背,殷坤龙,杜娟,2018. 基于时间序列与长短时记忆网络的滑坡位移动态预测模型[J]. 岩石力学与工程学报,37(10):2334-2343.
|
[58] |
袁维,孙瑞峰,钟辉亚,等,2023. 阶跃型滑坡综合变形预测及监测预警方法研究[J]. 水利学报,54(4):461-473.
|