地震地质 ›› 2022, Vol. 44 ›› Issue (5): 1333-1349.DOI: 10.3969/j.issn.0253-4967.2022.05.015

• 研究论文 • 上一篇    

水力压裂诱发地震活动中的b值时空异质性及其应用

姜丛1)(), 蒋长胜1),*(), 尹欣欣1),2), 王蕊嘉3), 翟鸿宇1), 张延保1), 来贵娟1), 尹凤玲1)   

  1. 1)中国地震局地球物理研究所, 北京 100081
    2)甘肃省地震局, 兰州 730000
    3)南方科技大学, 地球与空间科学系, 深圳 518055
  • 收稿日期:2021-07-30 修回日期:2022-02-28 出版日期:2022-10-20 发布日期:2022-11-28
  • 通讯作者: 蒋长胜
  • 作者简介:

    姜丛, 女, 1997年生, 2022年于中国地震局地球物理研究所获地球探测与信息技术硕士学位, 现为中国地震局地球物理研究所固体地球物理学在读博士研究生, 主要从事地震监测预测研究, E-mail:

  • 基金资助:
    国家科技基础资源调查专项(2018FY100504); 中国地震局地球物理研究所基本科研业务专项(DQJB20X11)

SPATIOTEMPORAL HETEROGENEITY AND APPLICATION OF b VALUES IN HYDRAULIC FRACTURING INDUCED SEISMICITY

JIANG Cong1)(), JIANG Chang-sheng1)(), YIN Xin-xin1),2), WANG Rui-jia3), ZHAI Hong-yu1), ZHANG Yan-bao1), LAI Gui-juan1), YIN Feng-ling1)   

  1. 1) Institute of Geophysics, China Earthquake Administration, Beijing 100081, China
    2) Gansu Earthquake Agency, Lanzhou 730000, China
    3) Department of Earth and Space Sciences, Southern University of Science and Technology, Shenzhen 518055, China
  • Received:2021-07-30 Revised:2022-02-28 Online:2022-10-20 Published:2022-11-28
  • Contact: JIANG Chang-sheng

摘要:

诱发地震及其灾害风险对非常规油气资源开发和废水回注等新型工业活动的顺利实施存在严重威胁, 并早已引发社会各界及相关科研人员的广泛关注, 由此也产生了许多需要进一步研究且亟待解决的科学问题。震级-频度分布(FMD)是描述地震活动的重要内容, 对诱发地震的b值开展系统研究有助于揭示区域应力累积状态、 地下构造特征以及诱发地震危险性。文中系统梳理了近十年来国内外工业开采水力压裂诱发地震b值的数值大小、 时空异质性、 物理机制解释、 受影响因素、 应用现状等相关的研究进展; 归纳了前人发现的可能影响b值大小的因素, 包括原位应力场、 断层几何参数、 断层成熟程度、 震源深度等客观因素, 以及注液体积、 注液速率等由具体施工情况引起的孔隙压力大小变化的主观因素; 总结了前人提出的包括孔隙压力控制论、 原位差应力大小控制论、 最大剪应力方向控制论、 地质条件非均匀性决定论等多种可能的物理机制。厘清其中已取得的共识性认识、 争议科学问题后, 我们认为诱发地震b值或许可以成为判断储层改造效果、 储层应力状态以及诱发地震灾害风险防治的入手点之一。文中对水力压裂诱发地震b值特征的总结分析可为从事新型能源开发、 地震安全监管和科学研究领域的企业、 管理人员及科学研究人员提供科学参考。

关键词: 水力压裂, 诱发地震, 时空异质性, b值, 地震危险性分析

Abstract:

Induced earthquakes and the corresponding seismic risk are rising concerns for the smooth implementation of new industrial activities such as the exploitation of unconventional oil and gas resources and has attracted broad attention from both the public and academia. As a result, many associated scientific problems need to be further examined. The magnitude-frequency distribution(FMD)is fundamental to seismicity characterization, where systematic study of b values for induced earthquakes could reveal the regional accumulated stress, subsurface structural characteristics, as well as seismic risks of induced seismicity.
In this study, we systematically reviewed the values, spatial-temporal heterogeneity, physical mechanism, dominating factors, as well as the application status of b values on hydraulic fracturing induced earthquakes in the past ten years. Multiple case analysis shows that the b value varies over a wide range(0.6~2.9)and exhibits large spatiotemporal heterogeneity. Felt earthquakes are often preceded by decreased b values and often occur in the regions with relatively low b backgrounds. In addition, fault activation or felt induced earthquakes are often accompanied by b values less than 1.0, despite that b>1.0 are commonly observed in the process of fracture expansion. Thus, the b value is promising for estimation of the state of faults(i.e., maturity and criticality).
This study further assessed the factors that may affect b values, including objective factors such as in-situ stress field, fault geometry, fault maturity and focal depth, as well as subjective factors associated specific construction conditions, such as injection volume and injection rate. We then summarized multiple possible physical mechanisms, including the pore pressure, in-situ differential stress, maximum shear stress, and the non-uniformity of geological conditions. Although the discussed factors and physical mechanisms imply the multiple complexities associated with the b value that may challenge the effectiveness of its utilization, the b value remains a preliminary but effective evaluation for first-order estimation induced-earthquake hazard. For example, in the cases of deep hypocenter, high differential stress, high fault maturity, developed initial fracture network or bedding, or when the pore pressure, fault geometry and in-situ stress field meet the conditions of high probability fault slip tendency, the b value is often less than 1.0. In fact, due to the limited understanding of the seismic risk induced by hydraulic fracturing, b value still serves as a key parameter for the seismic risk analysis such as earthquake rate prediction and maximum magnitude prediction, as well as risk control technologies such as “Traffic Light System”(TLS), the b value has been widely used in hydraulic fracturing.
Finally, we discussed the misunderstandings and challenges of b-value estimations for hydraulic fracturing induced earthquakes. For instance, b values calculated from different methods are less comparable and the quality of seismic catalogues, especially the reliability of magnitude measurement, also impact the accurate estimation and physical interpretation of b value. In addition, the mutation point of when the fault is about to reach its critical stage cannot be accurately identified through the temporal evolution of the b value alone. Even in cases where the mutation point is identified, the shut-down of current industrial operation does not guarantee the prevention of a subsequent felt event. Such challenges limit the effective utilizations of b values toward mitigating the seismic hazard associated with hydraulic fracturing induced earthquakes.
After clarifying the consensus and controversial scientific issues, we speculate that the b value for induced earthquakes may serve as one preliminary criterion for the evaluation of reservoir reconstruction, the estimation of reservoir stress state and the mitigation of induced seismicity hazard. Our study summarized and evaluated the b-value characteristics for hydraulic fracturing induced earthquakes. The paper could serve as a scientific reference for the industrial, regulating and research communities that are interested in non-conventional energy exploration and/or seismic safety supervision.

Key words: hydraulic fracturing, induced earthquake, spatiotemporal heterogeneity, b value, seismic risk analysis