当前位置:首页 -> 焦点新闻

利用地震波特征参数解释断层落差

2005/11/25 11:18:51       

摘要:断层落差是断层解释的一项重要内容。通过研究断层与地震波之间的关系,提出了一种利用地震波特征参数解释断层落差的方法。通过对实际资料的解释,证明该方法是可行的。
关键词:地震勘探 特征参数 断层 解释
中国图书资料分类法分类号:P631.44

THE INTERPRETATION OF FAULT THROW USING CHARACTERISTIC
PARAMETERS OF SEISMIC WAVE
  
Zhao Aihua (Geophysical institute of Chinese Academy of Sciences)
  Liu Tianfang Zhang Aimin (China University of Mining and Technology)

Abstract The fault throw is an important problem for fault interpretation.By studying the relationship between fault and seismic wave,a method,in which the fault throw is interpreted using the characteristic parameters of seismic wave,is presented.It is proved that this method is feasible based on the interpretation of practical data.
Keywords seismic exploration;characteristic parameters;fault;interpretation

1 引言
  断层构造是影响矿井生产效率、威胁矿井安全的主要因素之一。断层对矿井的影响程度与其大小有直接关系。一般来说,断层越大,其影响也越大。在矿井建设和生产过程中 ,根据断层的大小,会采取不同的设计方案和预处理措施。因此准确地确定断层落差具有十分重要的意义。
  目前断层落差h主要是根据断层上、下盘断点时差Δt及地震波速度v来确定,即利用公式

2001.gif (281 字节)

然而由于多种因素的影响,准确地确定Δt和v并非易事。事实上,地震波包含着丰富的断层构造信息,其中许多信息与断层落差有着密切的联系。断层落差越大,构造运动规模一般也越大,断层带内的岩石破碎程度也越高,对地震波的吸收也就越严重。同时反射波之间的相互干涉也越强。因此断层地震波从多方面反映着断层落差大小。正是基于这种认识,我们提出了一种利用地震波特征参数解释断层落差的方法。通过对实际资料的解释,证明该方法是可行的。
2 地震波特征参数
  用地震波特征参数解释断层落差,由于涉及到二者的数量关系,所以对特征参数的要求是比较高的,不仅要求参数反映断层地震波的特征,而且还要求参数与断层落差有稳定、确切的数量关系。更进一步,这种数量关系应具有一致性,也就是说参数与断层落差的关系是单调的。因此能够用于断层落差解释的地震特征参数并不多。这次利用了3个可以用于断层落差解释的特征参数:断层时差、波形和分维差。
2.1 断层时差DFT
  断层时差DFT是由主波峰相位时间tmax曲线得到的。如图1所示,将主波峰相位时间曲线用直线段拟合,波峰相位时间在断层处的阶跃值即为断层时差DFT。断层时差与断层落差的关系比较简单,基本上为线性关系:

h=K×DFT,(1)

2003.gif (1324 字节)

图1 参数断层时差DFT示意图
细实线——实际tmax曲线;虚线——tmax拟合直线段

式中 K——大于0的常数。
  图2是根据大屯龙东矿10T线资料绘制的DFT—h关系曲线。

2004.gif (1716 字节)

图2 DFT—h关系曲线
实线——实际资料曲线;虚线——理论拟合曲线

2.2 波形FORM

2002.gif (362 字节)(2)

式中 Amax——主波峰振幅;
   T*——地震波主周期。
  波形与断层落差的关系比较复杂,总的来说落差越大,FORM越小。图3是根据大屯龙东矿10T线资料绘制的FORM—h关系曲线。

2005.gif (1734 字节)

图3 FORM—h关系曲线
实线——实际资料曲线;虚线——理论拟合曲线

2.3 分维差DFRACT
  
分维差DFRACT是从短时分形维参数FRACT(FRACT的定义见参考文献[1])曲线上得到的。如图4所示,DFRACT是断层处拟合分形维与实际分形维之差。

2006.gif (1309 字节)

图4 参数分维差示意图
实线——实际FRACT曲线;虚线——FRACT拟合曲线

  DFRACT与落差h总体为单调递增关系,具体关系则比较复杂。图5是根据大屯龙东矿10T线资料绘制的DFRACT—h关系曲线。

2007.gif (2017 字节)

图5 DFRACT—h关系曲线
  实线——实际资料曲线;虚线——理论拟合曲线

3 断层落差解释方法
3.1 关系模型的建立
  由于对地震波激发、传播等认识上的局限性和众多影响因素的存在,精确地建立地震特征参数与断层落差的关系是困难的,但是可以用一个多项式近似描述二者的关系。考虑到最小二乘法的法方程组的系数矩阵往往是病态的,因此采用离散的正交函数来进行曲线拟合。
  设有限函数系{φ0(x),φ1(x),…,φn(x)}满足:

2008.gif (802 字节)(3)

则称{φ0(x),φ1(x),…,φn(x)}为点集{x1,x2,…,xm}上关于权w(x)的离散正交函数系。
  若在正交函数系生成的函数类φ=Span{φ0(x),φ1(x),…,φn(x)}中作最小二乘曲线拟合,则有方程组

A.B=C,(4)

式中

2009.gif (1745 字节)
  2010.gif (2124 字节)

其解为 2011.gif (893 字节)
  (k=0,1,2,…,n)。
  用正交多项式进行曲线拟合,首先根据已知断层(x0,y0),(x1,y1),…,(xm,ym)(x为地震特征参数,y为断层落差)及权函数w(x)>0(本文取w(x)=1),构造关于w(x)的离散正交多项式系{Pk(x)}。用递推关系表示(nzu.gif (74 字节)m):

 2012.gif (1194 字节)(5)
 (k=1,2,…,n-1),

式中

2013.gif (2493 字节)(6)

  这样的Pk(x)是首系数为1的k次多项式。利用数学归纳法可以证明{Pk(x)}在点集{x1,x1,…,xm}上关于权函数w(x)是正交的。用正交多项式{Pk(x)}的线性组合做最小二乘拟合,可根据公式(4)、(5)逐步求Pk(x),同时计算出系数

2014.gif (823 字节)

(k=0,1,…,n)。

则可得到所求的拟合曲线方程:

2015.gif (703 字节)(7)

  这里拟合多项式的次数n可事先给定,也可在计算过程中根据误差确定。
  如果参数与落差的关系比较复杂,则拟合多项式的次数n需要很大才能满足精度要求。有时已知断层较少,难以进一步提高拟合多项式的次数,这时可以采用变量代换的方法,用低次多项式来拟合比较复杂的曲线。
  令      2016.gif (387 字节)(8)

采用同样的方法,可求得x*、y*的拟合曲线方程:

2017.gif (514 字节)(9)

3.2 落差的确定
  关系模型建立以后,计算落差就容易了。首先对地震特征参数进行变化,然后利用(9)式求得y*,最后利用(8)式求得落差y。
4 实际资料解释
  我们用波形和分维差两个参数对淮南谢桥东一采区的16条测线进行了断层解释。解释结果如表1所示,断层落差是两种参数综合解释结果。

表1 淮南谢桥东一采区测线煤解释结果

序号 测线号 断层位置
/CDP
断层落差
/m
断层落差
隶属度/%
1 EX101 16 3.3 83.2
2 EX102 17 4.4 72.8
3 EX103 16 4.5 71.9
4 EX104 14 6.8 84.4
5 EX105 15 6.1 70.0
6 EX106 16 2.4 66.9
7 EX106 45 2.4 67.8
8 EX107 17 2.8 73.3
9 EX107 42 4.9 86.7
10 EX108 17 3.5 67.1
11 EX108 42 6.4 76.0
12 EX109 44 3.0 61.2
13 EX110 44 3.7 72.1
14 EX111 45 4.0 67.2
15 EX112 44 5.9 77.7
16 EX113 44 3.8 72.3
17 EX114 44 4.5 75.6
18 EX115 44 5.4 76.2
19 EX117 42 5.6 82.1
  从地质规律来看,表1的解释结果还是较为合理的。
5 结论

  地震波包含着多种反映断层落差大小的信息。因此利用地震波特征参数解释断层落差是可行的。与传统方法相比,对于中、小断层,该方法能够更准确地确定落差。利用多个参数综合解释,可进一步提高解释的精度和可靠性。
  对于地质构造复杂区域,则需建立曲面模型,利用拟合曲面进行落差解释。

作者简介:赵爱华 男 28岁 博士研究生 地球物理勘探

作者单位:赵爱华 中国科学院地球物理研究所 北京 100101
     刘天放 张爱敏 中国矿业大学 徐州 221008

参考文献

1 吴敏金.分形信息导论.上海:上海科学技术文献出版社,1994
2 朱儒楷等.高等工程数学(上册).徐州:中国矿业大学出版社,1993
3 杨锡禄,周国铨等.中国煤炭工业百科全书:地质.测量卷.北京:煤炭工业出版社,1995
4 王桂梁等.矿井构造预测.北京:煤炭工业出版社,1993

煤炭网版权与免责声明:

凡本网注明"来源:煤炭网www.coal.com.cn "的所有文字、图片和音视频稿件,版权均为"煤炭网www.coal.com.cn "独家所有,任何媒体、网站或个人在转载使用时必须注明"来源:煤炭网www.coal.com.cn ",违反者本网将依法追究责任。

本网转载并注明其他来源的稿件,是本着为读者传递更多信息的目的,并不意味着本网赞同其观点或证实其内容的真实性。其他媒体、网站或个人从本网转载使用时,必须保留本网注明的稿件来源,禁止擅自篡改稿件来源,并自负版权等法律责任。违反者本网也将依法追究责任。 如本网转载稿件涉及版权等问题,请作者在两周内尽快来电或来函联系。

  • 用手机也能做煤炭生意啦!
  • 中煤远大:煤炭贸易也有了“支付宝”
  • 中煤开启煤炭出口贸易人民币结算新时代
  • 下半年煤炭市场依然严峻
市场动态

网站技术运营:北京真石数字科技股份有限公司、喀什中煤远大供应链管理有限公司、喀什煤网数字科技有限公司

总部地址:北京市丰台区总部基地航丰路中航荣丰1层

京ICP备18023690号-1      京公网安备 11010602010109号


关注中煤远大微信
跟踪最新行业资讯