In-Situ Stress State in the Main Borehole of the Chinese Continental Scientific Drilling
-
摘要: 用钻孔崩落法确定了中国大陆科学钻探主钻孔5 047 m深度以上的现今地应力状态.由钻孔声波成像测井资料发现, 科学钻主钻孔在1 200 m深度以下出现了钻孔崩落现象.我们从1 216~5 047 m的深度范围内采集了143个钻孔成像测井图象资料, 对钻孔崩落椭圆长轴方位进行了统计, 结果表明崩落椭圆长轴平均方位为319.5°±3.5°, 最大水平主应力方位平均为49.5°±3.5°.利用崩落形状要素(崩落深度和崩落宽度) 以及岩石的内聚力和内摩擦角, 估算了1 269 m至5 047 m范围内52个深度上的最大和最小水平主应力的大小.结果表明, 在浅处1 216 m深度, 最大水平主应力为42 MPa, 最小水平主应力为30.3 MPa; 在深处5 000 mm深度, 最大水平主应力为160.5 MPa, 最小水平主应力为120 MPa; 地应力随深度近于线性增加.据岩石密度测井资料计算了各个深度上静负载应力.3个主应力的大小和方向反映出科学钻主孔位置的应力场处于走滑应力状态, 与临近地区地震震源机制解和其他方法得到的应力场一致.利用声发射法对岩心试件进行了声发射测量, 得到了最大水平主应力幅值, 并与崩落法测量结果进行了对比, 两者十分一致.Abstract: The breakouts began to occur under the depth of 1 200 m in the main borehole of Chinese Continental Scientific Drilling (CCSD), 143 breakout images are collected from acoustic borehole televiewer data between 1 216 m and 5 047 m. The average orientation of breakouts is 319.5°±3.5°. The average orientation of maximum horizontal stress is 49.5°±3.5°. Using breakout parameters measured from borehole breakouts (depth and width of breakout) and rock cohesion and inner friction angle determined from triaxial rock compression and deformation tests, the magnitudes of principal stress are calculated at 52 different depths from 1 269 m to 5 047 m. Overburden stress is calculated using the density logging data. The magnitude of maximum principal stress is determined using acoustic emissions. A comparison of measurement results of acoustic emission and borehole breakouts is carried out, and the results from the two methods are in good agreement. According to the orientations and magnitudes of three principal stresses, stress state (SH > Sv > Sh) in the area of CCSD indicates the strike-slip stress regime and is coincident with the strike-slip stress field observed from earthquake focal plane mechanisms.
-
0. 引言
中国大陆科学钻探的主孔位于江苏省东海县毛北村, 构造位置上处于郯庐断裂附近, 在苏鲁-大别高压超高压变质带内.了解钻孔地区区域地应力状态是该区地壳构造动力学研究不可缺少的内容.声波成像测井资料显示, 钻孔在1 200 m深度以下出现了钻孔崩落现象.理论分析和室内试验已经证明了这种崩落现象是地应力作用的结果.钻孔崩落出现在钻孔壁挤压应力集中处, 钻孔壁的压剪破裂导致孔壁剥落, 形成椭圆形钻孔.崩落椭圆的长轴平均方位垂直于最大水平主应力方向(Bell and Gough, 1979; Gough and Bell, 1981; Plumb, 1982; Blumling et al., 1983), 这已被许多钻井观测资料所证实(Gough and Bell, 1981; Hickman et al., 1985; Vernik and Zoback, 1992), 因此, 钻孔崩落可作为地应力主方向的指示器. Zoback et al. (1985)进一步发展了钻孔崩落理论, 指出孔壁附近出现压应力集中区并形成压剪破裂, 造成崩落, 崩落参数与地应力存在定量关系. Haimson and Herrick (1986)进行了室内实验, 结果表明, 钻孔崩落的位置和形状与Zoback et al. (1985) 的理论模型非常一致.证实了崩落椭圆的长轴平均方向垂直于最大水平主应力方向. Zoback et al. (1985)的理论模型不仅能利用钻井崩落来确定主应力方向, 而且还可确定水平主应力的大小.
钻孔崩落现象已被广泛用来确定地应力状态(Hickman et al., 1985; Plumb and Hickman, 1985; Zoback and Zoback, 1989; Vernik and Zoback, 1992; 王连捷等, 1993; 俞言祥和许忠淮, 1994, 1996; Huber et al., 1997; 崔军文等, 2004; 李明武等, 2005).
本文的目的就是利用中国大陆科学钻探的主孔1 216~5 047 m的钻孔崩落资料来确定最大和最小水平应力方向和大小, 从而评价主孔现今所处的构造应力环境, 为该区地壳构造动力学研究提供定量化信息.
1. 钻孔崩落资料的获取
以往多用四臂倾角测井记录获取钻孔崩落的资料, 但这种方法只能得到崩落椭圆的长轴和短轴的长度和方位, 不能得到崩落断面的形状.由声波成像测井记录可以得到钻孔崩落断面的详细资料.图 1为科学钻2 170~2 175 m深度的声波成像测井记录.图 1a中黑色竖直条带为崩落区, 正弦形黑色条带为与钻井斜交的破裂面, 上部2 170 m附近的多条黑色条带是小裂缝.图 1b是钻井崩落后的横断面, 由图可以看出崩落区的几何形状、崩落的宽度和深度.图 1b中的崩落区与图 1a中的黑色竖直条带相对应.根据声波成像测井记录, 每隔一定的深度取一个横断面, 量测出崩落区的几何形状参数——崩落的宽度和深度, 进而计算应力场的方向和大小.
2. 钻孔崩落的力学分析
假设远场最大和最小水平主应力分别为σH和σh, 井内泥浆流体对井壁的压力为Pw, 岩层内孔隙压力为P, 则井孔周围的应力状态表达式如下(Jaegr, 1961) :
(1) (2) (3) 其中, σr为径向应力, σθ为周向应力, τrθ为切向剪应力, ΔP为Pw-P, a为井孔半径, θ为由最大水平主应力σH方位向量测方位的夹角, 逆时针为正.
由式(1)、(2)、(3) 可知, 当r=a时, 即孔壁上,
(4) (5) (6) 由方程(4)、(5)、(6) 可知, 引起钻孔应力集中, 诱发孔壁破坏, 关键是周向应力σθ,
(7) (8) 由式(4)、(7)、(8) 可见, 当角度θ=π/2及3π/2时, 即在M点、N点(图 2a), 周向应力σθ值最大, 径向应力σr与周向应力σθ的差值σθ-σr也最大.可见, 在θ=π/2及3π/2的方向, 形成压应力及剪应力集中区(图 2b).当应力达到岩石的强度时, 产生孔壁崩落.崩落后形成的椭圆钻孔, 其长轴方向与最小水平主应力方向平行.
实际情况中, 主应力比值σH/σh为大于1, 小于3, 而ΔP很小.由式(7)、(8) 可见, 这时钻孔附近只有压应力, 没有张应力, 因而, 只有剪破坏, 没有张裂隙.
根据纳维-库仑破坏准则, 在给定的地应力条件下, 井壁发生破坏崩落的条件由下式确定:
(9) 当上式得到满足时, 即发生崩落, 即上式的右边如果小于左边, 井孔是稳定的, 若右边大于或等于左边, 则发生崩落.
将式(1)-式(3) 代入式(9) 并令θ=θb, r=a和r=rb, θ=π/2, 得到:
(10) (11) 式中:
当σθ-σr > 0时, a2、b2、c2和d2取正号, e、f取负号; σθ-σr < 0时, a2、b2、c2和d2取负号, e、f取正号.
由
, 可以解出2个水平主应力σH和σh (Zoback et al., 1985) :(12) (13) 3. 地应力的确定
一点的应力状态可由主应力的方向和大小来决定.确定一点的应力状态就意味着确定主应力的方向和大小.
3.1 主应力方向的确定
在钻孔近于铅直的情况下, 由于井壁崩落椭圆长轴方向平行于最小水平主应力方向, 因此, 可以利用成像测井资料得到的崩落椭圆长轴方向, 经统计平均后, 得到最大水平主应力方向.
3.2 主应力大小的确定
利用泥浆密度资料, 求出井孔压力与地层压力之差ΔP=Pw-P.通过岩石力学试验获得原地岩石的强度参数, 即岩石内聚力和内摩擦角, 在崩落椭圆断面上, 量取崩落深度(rb)和崩落角(θb), 根据这些数据, 可以利用式(12)、(13) 确定现场地应力的大小.
量取崩落深度(rb) 时, 要注意后续崩落扩展对崩落深度的影响.在计算地应力时, 要求使用初始崩落深度.实际情况中, 初始崩落后, 崩落经常继续扩展.崩落扩展对崩落角(θb) 影响较小, 但对崩落深度(rb) 的影响较大.因而, 量测初始崩落深度往往比较困难.为此, (1) 选择崩落扩展程度小的断面进行计算; (2) 崩落区两侧受后续崩落影响小, 可利用两侧的轮廓线大致恢复初始崩落区的形态, 确定崩落深度.
4. 大陆科学钻主钻孔现今地应力状态
4.1 最大主应力方向的确定
我们在1 216~5 047 m的范围内, 采集了143个声波成像测井横断面图像资料, 取得了崩落椭圆长轴方位角和相应的最大主应力方位.对这143个数据进行方向随深度分布的统计(图 3).
由图 3可以看出崩落椭圆长轴方位随深度略有变化.在2 600 m以上, 长轴方位随深度由北向西偏转约10°, 3 600 m深度以下, 方位平稳, 并少许向回偏转约3°~4°.最大水平主应力方位随深度亦相应偏转.总体来看, 方位随深度变化不大.在5 000 m深度范围内, 崩落椭圆长轴方位平均为319.5°±3.5°, 最大水平主应力方位平均为49.5°±3.5° (图 4), 可信度为90%.即90%的数据, 其角度变化范围为±3.5°.仅有10%的情况, 大于±3.5°.
由于2 600~3 600 m的成像测井记录质量较差, 未能得到该深度范围的崩落资料.
4.2 主应力大小的估算
4.2.1 岩石力学实验
估计主应力的大小时, 需要岩石强度资料.为此, 在主孔采集了有代表性的12个片麻岩岩石样品、12个角闪岩样品和12个榴辉岩样品进行了岩石力学试验.
根据岩石力学实验结果, 片麻岩岩石内聚力平均为19.20 MPa, 角闪岩为22.39 MPa, 榴辉岩为24.77 MPa, 相应地, 岩石内摩擦角分别为51.5°、48.6°和51.9°, 摩擦系数μ分别为1.258、1.134和1.275, 平均单轴破坏强度分别为110.03 MPa、118.56 MPa和143.47 MPa.
4.2.2 主应力大小的估算
我们对铅直应力大小进行了估算.由于井孔近乎直立, 其中一个主应力为铅直的.铅直应力可用上覆岩层载荷来计算.计算公式如下:
其中: ρ为岩石的平均密度, g为重力加速度, H为深度.
密度由测井资料获得, 取密度ρ为2.8 g/cm3, 得到的铅直应力如图 5.
4.2.3 水平主应力大小的估算
利用声波成像测井图像资料, 在1 269~5 047 m深度上, 测量了52钻孔崩落断面的初始崩落深度与崩落张开角.根据前述公式计算了最大(SH)和最小(Sh) 水平主应力的大小如表 1和图 5.图 5还给出了地应力随深度变化的回归曲线, 其表达式为:
表 1 最大和最小水平主应力的大小Table Supplementary Table Magnitudes of the maximum horizontal principle stress and minimum horizontal principle stress由图 5可以看出地应力随深度增加, 最大水平主应力在1 269 m深度为44.8 MPa, 在5 000 m深度为160.5 MPa, 最小水平主应力相应为31.6 MPa和117 MPa.3个主应力的大小关系为: SH > Sv > Sh, 表明主孔地应力为走滑应力状态.
4.2.4 钻孔崩落法与声发射法测量结果的比较
对中国大陆科学钻探主孔的岩心进行了声发射测量(王连捷等, 2005).沿钻孔5个深度即301 m, 598 m, 1 003 m, 1 258 m, 1 531 m分别取样, 每个点10个岩样, 进行声发射测量(丁原辰和张大伦, 1991; 丁原辰, 1992).将结果列入表 2, 同时放入图 5.由图 5可以看出, 2种方法所得结果非常一致.声发射法所得测值基本上落在钻孔崩落法测值的趋势线上.两者可以相互验证, 说明测试结果是可靠的.声发射法可以补充浅部资料, 声发射法与钻孔崩落法相结合使测量数据更加充分和完整.
表 2 声发射法测量结果Table Supplementary Table Result of acoustic emission measurement5. 大陆科学钻主孔附近地应力格局
将中国大陆科学钻主孔附近的震源机制解和其他方法所得结果(Zoback, 1992) 与大陆科学钻主孔井壁崩落确定的地应力数据一起, 编制了现今最大主应力方位分布图(图 6).可以看出科学钻主孔附近地区, 现今最大水平主应力方位在北东和北东东之间.震源机制解和钻孔崩落数据皆显示这一地区的区域地应力主要为走滑应力状态.
图 6 中国大陆科学钻探(CCSD) 主孔位置及及钻孔崩落和震源机制解显示的应力格局示意图(钻孔附近地区的地应力资料源于Zoback et al., 1992)BO为钻孔崩落指示的最大主应力方向; SS、TF、NF为震源机制解指示的最大主应力方向和应力格局; SS为走滑断层作用; TF为逆冲断层作用; NF正断层作用Fig. 6. Sketch map for the location of Chinese Continental Scientific Drill (CCSD) and stress regime indicated by breakouts and earthquake focal mechanisms周翠英等(1994)利用1969年至1991年的地震活动资料对大陆科学钻主孔所在位置的苏鲁-南黄海广大地区的应力场的动态变化进行了研究, 发现该区应力场总体特征为北东东向, P轴倾角小于10°. 华爱军等(1997)利用地震活动参数对沂沭断裂附近应力场动态变化进行了分析, 南段震源机制解P轴方位1978—1984年稳定在N72°E, 1985年后, 变为N50°E. 谢智等(2002)汇编了河南及邻区(山东、河北、安徽) 1965—2000年共52次ML≥3.5的地震震源机制解, 表明P轴方位多为北东到北东东向, 地震断层主要为走滑类型.
郯庐断裂带在燕山早期为左行平移运动, 在新生带晚近时期的活动则表现为右行扭动(王小风等, 2000). 刘光勋等(1989)、汪一鹏等(1989)也指出郯庐断裂的现今活动为水平右行扭动.中国大陆科学钻主孔地应力状态及附近地区地应力方向与郯庐断裂带晚近时期以来的扭动方向是协调的.
综上所述, 科学钻主孔附近地区地震震源机制解和其他方法得到的结果与钻孔崩落的结果基本一致, 与主要断裂带的活动性质一致.本区现今地应力格局主要为走滑状态, 最大水平主应力方向在北东和北东东之间.
6. 结论
(1) Zoback and Zoback (1989)对应力方向数据建立了可靠的定量指标.该指标规定, 对于A级质量的钻孔崩落数据, 要求崩落出现在钻孔两侧近似对称(相隔180°) 位置, 参加方位统计的数量大于10, 并贯穿井深300 m以上, 标准偏差≤12°.我们统计的钻孔崩落数量为143个, 选取的崩落资料满足两侧近似对称(相隔180°) 的要求, 覆盖深度由1 216 m至5 047 m, 贯穿井深3 800多m, 标准偏差≤3.5°.因此我们统计的最大主应力的方向是可信的.
(2) 地应力大小随深度增加, 在浅部1 269 m, 最大水平主应力为44.8 MPa, 最小主应力为31.6 MPa, 在深部5 047 m, 最大水平主应力为160.5 MPa, 最小主应力为117 MPa.最大水平主应力SH和最小主应力Sh大小随深度的增加率为:
主应力方向为N49.5±3.5°E, 可信度90%.
(3) 声发射法与钻孔崩落法测量结果的对比表明两者十分一致, 说明地应力的大小及其随深度的变化是可靠的.
(4) 综合分析主孔附近地区地震震源机制解等其他方法所得地应力资料, 可以看出本区区域现今地应力主方向为北东到北东东向压应力, 为走滑型应力状态格局.
(5) 所得地应力资料可作为构造动力学和地质灾害研究的基础资料.
致谢: 研究过程中得到了中国大陆科学钻探工程中心牛一群研究员、潘和平教授、胜利油田测井公司吴海燕高级工程师的大力支持和帮助, 并提供了相关资料, 在此表示衷心的感谢.在样品采集和加工过程中, 构造岩石编录组、物性研究组和现场蔡慈高级工程师、孙立文工程师和王尧工程师给予的大力支持和帮助, 在此也表示诚挚的感谢! -
图 6 中国大陆科学钻探(CCSD) 主孔位置及及钻孔崩落和震源机制解显示的应力格局示意图(钻孔附近地区的地应力资料源于Zoback et al., 1992)
BO为钻孔崩落指示的最大主应力方向; SS、TF、NF为震源机制解指示的最大主应力方向和应力格局; SS为走滑断层作用; TF为逆冲断层作用; NF正断层作用
Fig. 6. Sketch map for the location of Chinese Continental Scientific Drill (CCSD) and stress regime indicated by breakouts and earthquake focal mechanisms
表 1 最大和最小水平主应力的大小
Table 1. Magnitudes of the maximum horizontal principle stress and minimum horizontal principle stress
表 2 声发射法测量结果
Table 2. Result of acoustic emission measurement
-
[1] Bell, J.S., Gough, D.I., 1979. Northeast-southwest compressive stress in Alberta— Evidence from oil wells. Earth and Planetary Sciences Letters, 45: 475-482. doi: 10.1016/0012-821X(79)90146-8 [2] Blumling, P., Fuchs, K., Schneider, T., 1983. Orientation of the stress field from breakouts in a crystalline well in a seismic active area. Phys. Earth Planet. Inter. , 33: 250-254. doi: 10.1016/0031-9201(83)90042-0 [3] Cui, J. W., Wang, L. J., Li, P. W., et al., 2004. Tectonic stress field of brittle deformation within 2 000 m of the main borehole of Chinese Continental Scientific Drilling. Acta Petrologica Sinica, 20(1): 73-80(in Chinese with English abstract). [4] Ding, Y. C., 1992. Rock stress measureament by the AE method and its application in mine. Coal Engineer, (4): 50-56(in Chinese with English abstract). [5] Ding, Y.C., Zhang, D. L., 1991. Application of the incomplete erasion phenomenon in acoustic emission activities to the measurement of geostresses. Chinese Journal of Rock Mechanics and Engineering, 10(4): 313-326(in Chinese with English abstract). [6] Gough, D.I., Bell, J.S., 1981. Stress orientations from oil well fractures in Alberta and Texas. Can. J. Earth Sci. , 18: 638-645. doi: 10.1139/e81-056 [7] Haimson, B.C., Herrick, C.G., 1986. Borehole breakouts A new tool for estimation in situ stress. Proceedings of the international sympossium on rock stress and rock stress measurements, Stockholm, 1-3, Sept. . [8] Hickman, S.H., Healy, J.H., Zoback, M.D., 1985. In situ stress, natural fracture distribution, and borehole elongation in the Auburn geothermal well, Auburn, New York. Journal of Geophysics Research, 90: 5497-5512. doi: 10.1029/JB090iB07p05497 [9] Hua, A.J., Zhou, C.Y., Liu, X.L., 1997. Traced study on stress field region in Yishu fault zone and Cangshan earthquake of Ms 5.2. Journal of Seismology, (3): 18-23(in Chinese with English abstract). [10] Huber, K., Fuchs, K., Palmer, J., et al., 1997. Analysis of bore hole televiewer measurements in the Vorotilov drillhole, Russia— First results. Tectonophysics, 275(1-3): 261-272. [11] Jaeger, J.C., 1961. Elasticity, fracture and flow. Methuen, London, 212. [12] Li, P.W., Cui, J.W., Wang, L.J., et al., 2005. The determination of in-situ stress from wellbore breakouts in the main borehole of the Chinese Continental Scientific Drilling. Acta Petrologica Sinica, 21(2): 421-426(in Chinese with English abstract). [13] Liu, G.X., Ma, T.Z., Huang, P.Y., et al., 1989. Major active faults and their recent movement. In: Lithospheric dynamics atlas of China. China Cartographic Publishing House, Beijing(in Chinese). [14] Plumb, R.A., 1982. Breakouts in the geothermal well, Auburn, N.Y. . Eos Trans. AGU, 63: 1118. [15] Plumb, R. A., Hickman, S. H., 1985. Stress-induced hore hole elongation: A comparison between the four-arm dipmeter and the borehole televiewer in the Auburn geothermal well. Journal of Geophysical Research, 90(B7): 5513-5521. doi: 10.1029/JB090iB07p05513 [16] Vernik, L., Zoback, M.D., 1992. Estimation of the maximum horizontal principal stress magnitude from stress induced wellbore breakout in the Cajon pass scientific research borehole. Journal of Geophysical Research, 97(b4): 5109-5119. doi: 10.1029/91JB01673 [17] Wang, L.J., Li, P.W., Cui, J.W., et al., 2005. The determination of in-situ stress in the main borehole of the Chinese Continental Scientific Drilling by using acoustic emission. Geology in China, 32(2): 259-264(in Chinese with English abstract). [18] Wang, L. J., Zhang, L. R., Wang, W., et al., 1993. Deep crustal stress measurement by borehole breakouts. In: Science papers on hydrogelogy, engineering geology and environmental geology. Seismological Press, Beijing, 63-69(in Chinese with English abstract). [19] Wang, X.F., Li, Z.J., Chen, B.L., et al., 2000. Tanlu fault zone. Geological Publishing House, Beijing, China(in Chinese). [20] Wang, Y.P., Deng, Q.D., Zhu, S.L., 1989. Lithospheric dynamics of North China. In: Lithospheric dynamics atlas of China. China Cartographic Publishing House, Beijing(in Chinese). [21] Xie, Z., Liu, Y.X., Hu, W.J., et al., 2002. The focal mechanism of earthquake in Henan and its adjacent region. Northwestern Seismological Journal, 24(3): 283-286(in Chinese with English abstract). [22] Yu, Y.X., Xu, Z.H., 1994. A study on orientations of horizontal principal stress in Jizhong depression using borehole break out data. Petroleum Exploration and Development, 21(2): 48-55(in Chinese with English abstract). [23] Yu, Y.X., Xu, Z.H., 1996. The inverse of the stress state of upper crust in Jizhong depression using borehole break out data from inclined well. Acta Seismologica Sinica, 18(2): 246-253(in Chinese with English abstract). [24] Zhou, C.Y., Jiang, Z.F., Wang, H.W., et al., 1994. Dynamic variation of stress field and seismic activity in juncture area of Jiangsu, Shandong and Yellow Sea. Seismo logical Research of Northeast China, 10(1): 10-17(in Chinese with English abstract). [25] Zoback, M.L., 1992. First and second-order patterns of stress in the lithosphere: The world stress map project. Journal of Geophysical Research, 97(B8): 11703-11728. doi: 10.1029/92JB00132 [26] Zoback, M.D., Moos, D., Mastin, L., et al., 1985. Well bore breakouts and in situ stress. Journal of Geophysical Research, 90(B7): 5523-5530. doi: 10.1029/JB090iB07p05523 [27] Zoback, M.L., Zoback, M.D., 1989. Tectonic stress field of the Continental United States. In: Paiser, L.C., Mooney, W.D., eds., Geophysical framework of the Continental United States: Boulder, Colorado. Geological Society of America Memoir, 172: 523-539. [28] 崔军文, 王连捷, 李朋武, 等, 2004. 中国大陆科学钻探主孔2 000 m以上脆性变形构造应力场. 岩石学报, 20(1): 73-80. https://www.cnki.com.cn/Article/CJFDTOTAL-YSXB200401005.htm [29] 丁原辰, 1992. 矿区地应力状态声发射粗估法及其应用. 煤炭工程师, (4): 50-56. https://www.cnki.com.cn/Article/CJFDTOTAL-ENER199204010.htm [30] 丁原辰, 张大伦, 1991. 声发射抹录不净现象在地应力测量中的应用. 岩石力学与工程学报, 10(4): 313-326. doi: 10.3321/j.issn:1000-6915.1991.04.001 [31] 华爱军, 周翠英, 刘西林, 1997. 沂沭带应力场动态追踪研究与苍山5.1级地震. 地震学刊, (3): 18-23. https://www.cnki.com.cn/Article/CJFDTOTAL-DZXK703.003.htm [32] 李朋武, 崔军文, 王连捷, 等, 2005. 中国大陆科学钻探主孔钻孔崩落与现场应力状态的确定. 岩石学报, 2(2): 421-426. https://www.cnki.com.cn/Article/CJFDTOTAL-YSXB200502016.htm [33] 刘光勋, 马廷著, 黄佩玉, 等, 1989. 主要活动断裂及现今运动. 中国岩石圈动力学地图集. 北京: 中国地图出版社. [34] 王连捷, 李朋武, 崔军文, 等, 2005. 中国大陆科学钻探主孔声发射法现今地应力状态的确定. 中国地质, 32(2): 259-264. doi: 10.3969/j.issn.1000-3657.2005.02.009 [35] 王连捷, 张利容, 王薇, 等, 1993. 利用钻孔崩落法确定深部地应力状态. 工程地质, 水文地质, 环境地质论文集. 北京: 地震出版社, 63-69. [36] 王小风, 李中坚, 陈柏林, 等, 2000. 郯庐断裂带. 北京: 地质出版社, 322-331. [37] 汪一鹏, 邓起东, 朱世龙, 1989. 华北地区岩石圈动力学特征. 中国岩石圈动力学地图集. 北京: 中国地图出版社. [38] 谢智, 刘尧兴, 胡卫健, 等, 2002. 河南及邻区地震的震源机制. 西北地震学报, 24(3): 283-286. https://www.cnki.com.cn/Article/CJFDTOTAL-ZBDZ200203016.htm [39] 俞言祥, 许忠淮, 1994. 用钻孔崩落法研究冀中坳陷水平主应力方向. 石油勘探与开发, 21(2): 48-55. https://www.cnki.com.cn/Article/CJFDTOTAL-SKYK402.008.htm [40] 俞言祥, 许忠淮, 1996. 利用斜井钻孔崩落资料反演冀中坳陷上地壳应力状态. 地震学报, 118(2): 246-253. https://www.cnki.com.cn/Article/CJFDTOTAL-DZXB602.012.htm [41] 周翠英, 姜早峰, 王红卫, 等, 1994. 苏鲁交界至南黄海地区应力场的动态变化与地震活动. 东北地震研究, 10(1): 10-17. https://www.cnki.com.cn/Article/CJFDTOTAL-DDYJ401.002.htm -