川东地区下寒武统龙王庙组储层地震响应特征
郗诚1, 和源1, 吕龑1, 李益2, 白晓亮1, 吴仕虎1, 李虹3
1.中国石油西南油气田公司勘探开发研究院
2.中国石油昆仑燃气四川分公司
3.中国石油西南油气田公司重庆气矿

作者简介:郗诚,1982年生,工程师,博士;主要从事油气综合地质、地球物理方法与储层预测等方面的研究工作。地址:(610041)四川省成都市天府大道北段12号。ORCID: 0000-0003-4167-8267。E-mail: xicheng@petrochina.com.cn

摘要

四川盆地中部下寒武统龙王庙组特大型气藏的发现展示了乐山—龙女寺古隆起外围巨大的天然气勘探潜力,川东地区在龙王庙组也发现了滩相储层,预示着该天然气勘探领域有可能向四川盆地东部扩展。为了弄清川东地区龙王庙组的天然气勘探潜力,有必要针对该区的储层识别和预测开展进一步的研究工作。为此,以岩石物理实验为基础,开展了基于波动方程的正演研究,归纳总结川东地区龙王庙组储层的地震反射特征,并与川中地区典型井龙王庙组反射特征进行对比,进而评价川东地区龙王庙组的天然气勘探潜力。研究结果表明:①川东地区龙王庙组储层发育在顶部及中部时具有较为明显的地震反射异常,而发育在中下部时则无明显的地震反射异常;②川东地区龙王庙组储层发育在顶部时龙王庙组顶界为复波反射,而发育在中上部时龙王庙组内部发育亮点;③川东地区龙王庙组与川中地区龙王庙组具有相似的地震反射特征。结论认为,五百梯西—巫山坎—正坝南地区及五宝场构造的东部应为下一步天然气勘探的首选方向,上述地区有望在川东地区龙王庙组取得天然气勘探突破。

关键词: 四川盆地东部; 早寒武世; 龙王庙期; 岩石物理实验; 波动方程; 地震相; 天然气勘探潜力; 有利勘探区
Seismic response characteristics of the Lower Cambrian Longwangmiao Formation reservoirs in the eastern Sichuan Basin
XI Cheng1, HE Yuan1, LYU Yan1, LI Yi2, BAI Xiaoliang1, WU Shihu1, LI Hong3
1. Exploration and Development Research Institute, PetroChina Southwest Oil & Gasfield Company, Chengdu, Sichuan 610051, China
2. Sichuan Branch, PetroChina Kunlun Gas Co., Ltd., Chengdu, Sichuan 610199, China
3. Chongqing Division of PetroChina Southwest Oil & Gasfield Company, Chongqing 402160, China
Abstract

The discovery of the supergiant gas reservoirs in the Lower Cambrian Longwangmiao Formation in the central Sichuan Basin reveals a great natural gas exploration potential in the periphery of Leshan-Longnüsi paleouplift. Shoal-facies reservoirs are also discovered in the Longwangmiao Formation of eastern Sichuan Basin, which indicates that the natural gas exploration area may expand to the eastern basin. In order to determine the natural gas exploration potential of Longwangmiao Formation there, it is necessary to make further studies on the reservoir identification and prediction of this area. Through petrophysical experiments, this paper carried out forward modeling researches based on wave equation. Then, the seismic reflection characteristics of the Longwangmiao Formation reservoirs in the eastern Sichuan Basin were summarized and compared with those in the typical well of the central Sichuan Basin. Finally, the natural gas exploration potential of Longwangmiao Formation in the eastern Sichuan Basin was evaluated. And the following research results were obtained. First, the Longwangmiao Formation reservoir in the eastern Sichuan Basin has an obvious seismic reflection anomaly when it is developed at the top and in the middle parts, but no obvious seismic reflection anomaly when it is developed in the middle-lower parts. Second, when the Longwangmiao Formation reservoir in the eastern Sichuan Basin is developed at the top, there is a complex wave reflection at the top boundary of Longwangmiao Formation; and when it is developed in the middle-upper parts, bright spots develop inside the Longwangmiao Formation. Third, the seismic reflection characteristics of the Longwangmiao Formation in the eastern Sichuan Basin are similar to those in the central Sichuan Basin. In conclusion, the Wubaitixi-Wushankan-Zhengbanan area and the eastern area of Wubaochang structure shall be taken as the first choice for natural gas exploration in the next step. The above areas are expected to make a breakthrough in natural gas exploration in Longwangmiao Formation in the eastern Sichuan Basin.

Keyword: Natural gas eastern Sichuan Basin; Earlier Cambrian; Longwangmiao Formation; Petrophysical experiment; Wave equation; Seismic facies; Natural gas exploration potential; Favorable exploration area
0 引言

近年来, 针对四川盆地震旦系— 寒武系持续开展了综合地质研究[1, 2, 3, 4, 5, 6], 积极实施油气风险勘探。2012年, 川中古隆起MX8井在下寒武统龙王庙组获得日产百万立方米高产工业气流; 2013年, 安岳气田龙王庙组气藏提交天然气探明储量4 403.83× 108 m3, 同时位于该古隆起斜坡区荷包场构造的HS1井在龙王庙组钻遇优质储层; 2013年, 在该古隆起斜坡带广安构造钻探的GT2井在上寒武统洗象池组钻遇厚层优质白云岩储层。上述勘探成果证实了四川盆地下古生界在川中古隆起及其外围具有良好的含气性和巨大的天然气勘探潜力[7, 8, 9]

震旦纪前, 受晋宁运动影响, 四川盆地东部达州— 开江地区古隆起已见雏形; 早寒武世, 受兴凯运动早期稳定抬升影响, 达州— 开江古隆起继承性演化。川中— 川东地区的格架大剖面显示, 达州— 开江古隆起与川中古隆起在隆起幅度上相当, 具有形成连片滩体发育的地质条件。在上述地质认识的指导下, 在川东地区部署了两口探井— — WT1井和LT1井, 它们均钻遇龙王庙组滩相细— 粉晶颗粒云岩储层, 测井解释储层孔隙度较低, 介于2%~4%, 为低能滩沉积[10], 表明川东地区龙王庙组也具有一定的油气勘探潜力。为了降低川东地区下古生界油气勘探的风险, 笔者基于前期地质研究的成果和新钻井资料, 以岩石物理实验为基础, 开展了基于波动方程的正演模拟, 总结了川东地区龙王庙组的地震反射特征, 并与川中高石梯— 磨溪地区典型井龙王庙组地震响应特征进行对比; 在此基础上, 完成了川东地区9 700 km2三维区的地震相刻画工作, 进而评价了川东地区龙王庙组天然气勘探的有利潜力区, 以期为该区下一步的油气勘探部署提供技术支撑。

1 地质背景

川东地区位于四川盆地东部(图1), 是我国最为典型的格挡式构造区, 地表主要由北东— 南西向条形背斜山和开阔的丘陵平坝相间组成, 背斜轴部多出露二叠系、三叠系, 开阔的向斜部位出露侏罗系砂泥岩地层。下构造层(中寒武统以下地层)变形相对较弱, 构造形态相对简单, 断层不发育。

图1 四川盆地区域地质图

川东地区寒武系地层自下而上分别为筇竹寺组(Є 1q)、沧浪铺组(Є 1c)、龙王庙组(Є 1l)、高台组(Є 2g)及洗象池组(Є 3x)。筇竹寺组岩性以黑灰色— 黑色页岩为主, 为下古生界重要的烃源岩层, 沧浪铺组以砂岩、粉砂岩及页岩为主, 龙王庙组以碳酸盐岩沉积为主, 在局部高能环境易形成高能滩体储层, 高台组以膏盐岩为主, 膏盐岩可以作为龙王庙组良好的区域盖层。因此, 川东地区龙王庙组具有形成大型油气田的地质条件。

2 基于岩石物理的地震响应特征
2.1 岩石物理特征

川东地区已完成WT1和LT1两口深层探井(图1)。WT1井龙王庙组储层物性相对较差, 储层位于龙王庙组下段, 发育细— 粉晶云岩, 储集空间主要以晶间溶孔为主, 测井解释储层厚度为12.0 m, 平均孔隙度为3.5%, 测井解释为水层。LT1井龙王庙组滩相发育, 为残余生屑含灰质云岩及残余粒屑云岩, 但储层厚度较薄, 测井解释储层厚度为5.4 m, 平均孔隙度为2.7%, 测井解释为差气层。WT1井及LT1井均未钻遇优质储层。

笔者在充分利用WT1井及LT1井的钻井资料的基础上, 补充了13个野外露头样品, 开展了岩石物理实验, 通过实验获得了纵波速度、横波速度、孔隙度、密度及渗透率等数据资料(表1), 将测得的实验数据作为龙王庙组地质模型建立的基础。

表1 川东地区龙王庙组野外柱塞样岩石物理实验数据表

实验数据表明, 川东地区龙王庙组岩石孔隙度普遍较低, 渗透率也较低, 仅6号及9号样品渗透率高, 岩样观察发现这两块样品均发育不同程度的裂缝, 渗透率得到了明显改善。川东地区龙王庙组基质孔隙较低, 裂缝及岩溶作用对储层的改造较为重要[7, 11, 12]

图2为排除6号及9号样品后得到的岩石物性交汇图。图2-a、b表明, 密度与纵/横波速度具有较好的正相关性, 纵/横波速度随着密度的增大而增大; 图2-c表明, 孔隙度与纵波速度呈负相关关系, 随着孔隙度的增大, 纵波速度减小, 两者存在一定的线性关系, 为后续不同孔隙度的地质模型的正演提供了数据支撑; 图2-d反映出孔隙度与密度呈负相关关系。

图2 川东地区龙王庙组岩石物理分析图

2.2 地质模型与地震响应特征

地震波正演模拟是研究地球介质中地震波传播的运动学和动力学特征的重要手段, 也是地震资料偏移成像的基础[13, 14, 15]。地震波正演模拟常常被用于评价地震处理的效果、检验地震解释的结论, 甚至直接被用于含油气性检测。笔者在传统褶积模型的基础上, 采用了基于波动方程的正演模拟, 使得结果更加符合地下真实的波场特征[16, 17]

在获取岩石物理实验数据的基础上, 建立了地层格架模型, 速度背景主要由WT1井的实钻曲线得到(图3)。WT1井高台组以膏质云岩及云质膏盐岩为主, 纯盐岩不发育, 以低速层发育为主, 平均速度为5 747 m/s; 龙王庙组主要发育白云岩, 平均速度为6 850 m/s; 沧浪铺组地层较为复杂, 为了真实还原地下地层特征, 将沧浪铺组划分为四段, 平均速度由上至下分别为5 263 m/s、6 290 m/s、5 128 m/s、6 490 m/s; 筇竹寺组作为龙王庙组的烃源岩, 主要以低速为主, 平均速度为5 320 m/s。

图3 川东地区WT1井龙王庙组地层格架模型图

根据WT1井测井曲线建立的正演地质模型(图4-a), 分别在龙王庙组的顶部、中部及底部建立孔隙度为2%的储层(楔状储层厚度介于0~20 m), 正演采用的地震主频为25 Hz, 与川东地区深层实际地震的主频保持一致, 真实还原波场的反射。根据岩石物理实验交汇得到的图板, 储层孔隙度为2%对应的纵波速度及密度分别为5 500 m/s、2.81 g/cm3, 以此为储层的填充速度及密度, 正演模拟结果(图4-b)表明当储层发育龙王庙组顶部时, 龙王庙组顶界波峰出现一定的下拉; 储层发育在龙王庙组中部时, 龙王庙组顶界减弱, 呈复波反射特征, 且储层底部为波峰反射; 当储层发育在龙王庙组中下部时, 地震反射无明显响应。

图4 孔隙度为2%的储层在不同发育位置的速度模型及正演模拟结果图

同时, 为了论证物性变化对反射特征的影响, 设计了不同孔隙度的地质模型, 储层均发育在龙王庙组顶部, 储层孔隙度分别为1%、2%、3%, 由岩石物理实验图板统计对应的纵波速度分别为6 200 m/s、5 600 m/s及5 000 m/s(图5-a), 密度分别为2.84 g/cm3、2.81 g/cm3、2.78 g/cm3, 将上述纵波速度及密度作为储层的约束参数, 通过波动方程正演得到相应的波场记录(图5-b)。正演模拟结果表明, 随着孔隙度的变化, 地震反射也出现相应的差异。储层孔隙度为1%时, 地震反射没有明显的变化; 储层孔隙度为2%时, 龙王庙组顶界出现下拉; 储层孔隙度为3%时, 龙王庙组顶界出现明显减弱, 储层下部发育强亮点。

图5 不同孔隙度储层地质模型及正演模拟结果图

通过两个地质模型的正演模拟, 得到了两个认识:①当储层发育在龙王庙组顶部及中部时, 利用主频为25 Hz的地震可以进行识别, 而当储层发育在底部时, 储层没有明显响应。WT1井的储层发育在中下部, 实际地震剖面没有出现明显响应, 证实了正演结果的可靠性; ②随着孔隙度的增加(储层物性变好), 储层引起的亮点反射逐渐增强。

3 川东地区龙王庙组地震反射特征

根据波动方程正演模拟取得的认识, 对川东地区9 700 Km2三维地震开展龙王庙组地震反射识别, 定性解释了龙王庙组顶界复波、层间亮点(顶界弱反射或顶界空白反射)、顶界强反射— 龙王庙组内部空白反射等3种地震相。并与高石梯— 磨溪地区典型井的地震反射进行了对比, 认为川东地区龙王庙组与高石梯— 磨溪地区龙王庙组储层具有相似的反射特征, 可以有效地对川东地区龙王庙组有利储层进行定性预测, 寻找川东地区龙王庙组下一步的油气勘探潜力区。

图6为川东地区五百梯西三维工区line605线地震剖面图, 图中红色虚线为解释的亮点, 龙王庙组顶界呈弱波峰反射, 正演模拟结果表明此类地震反射特征储层较为发育。图7为川东地区铁山— 龙门连片三维trace2400线地震剖面图, 该剖面上解释有复波反射、亮点反射(龙王庙组顶界反射减弱), 正演模拟结果表明该线在龙王庙组中上部储层应较发育。

图6 五百梯构造西line605线地震剖面图

图7 铁山— 龙门连片三维trace2400线地震剖面图

表2为川东地区与川中地区龙王庙组地震反射特征对比表, 归纳解释了4类反射特征。其中Ⅳ 类反射特征主要以龙王庙组顶界强峰、层间空白反射为主, WT1井即为此种反射, 该井测井解释储层较薄, 测井解释为水层, Ⅳ 类反射特征类似于高石梯— 磨溪地区的GS2井, GS2井储层仅有4.5 m, 储层欠发育, 该类地震反射特征主要集中分布在川东地区檀木场区块及龙岗区块; Ⅲ 类反射特征主要以龙王庙组顶界复波反射为主, 该类反射特征与MX23井类似, MX23井测井解释储层厚度为20.0 m, 测试日产气量为110.8× 104 m3, 该类反射特征主要集中分布在川东地区五宝场区块及铁山— 龙门地区; Ⅱ 类反射特征为龙王庙组顶界弱波峰反射、层间发育亮点, 类似于川中地区的MX13井, MX13井储层厚度为41.7 m, 测试日产气量为128.84× 104 m3, 该类反射特征主要集中分布在正坝南及大猫坪西区块; Ⅰ 类反射特征以龙王庙组顶界空白反射, 层间亮点发育为主, 类似于川中地区MX204井, MX204井储层厚度为54.0 m, 测试日产气量为115.62× 104 m3, 该类反射特征主要分布在五百梯西地区。

表2 川东地区与川中地区龙王庙组地震反射特征对比表
4 川东地区龙王庙组勘探领域

图8为川东地区龙王庙组地震相刻画平面图。图8中蓝色对应Ⅳ 类反射特征, 主要以顶界强波峰, 内部空白反射; 黄色对应Ⅲ 类反射特征, 主要以复波发育为主; 橙色对应Ⅱ 类及Ⅰ 类反射特征, 以亮点发育为主。需要说明的是, 图8中的黑色区域, 由于川东地区整体呈北东— 南西向的高陡格挡式构造格局, 高陡区构造成像较差, 在深层下古生界基本没有成像, 地震反射较为杂乱, 未作地震相刻画。

图8 川东地区三维工区龙王庙组地震相刻画平面图

从图8可以看出, 川东地区龙王庙组地震反射特征主要以Ⅳ 类反射为主, 储层不发育(或储层纵向厚度较小)。Ⅲ 类反射特征主要发育在在巫山坎、正坝南、大猫坪西地区。Ⅰ 类及Ⅱ 类反射特征主要发育在五百梯西及五宝场地区。通过本次地震相的刻画, 认为下一步川东龙王庙组的勘探潜力区应首先优选五百梯西— 巫山坎— 正坝南地区及五宝场区块东部, 以期在川东地区下古生界龙王庙组获得一定的突破。

5 结论

1)储层发育在龙王庙组顶部及中部时, 25 Hz频率的地震可以识别, 当储层发育在龙王庙组中下部时, 储层没有明显响应; 随着孔隙度的增大, 储层引起的亮点反射逐渐增强。

2)川东地区与川中地区的龙王庙组储层地震响应特征对比结果表明, 川东地区龙王庙组储层的地震响应特征与川中地区高石梯— 磨溪地区相同。

3)通过本次地震相的研究成果, 认为应首选五百梯西— 巫山坎— 正坝南地区作为川东地区龙王庙组储层发育的有利潜力区, 其次为五宝场区块东部地区。

编 辑 陈古明

参考文献
[1] 杨跃明, 杨雨, 杨光, . 安岳气田震旦系、寒武系气藏成藏条件及勘探开发关键技术[J]. 石油学报, 2019, 40(4): 493-508.
YANG Yueming, YANG Yu, YANG Guang, et al. Gas accumulation conditions and key exploration & development technologies of Sinian and Cambrian gas reservoirs in Anyue gas field[J]. Acta Petrolei Sinica, 2019, 40(4): 493-508. [本文引用:]
[2] 徐春春, 沈平, 杨跃明, . 乐山—龙女寺古隆起震旦系—下寒武统龙王庙组天然气成藏条件与富集规律[J]. 天然气工业, 2014, 34(3): 1-7.
XU Chunchun, SHEN Ping, YANG Yueming, et al. Accumulation conditions and enrichment patterns of natural gas in the Lower Cambrian Longwangmiao Fm reservoirs of the Leshan-Longnüsi Paleohigh, Sichuan Basin[J]. Natural Gas Industry, 2014, 34(3): 1-7. [本文引用:]
[3] 邹才能, 杜金虎, 徐春春, . 四川盆地震旦系—寒武系特大型气田形成分布、资源潜力及勘探发现[J]. 石油勘探与开发, 2014, 41(3): 278-293.
ZOU Caineng, DU Jinhu, XU Chunchun, et al. Formation, distribution, resource potential and discovery of the Sinian-Cambrian giant gas field, Sichuan Basin, SW China[J]. Petroleum Exploration and Development, 2014, 41(3): 278-293. [本文引用:]
[4] 李辉, 朱祥, 李毕松, . 四川盆地北部震旦系天然气勘探潜力分析[J]. 特种油气藏, 2018, 25(1): 52-57.
LI Hui, ZHU Xiang, LI Bisong, et al. Sinian natural gas exploration potential analysis in the northern Sichuan Basin[J]. Special Oil & Gas Reservoirs, 2018, 25(1): 52-57. [本文引用:]
[5] 魏国齐, 杨威, 谢武仁, . 四川盆地震旦系—寒武系天然气成藏模式与勘探领域[J]. 石油学报, 2018, 39(12): 1317-1327.
WEI Guoqi, YANG Wei, XIE Wuren, et al. Accumulation modes and exploration domains of Sinian-Cambrian natural gas in Sichuan Basin[J]. Acta Petrolei Sinica, 2018, 39(12): 1317-1327. [本文引用:]
[6] 谷志东, 殷积峰, 袁苗, . 四川盆地东部深层盐下震旦系—寒武系天然气成藏条件与勘探方向[J]. 石油勘探与开发, 2015, 42(2): 137-149.
GU Zhidong, YIN Jifeng, YUAN Miao, et al. Accumulation conditions and exploration directions of natural gas in deep subsalt Sinian-Cambrian System in the eastern Sichuan Basin, SW China[J]. Petroleum Exploration and Development, 2015, 42(2): 137-149. [本文引用:]
[7] 肖富森, 冉崎, 唐玉林, . 乐山—龙女寺古隆起深层海相碳酸盐岩地震勘探关键技术及其应用[J]. 天然气工业, 2014, 34(3): 67-73.
XIAO Fusen, RAN Qi, TANG Yulin, et al. Key technologies and their application to seismic exploration of the deep marine carbonate reservoirs in the Leshan-Longnüsi paleouplift, Sichuan Basin[J]. Natural Gas Industry, 2014, 34(3): 67-73. [本文引用:]
[8] 周路, 张欣吉, 钟克修, . 四川盆地东部五百梯地区构造变形特征[J]. 新疆石油地质, 2018, 39(3): 251-256.
ZHOU Lu, ZHANG Xinji, ZHONG Kexiu, et al. Structure deformation characteristics of Wubaiti area in the eastern Sichuan Basin[J]. Xinjiang Petroleum Geology, 2018, 39(3): 251-256. [本文引用:]
[9] 冯伟明, 谢渊, 刘建清, . 上扬子下寒武统龙王庙组沉积模式与油气勘探方向[J]. 地质科技情报, 2014, 33(3): 106-111.
FENG Weiming, XIE Yuan, LIU Jianqing, et al. Sedimentary model and hydrocarbon exploration targets of the Lower Cambrian Longwangmiao Formation in the Upper Yangtze area[J]. Geological Science and Technology Information, 2014, 33(3): 106-111. [本文引用:]
[10] 杜金虎, 张宝民, 汪泽成, . 四川盆地下寒武统龙王庙组碳酸盐缓坡双颗粒滩沉积模式及储层成因[J]. 天然气工业, 2016, 36(6): 1-10.
DU Jinhu, ZHANG Baomin, WANG Zecheng, et al. Sedimentary model and reservoir genesis of dual grain banks at the Lower Cambrian Longwangmiao Fm carbonate ramp in the Sichuan Basin[J]. Natural Gas Industry, 2016, 36(6): 1-10. [本文引用:]
[11] 李伟, 余华琪, 邓鸿斌. 四川盆地中南部寒武系地层划分对比与沉积演化特征[J]. 石油勘探与开发, 2012, 39(6): 681-690.
LI Wei, YU Huaqi, DENG Hongbin. Stratigraphic division and correlation and sedimentary characteristics of the Cambrian in central-southern Sichuan Basin[J]. Petroleum Exploration and Development, 2012, 39(6): 681-690. [本文引用:]
[12] 周进高, 徐春春, 姚根顺, . 四川盆地下寒武统龙王庙组储集层形成与演化[J]. 石油勘探与开发, 2015, 42(2): 158-166.
ZHOU Jingao, XU Chunchun, YAO Genshun, et al. Genesis and evolution of Lower Cambrian Longwangmiao Formation reservoirs, Sichuan Basin, SW China[J]. Petroleum Exploration and Development, 2015, 42(2): 158-166. [本文引用:]
[13] 陈可洋. 三维随机建模方法及其波场模拟分析[J]. 勘探地球物理进展, 2009, 32(5): 315-320.
CHEN Keyang. 3-D rand om modeling scheme and wavefield simulation analysis[J]. Progress in Exploration Geophysics, 2009, 32(5): 315-320. [本文引用:]
[14] 李振春, 张华, 刘庆敏, . 弹性波交错网格高阶有限差分法波场分离数值模拟[J]. 石油地球物理勘探, 2007, 42(5): 510-515.
LI Zhenchun, ZHANG Hua, LIU Qingmin, et al. Numeric simulation of elastic wavefield separation by staggering grid high-order finite-difference algorithm[J]. Oil Geophysical Prospecting, 2007, 42(5): 510-515. [本文引用:]
[15] 刘洪, 刘国峰, 武威, . 多维波动方程逆散射的基础理论研究[J]. 石油物探, 2007, 46(6): 569-581.
LIU Hong, LIU Guofeng, WU Wei, et al. Study into fundamental theory of inverse scattering of multi-dimensional wave equation[J]. Geophysical Prospecting for Petroleum, 2007, 46(6): 569-581. [本文引用:]
[16] 奚先, 姚姚. 二维粘弹性随机介质中的波场特征分析[J]. 地球物理学进展, 2004, 19(3): 608-615.
XI Xian, YAO Yao. The analysis of the wave field characteristics in 2-D viscoelastic rand om medium[J]. Progress in Geophysics, 2004, 19(3): 608-615. [本文引用:]
[17] 董良国, 马在田, 曹景忠. 一阶弹性波方程交错网格高阶差分解法稳定性研究[J]. 地球物理学报, 2000, 43(6): 856-864.
DONG Liangguo, MA Zaitian, CAO Jingzhong. A study on stability of the staggered-grid high-order difference method of first-order elastic wave equation[J]. Chinese Journal of Geophysics, 2000, 43(6): 856-864. [本文引用:]