CN107807388B - 一种基于多普勒效应的地震断层滑动速度计算方法 - Google Patents
一种基于多普勒效应的地震断层滑动速度计算方法 Download PDFInfo
- Publication number
- CN107807388B CN107807388B CN201711062435.9A CN201711062435A CN107807388B CN 107807388 B CN107807388 B CN 107807388B CN 201711062435 A CN201711062435 A CN 201711062435A CN 107807388 B CN107807388 B CN 107807388B
- Authority
- CN
- China
- Prior art keywords
- time
- station
- fault
- earthquake
- similitude
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 238000004364 calculation method Methods 0.000 title claims abstract description 20
- 230000000694 effects Effects 0.000 title claims abstract description 19
- 230000009466 transformation Effects 0.000 claims abstract description 24
- 238000002591 computed tomography Methods 0.000 claims abstract description 8
- 238000001228 spectrum Methods 0.000 claims description 140
- 238000003325 tomography Methods 0.000 claims description 31
- 230000011218 segmentation Effects 0.000 claims description 15
- 238000000034 method Methods 0.000 claims description 13
- 230000007423 decrease Effects 0.000 abstract description 2
- 230000002708 enhancing effect Effects 0.000 description 2
- 230000035772 mutation Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/622—Velocity, density or impedance
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开一种基于多普勒效应的地震断层滑动速度计算方法。包括:获取两台站的地震断层的滑动位错和滑动时间;在断层滑动时间内选择时间窗提取地震波,并对其进行短时傅里叶变换;确定具有相似性的两台站的短时傅里叶变换的变换对;确定当前具有相似性的两台站的短时傅里叶变换的变换对的断层滑动速度;确定断层位错与速度曲线下包围面积之间的相对误差;得到断层滑动速度函数。本发明无需掌握地壳介质参数,只需掌握震源位置、地震记录、记录台站位置信息,因此物理意义清晰。本发明增强了滑动速度时间定位、计算断层滑动速度的准确性,得到的地震断层的滑动速度具有突然上升,迅速下降的特点。
Description
技术领域
本发明属于地震学、地震工程学技术领域,尤其涉及一种基于多普勒效应的地震断层滑动速度计算方法。
背景技术
目前计算理论地震图经常使用的断层滑动速度函数主要有Haskell函数、钟形函数、指数函数、三角形函数等,Hisada的理论研究认为滑动速度函数具有迅速上升,相对缓慢地下降的特点。以上这些滑动速度函数的共同特点都是理论研究结果,都没有得到实测资料证实。
没有实测资料的原因是无法直接测定。如果用地面测得的资料反推地震断层地震时的滑动速度,需要知道地壳介质参数,而地球的不可入性阻碍了人类对地壳细节信息的掌握,欲得到较准确的地震时断层滑动速度,需另辟蹊径。
2014年在《地球物理学进展》发表的“用多普勒效应计算汶川余震断层滑动速度函数”以及2016年在《地球物理学进展》发表的“断层滑动速度的研究”两篇文章中提出可以用多普勒效应可以计算地震时断层滑动速度。该方法仅仅需要地震台站的记录、位置、震源位置和地震波传播速度等信息,无需掌握地壳介质信息,是值得研究的方法。
多普勒效应是当波源和接收器之间发生相对运动时,接收器接收到波的频率与波源发出频率不同的现象。日常生活中交通工具行驶速度的测定就是利用了多普勒效应。
多普勒效应可以这样进一步说明,当波源向接收器运动时,接收器接收到的频率比实际的高;当波源远离接收器运动时,接收器接收到的频率比实际低。实际上波源发出的频率是一定的,根据接收器接到的频率和波源发出的频率就可以计算出波源运动速度。
地震发生时,断层发出的某一频率地震波,由于断层的滑动使得接收到的地震波频率不同于实际发出的频率,根据多台站接收到的频率,可以计算出断层滑动速度。
采用多普勒效应计算断层滑动速度,具体计算了某余震断层滑动速度,证实了用多普勒效应确定断层滑动速度的有效性。但对该方法的运用上存在四方面的问题,首先对震源发出的某一频段、被二个接收台站接收到的地震波,在两台站表现出的特征缺乏理论分析;其次是如何具体确定二个接收台站接收到的地震波是由同一频段地震波发出的;再次是如何减少反射波和折射波干扰的问题;最后对断层滑动速度时间定位的准确性需要提高。
发明内容
针对现有技术存在的问题,本发明提供一种基于多普勒效应的地震断层滑动速度计算方法。
本发明的具体技术方案是:
一种基于多普勒效应的地震断层滑动速度计算方法,包括以下步骤:
步骤1:获取两台站的地震断层的滑动位错和滑动时间;
步骤2:在地震断层滑动时间内提取地震波,按时间分段进行短时傅里叶变换:在两台站的地震记录中,在断层滑动时间内选择时间窗提取地震波,并对全部断层滑动时间内的地震记录按时间顺序进行短时傅里叶变换,将两台站短时傅里叶变换结果成对按时间分段次序进行排列;
步骤3:按时间分段次序成对选定两台站的地震波的短时傅里叶变换,在频域内选定窗宽,确定时间分段上的两台站地震波在选定窗宽内的相似性,确定当前具有相似性的两台站的短时傅里叶变换的变换对;
步骤4:确定当前具有相似性的两台站的短时傅里叶变换的变换对的断层滑动速度;
所述当前具有相似性的两台站的短时傅里叶变换的变换对的断层滑动速度的计算公式如下所示:
其中,V是当前具有相似性的时间分段的断层滑动速率,fA是当前具有相似性的时间分段内台站A接收到的波源发出的地震波频率,∠A是当前具有相似性的时间分段内断层上的波源与接收台站A连线与断层上的波源运动方向之间的夹角;fB是当前具有相似性的时间分段内台站B接收到的波源发出的同一地震波频率,∠B是当前具有相似性的时间分段内断层上的波源与接收台站B连线与断层上的波源运动方向之间的夹角,u为地震波传播速率;
步骤5:判断是否计算完成全部地震记录时间分段的断层滑动速度,若是,执行步骤6,否则进行下一时间分段的两台站地震波在选定窗宽内的相似性判断,返回步骤3;
步骤6:以断层滑动时间和短时傅里叶变换时间为约束,在相应的时间内随机排列计算得到的全部断层滑动速度,得到断层速度-时间曲线,计算断层滑动速度-时间曲线下包围的面积S,以及断层位错与速度曲线下包围面积之间的相对误差e;
步骤7:若断层位错与速度曲线下包围面积之间的相对误差e大于相对误差的阈值,则返回步骤6,若断层位错与速度曲线下包围面积之间的相对误差e小于等于相对误差的阈值,得到断层滑动速度-时间曲线,即断层滑动速度函数。
所述步骤1中地震断层的滑动位错D和滑动时间Tf的计算公式如下所示:
Tf=2.03×10-9×M0 1/3;
其中,D是断层的滑动位错,Tf是断层的滑动时间,M0是地震矩。
所述按时间分段次序成对选定两台站的地震波的短时傅里叶变换,在频域内选定窗宽,确定时间分段上的两台站地震波在选定窗宽内的相似性,确定当前具有相似性的两台站的短时傅里叶变换的变换对的具体过程如下:
按时间分段次序成对选定两台站的地震波的短时傅里叶变换,在频域内选定窗宽,固定在第一个台站的短时傅里叶变换频谱上;在第二个台站的短时傅里叶变换频谱上以同样的窗宽,同样的步长进行滑动,计算两台站地震在该时间分段上的记录短时傅里叶变换频谱的相似性序列sm,若当前时间分段上两台站地震波的相似性序列sm的值小于等于相似性阈值时,则当前时间分段上两台站接收到的地震波地震记录有相似性;若当前时间分段上两台站地震波的相似性序列sm的值大于相似性阈值时,以同样的窗宽,同样的步长继续在第二个台站短时傅里叶变换频谱上滑动,计算两台站地震记录短时傅里叶变换频谱的相似性序列sm,直到在当前时间分段上第二个台站短时傅里叶变换频谱的全部滑动完成,若依然不具有相似性,选定下一时间分段的两台站的短时傅里叶变换的变换对,直至两台站地震记录短时傅里叶频谱的相似性序列sm的值小于等于相似性阈值,即确定当前两台站地震波在选定窗宽内时间分段次序中是具有相似性的短时傅里叶变换的变换对;
所述当前时间分段上两台站地震波的相似性序列sm的计算公式如下所示:
其中,Ai是台站A地震记录短时傅里叶频谱的第i个幅值,Bj是台站B地震记录短时傅里叶频谱的第j个幅值,p是台站A在选定窗宽内短时傅里叶频谱的序号,q是台站B在选定窗宽内短时傅里叶频谱的序号,p=q,n是选定窗宽内地震记录短时傅里叶频谱的幅值的总数目。
所述以断层滑动时间和短时傅里叶变换时间为约束,是指断层滑动速度对应的时间应限制在整个断层滑动时间内,短时傅里叶变换得到的断层滑动速度与短时傅里叶变换时间相对应。
所述计算断层位错于速度曲线下包围面积的相对误差e的计算公式如下所示:
其中,D是断层的滑动位错;S是断层速度-时间曲线下包围的面积。
本发明的有益效果:
本发明提出一种基于多普勒效应的地震断层滑动速度计算方法,该方法仅仅需要地震台站的记录、位置、震源位置和地震波传播速度等信息,无需掌握地壳介质信息,用该方法计算得到的地震断层的滑动速度具有突然上升,迅速下降的特点。这与以往使用的理想化的“三角函数”、“钟形函数”明显不同。并且对断层滑动时间内的地震波记录做短时傅里叶变换,基本上可消除反射波和折射波的影响,同时增强了滑动速度时间定位的准确性。该方法提出的用两台站地震在时间分段上的记录短时傅里叶变换频谱的相似性序列确定两台站地震波在选定窗宽内时间分段次序中是否是具有相似性的短时傅里叶变换的变换对,即确定两台站所接收到的地震波是否由同一频段地震波发出的,提高了计算断层滑动速度的准确性。
附图说明
图1本发明具体实施例中基于多普勒效应的地震断层滑动速度计算方法的流程图;
图2为本发明具体实施例中地震台站A和地震台站B记录到的某次大地震在11HZ附近南北方向的傅里叶幅值谱;
其中,(a)为地震台站A记录到的某次大地震在11HZ附近南北方向的傅里叶幅值谱;
(b)为地震台站B记录到的某次大地震在11HZ附近南北方向的傅里叶幅值谱;
图3为本发明具体实施例中地震台站A和地震台站B记录到的某次大地震在21HZ附近南北方向的傅里叶幅值谱;
其中,(a)为地震台站A记录到的某次大地震在21HZ附近南北方向的傅里叶幅值谱;
(b)为地震台站B记录到的某次大地震在21HZ附近南北方向的傅里叶幅值谱;
图4为本发明具体实施例中地震台站A和地震台站B记录到的某次大地震在23HZ附近南北方向的傅里叶幅值谱;
其中,(a)为地震台站A记录到的某次大地震在23HZ附近南北方向的傅里叶幅值谱;
(b)为地震台站B记录到的某次大地震在23HZ附近南北方向的傅里叶幅值谱;
图5为本发明具体实施例中地震台站A和地震台站B记录到的某次大地震在26HZ附近南北方向的傅里叶幅值谱;
其中,(a)为地震台站A记录到的某次大地震在26HZ附近南北方向的傅里叶幅值谱;
(b)为地震台站B记录到的某次大地震在26HZ附近南北方向的傅里叶幅值谱;
图6为本发明具体实施例中地震台站A和地震台站B记录到的某次大地震在28HZ附近南北方向的傅里叶幅值谱;
其中,(a)为地震台站A记录到的某次大地震在28HZ附近南北方向的傅里叶幅值谱;
(b)为地震台站B记录到的某次大地震在28HZ附近南北方向的傅里叶幅值谱;
图7为本发明具体实施例中地震台站A和地震台站B记录到的某次大地震在28.8HZ附近南北方向的傅里叶幅值谱;
其中,(a)为地震台站A记录到的某次大地震在28.8HZ附近南北方向的傅里叶幅值谱;
(b)为地震台站B记录到的某次大地震在28.8HZ附近南北方向的傅里叶幅值谱;
图8为本发明具体实施例中地震台站A和地震台站B记录到的某次大地震在29HZ附近南北方向的傅里叶幅值谱;
其中,(a)为地震台站A记录到的某次大地震在29HZ附近南北方向的傅里叶幅值谱;
(b)为地震台站B记录到的某次大地震在29HZ附近南北方向的傅里叶幅值谱;
图9为本发明具体实施例中地震台站A和地震台站B记录到的某次大地震在30HZ附近南北方向的傅里叶幅值谱;
其中,(a)为地震台站A记录到的某次大地震在30HZ附近南北方向的傅里叶幅值谱;
(b)为地震台站B记录到的某次大地震在30HZ附近南北方向的傅里叶幅值谱;
图10为本发明具体实施例中地震台站A和地震台站B记录到的某次大地震在31HZ附近南北方向的傅里叶幅值谱;
其中,(a)为地震台站A记录到的某次大地震在31HZ附近南北方向的傅里叶幅值谱;
(b)为地震台站B记录到的某次大地震在31HZ附近南北方向的傅里叶幅值谱;
图11为本发明具体实施例中地震台站A和地震台站B记录到的某次大地震在32HZ附近南北方向的傅里叶幅值谱;
其中,(a)为地震台站A记录到的某次大地震在32HZ附近南北方向的傅里叶幅值谱;
(b)为地震台站B记录到的某次大地震在32HZ附近南北方向的傅里叶幅值谱;
图12为本发明具体实施例中地震台站A和地震台站B记录到的某次大地震在35HZ附近南北方向的傅里叶幅值谱;
其中,(a)为地震台站A记录到的某次大地震在35HZ附近南北方向的傅里叶幅值谱;
(b)为地震台站B记录到的某次大地震在35HZ附近南北方向的傅里叶幅值谱;
图13为本发明具体实施例中地震台站A和地震台站B记录到的某次大地震在40HZ附近南北方向的傅里叶幅值谱;
其中,(a)为地震台站A记录到的某次大地震在40HZ附近南北方向的傅里叶幅值谱;
(b)为地震台站B记录到的某次大地震在40HZ附近南北方向的傅里叶幅值谱;
图14为本发明具体实施例中地震台站A和地震台站B记录到的某次大地震在43HZ附近南北方向的傅里叶幅值谱;
其中,(a)为地震台站A记录到的某次大地震在43HZ附近南北方向的傅里叶幅值谱;
(b)为地震台站B记录到的某次大地震在43HZ附近南北方向的傅里叶幅值谱;
图15为本发明具体实施例中地震台站A和地震台站B记录到的某次大地震在46HZ附近南北方向的傅里叶幅值谱;
图16为本发明具体实施例中地震台站A和地震台站B记录到的某次大地震在48HZ附近南北方向的傅里叶幅值谱;
其中,(a)为地震台站A记录到的某次大地震在48HZ附近南北方向的傅里叶幅值谱;
(b)为地震台站B记录到的某次大地震在48HZ附近南北方向的傅里叶幅值谱;
图17为本发明具体实施例中地震台站A和地震台站B记录到的某次大地震在50HZ附近南北方向的傅里叶幅值谱;
其中,(a)为地震台站A记录到的某次大地震在50HZ附近南北方向的傅里叶幅值谱;
(b)为地震台站B记录到的某次大地震在50HZ附近南北方向的傅里叶幅值谱;
图18为本发明具体实施例中地震台站A和地震台站B记录到的某次大地震在55HZ附近南北方向的傅里叶幅值谱;
其中,(a)为地震台站A记录到的某次大地震在55HZ附近南北方向的傅里叶幅值谱;
(b)为地震台站B记录到的某次大地震在55HZ附近南北方向的傅里叶幅值谱;
图19为本发明具体实施例中地震台站A和地震台站B记录到的某次大地震在58HZ附近南北方向的傅里叶幅值谱;
其中,(a)为地震台站A记录到的某次大地震在58HZ附近南北方向的傅里叶幅值谱;
(b)为地震台站B记录到的某次大地震在58HZ附近南北方向的傅里叶幅值谱;
图20为本发明具体实施例中地震台站A和地震台站B记录到的某次大地震在59HZ附近南北方向的傅里叶幅值谱;
其中,(a)为地震台站A记录到的某次大地震在59HZ附近南北方向的傅里叶幅值谱;
(b)为地震台站B记录到的某次大地震在59HZ附近南北方向的傅里叶幅值谱;
图21为本发明具体实施例中地震台站A和地震台站B记录到的某次大地震在63HZ附近南北方向的傅里叶幅值谱;
其中,(a)为地震台站A记录到的某次大地震在63HZ附近南北方向的傅里叶幅值谱;
(b)为地震台站B记录到的某次大地震在63HZ附近南北方向的傅里叶幅值谱;
图22为本发明具体实施例中地震台站A和地震台站B记录到的某次大地震在65HZ附近南北方向的傅里叶幅值谱;
其中,(a)为地震台站A记录到的某次大地震在65HZ附近南北方向的傅里叶幅值谱;
(b)为地震台站B记录到的某次大地震在65HZ附近南北方向的傅里叶幅值谱;
图23为本发明具体实施例中地震台站A和地震台站B记录到的某次大地震在68HZ附近南北方向的傅里叶幅值谱;
其中,(a)为地震台站A记录到的某次大地震在68HZ附近南北方向的傅里叶幅值谱;
(b)为地震台站B记录到的某次大地震在68HZ附近南北方向的傅里叶幅值谱;
图24为本发明具体实施例中地震台站A和地震台站B记录到的某次大地震在73HZ附近南北方向的傅里叶幅值谱;
其中,(a)为地震台站A记录到的某次大地震在73HZ附近南北方向的傅里叶幅值谱;
(b)为地震台站B记录到的某次大地震在73HZ附近南北方向的傅里叶幅值谱;
图25为本发明具体实施例中地震台站A和地震台站B记录到的某次大地震在76HZ附近南北方向的傅里叶幅值谱;
其中,(a)为地震台站A记录到的某次大地震在76HZ附近南北方向的傅里叶幅值谱;
(b)为地震台站B记录到的某次大地震在76HZ附近南北方向的傅里叶幅值谱;
图26为本发明具体实施例中地震台站A和地震台站B记录到的某次大地震在80HZ附近南北方向的傅里叶幅值谱;
其中,(a)为地震台站A记录到的某次大地震在80HZ附近南北方向的傅里叶幅值谱;
(b)为地震台站B记录到的某次大地震在80HZ附近南北方向的傅里叶幅值谱;
图27为本发明具体实施例中地震台站A和地震台站B记录到的某次大地震在87HZ附近南北方向的傅里叶幅值谱;
其中,(a)为地震台站A记录到的某次大地震在87HZ附近南北方向的傅里叶幅值谱;
(b)为地震台站B记录到的某次大地震在87HZ附近南北方向的傅里叶幅值谱;
图28为本发明具体实施例中某次地震断层向北方向滑动速度随时间的变化;
其中,(a)为0-5秒内某次地震断层向北方向滑动速度滑动速度随时间的变化;
(b)为6-17秒内某次地震断层向北方向滑动速度滑动速度随时间的变化;
(c)为18-40秒内某次地震断层向北方向滑动速度滑动速度随时间的变化;
(d)为43-55秒内某次地震断层向北方向滑动速度滑动速度随时间的变化;
图29为本发明具体实施例中某次地震断层向南方向滑动速度随时间的变化;
其中,(a)为0-5秒内某次地震断层向南方向滑动速度滑动速度随时间的变化;
(b)为6-17秒内某次地震断层向南方向滑动速度滑动速度随时间的变化;
(c)为18-40秒内某次地震断层向南方向滑动速度滑动速度随时间的变化;
(d)为43-55秒内某次地震断层向南方向滑动速度滑动速度随时间的变化。
具体实施方式
下面结合附图和实例对本发明做进一步说明。以下实施例用于说明本发明,但不用来限制本发明的范围。
一种基于多普勒效应的地震断层滑动速度计算方法,如图1所示,具体方法如下所述。
步骤1:获取得到两台站的地震断层的滑动位错和滑动时间;
本实施方式中,由相关文献获取某境内发生Ms8.0地震,震中位于31.0°N,103.4°E,震源深度为15公里。在地震台站A(31.54°N,103.69°E)和地震台站B(31.89°N,105.26°E)均记录到了本次地震。
地震破裂滑动持续了90S,主要位错发生在近60S的四个时间段内,地震过程中释放的标量地震矩为9.4×1020Nm。
步骤2:在地震断层滑动时间内提取地震波,按时间分段进行短时傅里叶变换:在两台站的地震记录中,在断层滑动时间内选择时间窗提取地震波,并对全部断层滑动时间内的地震记录按时间顺序进行短时傅里叶变换,将两台站短时傅里叶变换结果成对按时间分段次序进行排列;
此处对断层滑动时间内取s波记录做短时傅里叶变换,选择时间窗为3S,基本上可消除反射波和折射波的影响,同时增强了滑动速度时间定位的准确性。详细信息和短时傅里叶变换分段信息见表1。
表1短时傅里叶变换分段信息
图2到图27为此次大地震s波到达至90S时间内地震台站A和地震台站B记录短时傅里叶变换谱。
步骤3:按时间分段次序成对选定两台站的地震波的短时傅里叶变换,在频域内选定窗宽,确定时间分段上的两台站地震波在选定窗宽内的相似性,确定当前具有相似性的两台站的短时傅里叶变换的变换对;
在本实施方式中,按时间分段次序成对选定两台站的地震波的短时傅里叶变换,在频域内选定窗宽,固定在第一个台站的短时傅里叶变换频谱上;在第二个台站的短时傅里叶变换频谱上以同样的窗宽,同样的步长进行滑动,计算两台站地震在该时间分段上的记录短时傅里叶变换频谱的相似性序列sm,若当前时间分段上两台站地震波的相似性序列sm的值小于等于相似性阈值时,则当前时间分段上两台站接收到的地震波地震记录有相似性;若当前时间分段上两台站地震波的相似性序列sm的值大于相似性阈值时,以同样的窗宽,同样的步长继续在第二个台站短时傅里叶变换频谱上滑动,计算两台站地震记录短时傅里叶变换频谱的相似性序列sm,直到在当前时间分段上第二个台站短时傅里叶变换频谱的全部滑动完成,若依然不具有相似性,选定下一时间分段的两台站的短时傅里叶变换的变换对,直至两台站地震记录短时傅里叶频谱的相似性序列sm的值小于等于相似性阈值,即确定当前两台站地震波在选定窗宽内时间分段次序中是具有相似性的短时傅里叶变换的变换对;
所述当前时间分段上两台站地震波的相似性序列sm的计算公式如公式(1)所示:
其中,Ai是台站A地震记录短时傅里叶频谱的第i个幅值,Bj是台站B地震记录短时傅里叶频谱的第j个幅值,p是台站A在选定窗宽内短时傅里叶频谱的序号,q是台站B在选定窗宽内短时傅里叶频谱的序号,p=q,n是选定窗宽内地震记录短时傅里叶频谱的幅值的总数目。
当序列sm取最小值,且小于等于相似性阈值时,两台站接收到的地震记录有相似性,即两台站的接收到的地震记录是由相同频段地震波发出的,其中,相似性阈值为8‰。
步骤4:确定当前具有相似性的两台站的短时傅里叶变换的变换对的断层滑动速度;
当前具有相似性的两台站的短时傅里叶变换的变换对的断层滑动速度的计算公式如式(2)所示:
其中,V是当前具有相似性的时间分段的断层滑动速率,fA是当前具有相似性的时间分段内台站A接收到的波源发出的地震波频率,∠A是当前具有相似性的时间分段内断层上的波源与接收台站A连线与断层上的波源运动方向之间的夹角;fB是当前具有相似性的时间分段内台站B接收到的波源发出的同一地震波频率,∠B是当前具有相似性的时间分段内断层上的波源与接收台站B连线与断层上的波源运动方向之间的夹角,u为地震波传播速率。
本实施例方式中,图2到图27为地震台站A和地震台站B具有相似性的傅里叶谱。
图2为地震台站A和地震台站B记录到的某次大地震在11HZ附近南北方向的傅里叶幅值谱,根据公式(2),两者是同一波段的地震波发出的。根据公式(2)计算得到对应的滑动速度为-60.1m/s。
下面的计算都得到公式(1)验证,故不赘述。
图3为地震台站A和地震台站B记录到的某次大地震在21HZ附近南北方向的傅里叶幅值谱,根据公式(2)计算得到对应的滑动速度为-143.6m/s。
图4为地震台站A和地震台站B记录到的此次大地震在23HZ附近南北方向的傅里叶幅值谱,根据公式(2)计算得到对应的滑动速度为0m/s。
图5为地震台站A和地震台站B记录到的此次大地震在26HZ附近南北方向的傅里叶幅值谱,根据公式(2)计算得到对应的滑动速度为0m/s。
图6为地震台站A和地震台站B记录到的此次大地震在28HZ附近南北方向的傅里叶幅值谱,根据公式(2)计算得到对应的滑动速度为1.3m/s。
图7为地震台站A和地震台站B记录到的此次大地震在28.8HZ附近南北方向的傅里叶幅值谱,根据公式(2)计算得到对应的滑动速度为2.7m/s。
图8为地震台站A和地震台站B记录到的此次大地震在29HZ附近南北方向的傅里叶幅值谱,根据公式(2)计算得到对应的滑动速度为0m/s。
图9为地震台站A和地震台站B记录到的此次大地震在30HZ附近南北方向的傅里叶幅值谱,根据公式(2)计算得到对应的滑动速度为0m/s。
图10为地震台站A和地震台站B记录到的此次大地震在30HZ附近南北方向的傅里叶幅值谱,根据公式(2)计算得到对应的滑动速度为-7.9m/s。
图11为地震台站A和地震台站B记录到的此次大地震在32HZ附近南北方向的傅里叶幅值谱,根据公式(2)计算得到对应的滑动速度为-0.01m/s。
图12为地震台站A和地震台站B记录到的此次大地震在35HZ附近南北方向的傅里叶幅值谱,根据公式(2)计算得到对应的滑动速度为13.8m/s。
图13为地震台站A和地震台站B记录到的此次大地震在40HZ附近南北方向的傅里叶幅值谱,根据公式(2)计算得到对应的滑动速度为-1.2m/s。
图14为地震台站A和地震台站B记录到的此次大地震在40HZ附近南北方向的傅里叶幅值谱,根据公式(2)计算得到对应的滑动速度为2.3m/s。
图15为地震台站A和地震台站B记录到的此次大地震在46HZ附近南北方向的傅里叶幅值谱,根据公式(2)计算得到对应的滑动速度为-16.7m/s。
图16为地震台站A和地震台站B记录到的此次大地震在48HZ附近南北方向的傅里叶幅值谱,根据公式(2)计算得到对应的滑动速度为-10.9m/s。
图17为地震台站A和地震台站B记录到的此次大地震在50HZ附近南北方向的傅里叶幅值谱,根据公式(2)计算得到对应的滑动速度为30.4m/s。
图18为地震台站A和地震台站B记录到的此次大地震在55HZ附近南北方向的傅里叶幅值谱,根据公式(2)计算得到对应的滑动速度为2.7m/s。
图19为地震台站A和地震台站B记录到的此次大地震在58HZ附近南北方向的傅里叶幅值谱,根据公式(2)计算得到对应的滑动速度为-148.0m/s。
图20为地震台站A和地震台站B记录到的此次大地震在59HZ附近南北方向的傅里叶幅值谱,根据公式(2)计算得到对应的滑动速度为-4.4m/s。
图21为地震台站A和地震台站B记录到的此次大地震在63HZ附近南北方向的傅里叶幅值谱,根据公式(2)计算得到对应的滑动速度为113.2m/s。
图22为地震台站A和地震台站B记录到的此次大地震在65HZ附近南北方向的傅里叶幅值谱,根据公式(2)计算得到对应的滑动速度为-15.9m/s。
图23为地震台站A和地震台站B记录到的此次大地震在68HZ附近南北方向的傅里叶幅值谱,根据公式(2)计算得到对应的滑动速度为189.4m/s。
图24为地震台站A和地震台站B记录到的此次大地震在73HZ附近南北方向的傅里叶幅值谱,根据公式(2)计算得到对应的滑动速度为16.4m/s。
图25为地震台站A和地震台站B记录到的此次大地震在76HZ附近南北方向的傅里叶幅值谱,根据公式(2)计算得到对应的滑动速度为106.1m/s。
图26为地震台站A和地震台站B记录到的此次大地震在80HZ附近南北方向的傅里叶幅值谱,根据公式(2)计算得到对应的滑动速度为87.4m/s。
图27为地震台站A和地震台站B记录到的此次大地震在87HZ附近南北方向的傅里叶幅值谱,根据公式(2)计算得到对应的滑动速度为-92.3m/s。
步骤5:判断是否计算完成全部地震记录时间分段的断层滑动速度,若是,执行步骤6,否则进行下一时间分段的两台站地震波在选定窗宽内的相似性判断,返回步骤3;
步骤6:以断层滑动时间和短时傅里叶变换时间为约束,在相应的时间内随机排列计算得到的全部断层滑动速度,得到断层速度-时间曲线,计算断层滑动速度-时间曲线下包围的面积S,以及断层位错与速度曲线下包围面积之间的相对误差e;其计算公式如式(3)所示:
其中,D是断层的滑动位错;S是断层速度-时间曲线下包围的面积。
步骤7:若断层位错与速度曲线下包围面积之间的相对误差e大于相对误差的阈值,则返回步骤6,若断层位错与速度曲线下包围面积之间的相对误差e小于等于相对误差的阈值,得到断层滑动速度-时间曲线,即断层滑动速度函数,其中相对误差的阈值为5%。
图28为此次地震断层沿向北方向滑动速度随时间的变化,其中,(a)为0-5秒内此次地震断层向北方向滑动速度滑动速度随时间的变化、(b)为6-17秒内此次地震断层向北方向滑动速度滑动速度随时间的变化、(c)为18-40秒内此次地震断层向北方向滑动速度滑动速度随时间的变化、(d)为43-55秒内此次地震断层向北方向滑动速度滑动速度随时间的变化,其中滑动速度较小的时间段没有画出。滑动速度的显著特点是速度的突变,绝大多数时间滑动速度都不大,有些时间段断层没有滑动(如图4,5,8,9)。
图29为此次地震断层向南方向滑动速度随时间的变化,其中,(a)为0-5秒内此次地震断层向南方向滑动速度滑动速度随时间的变化、(b)为6-17秒内此次地震断层向南方向滑动速度滑动速度随时间的变化、(c)为18-40秒内此次地震断层向南方向滑动速度滑动速度随时间的变化、(d)为43-55秒内此次地震断层向南方向滑动速度滑动速度随时间的变化,其中滑动速度较小的时间段没有画出。滑动速度的显著特点是速度的突变,绝大多数时间滑动速度都不大,有些时间段断层没有滑动。
根据张勇等在2008年《中国科学D辑:地球科学》中发表的文章“某次大地震的时空破裂过程”的反演结果,此次地震发生后的0-14S内大约释放了全部地震矩的9%;15-34S内释放了全部地震矩的60%;34-43S内释放了全部地震矩的8%;43-58S内释放了全部地震矩的17%。我们计算结果在对应时间内都有一次滑动速度峰值出现,说明了计算结果的可靠性。
此次地震断层沿东西方向的滑动速度与图28和图29的情况类似。限于篇幅,不再列出。
Claims (3)
1.一种基于多普勒效应的地震断层滑动速度计算方法,其特征在于,包括以下步骤:
步骤1:获取两台站的地震断层的滑动位错和滑动时间;
步骤2:在地震断层滑动时间内提取地震波,按时间分段进行短时傅里叶变换:在两台站的地震记录中,在断层滑动时间内选择时间窗提取地震波,并对全部断层滑动时间内的地震记录按时间顺序进行短时傅里叶变换,将两台站短时傅里叶变换结果成对按时间分段次序进行排列;
步骤3:按时间分段次序成对选定两台站的地震波的短时傅里叶变换,在频域内选定窗宽,确定时间分段上的两台站地震波在选定窗宽内的相似性,确定当前具有相似性的两台站的短时傅里叶变换的变换对,具体过程如下:
按时间分段次序成对选定两台站的地震波的短时傅里叶变换,在频域内选定窗宽,固定在第一个台站的短时傅里叶变换频谱上;在第二个台站的短时傅里叶变换频谱上以同样的窗宽,同样的步长进行滑动,计算两台站地震在该时间分段上的记录短时傅里叶变换频谱的相似性序列sm,若当前时间分段上两台站地震波的相似性序列sm的值小于等于相似性阈值时,则当前时间分段上两台站接收到的地震波地震记录有相似性;若当前时间分段上两台站地震波的相似性序列sm的值大于相似性阈值时,以同样的窗宽,同样的步长继续在第二个台站短时傅里叶变换频谱上滑动,计算两台站地震记录短时傅里叶变换频谱的相似性序列sm,直到在当前时间分段上第二个台站短时傅里叶变换频谱的全部滑动完成,若依然不具有相似性,选定下一时间分段的两台站的短时傅里叶变换的变换对,直至两台站地震记录短时傅里叶频谱的相似性序列sm的值小于等于相似性阈值,即确定当前两台站地震波在选定窗宽内时间分段次序中是具有相似性的短时傅里叶变换的变换对;
所述当前时间分段上两台站地震波的相似性序列sm的计算公式如下所示:
其中,Ai是台站A地震记录短时傅里叶频谱的第i个幅值,Bj是台站B地震记录短时傅里叶频谱的第j个幅值,p是台站A在选定窗宽内短时傅里叶频谱的序号,q是台站B在选定窗宽内短时傅里叶频谱的序号,p=q,n是选定窗宽内地震记录短时傅里叶频谱的幅值的总数目;
步骤4:确定当前具有相似性的两台站的短时傅里叶变换的变换对的断层滑动速度;
所述当前具有相似性的两台站的短时傅里叶变换的变换对的断层滑动速度的计算公式如下所示:
其中,V是当前具有相似性的时间分段的断层滑动速率,fA是当前具有相似性的时间分段内台站A接收到的波源发出的地震波频率,∠A是当前具有相似性的时间分段内断层上的波源与接收台站A连线与断层上的波源运动方向之间的夹角;fB是当前具有相似性的时间分段内台站B接收到的波源发出的同一地震波频率,∠B是当前具有相似性的时间分段内断层上的波源与接收台站B连线与断层上的波源运动方向之间的夹角,u为地震波传播速率;
步骤5:判断是否计算完成全部地震记录时间分段的断层滑动速度,若是,执行步骤6,否则进行下一时间分段的两台站地震波在选定窗宽内的相似性判断,返回步骤3;
步骤6:以断层滑动时间和短时傅里叶变换时间为约束,在相应的时间内随机排列计算得到的全部断层滑动速度,得到断层速度-时间曲线,计算断层滑动速度-时间曲线下包围的面积S,以及断层位错与速度曲线下包围面积之间的相对误差e;
所述以断层滑动时间和短时傅里叶变换时间为约束,是指断层滑动速度对应的时间应限制在整个断层滑动时间内,短时傅里叶变换得到的断层滑动速度与短时傅里叶变换时间相对应;
步骤7:若断层位错与速度曲线下包围面积之间的相对误差e大于相对误差的阈值,则返回步骤6,若断层位错与速度曲线下包围面积之间的相对误差e小于等于相对误差的阈值,得到断层滑动速度-时间曲线,即断层滑动速度函数。
2.根据权利要求1所述的基于多普勒效应的地震断层滑动速度计算方法,其特征在于,所述步骤1中地震断层的滑动位错D和滑动时间Tf的计算公式如下所示:
Tf=2.03×10-9×M0 1/3;
其中,D是断层的滑动位错,Tf是断层的滑动时间,M0是地震矩。
3.根据权利要求1所述的基于多普勒效应的地震断层滑动速度计算方法,其特征在于,所述计算断层位错于速度曲线下包围面积的相对误差e的计算公式如下所示:
其中,D是断层的滑动位错,S是断层速度-时间曲线下包围的面积。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711062435.9A CN107807388B (zh) | 2017-11-02 | 2017-11-02 | 一种基于多普勒效应的地震断层滑动速度计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711062435.9A CN107807388B (zh) | 2017-11-02 | 2017-11-02 | 一种基于多普勒效应的地震断层滑动速度计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107807388A CN107807388A (zh) | 2018-03-16 |
CN107807388B true CN107807388B (zh) | 2019-06-07 |
Family
ID=61591581
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711062435.9A Expired - Fee Related CN107807388B (zh) | 2017-11-02 | 2017-11-02 | 一种基于多普勒效应的地震断层滑动速度计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107807388B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109188523B (zh) * | 2018-10-26 | 2020-03-17 | 辽宁工程技术大学 | 考虑介质对地震波吸收衰减的地震中场地反应计算方法 |
CN110244355B (zh) * | 2019-07-25 | 2021-06-08 | 西南交通大学 | 一种基于震源断层模型的脉冲地震动模拟方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1139973A (zh) * | 1993-11-18 | 1997-01-08 | 鲁阿科勒股份公司 | 尤其在硬煤煤层中优化定位回采工作面的方法 |
CN102798891A (zh) * | 2012-08-22 | 2012-11-28 | 电子科技大学 | 基于短时分数阶傅里叶变换的地震信号时频分解方法 |
CN102879813A (zh) * | 2012-09-25 | 2013-01-16 | 中国科学院地质与地球物理研究所 | 一种微地震信号到时自动拾取的方法及装置 |
CN103412331A (zh) * | 2013-08-30 | 2013-11-27 | 电子科技大学 | 一种三维地震断层自动提取方法 |
CN104335072A (zh) * | 2012-02-06 | 2015-02-04 | 离子地球物理公司 | 利用多个阵列的集成的被动和主动地震勘测 |
CN105467428A (zh) * | 2015-11-17 | 2016-04-06 | 南京航空航天大学 | 一种基于短时能量检测和频谱特征分析的地震波预警方法 |
CN105866835A (zh) * | 2016-03-28 | 2016-08-17 | 中国石油大学(华东) | 一种基于地应力分布的断层三维封闭性定量评价方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9507042B2 (en) * | 2012-08-31 | 2016-11-29 | Lumina Geophysical LLC | System and method for constrained least-squares spectral processing and analysis of seismic data |
-
2017
- 2017-11-02 CN CN201711062435.9A patent/CN107807388B/zh not_active Expired - Fee Related
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1139973A (zh) * | 1993-11-18 | 1997-01-08 | 鲁阿科勒股份公司 | 尤其在硬煤煤层中优化定位回采工作面的方法 |
CN104335072A (zh) * | 2012-02-06 | 2015-02-04 | 离子地球物理公司 | 利用多个阵列的集成的被动和主动地震勘测 |
CN102798891A (zh) * | 2012-08-22 | 2012-11-28 | 电子科技大学 | 基于短时分数阶傅里叶变换的地震信号时频分解方法 |
CN102879813A (zh) * | 2012-09-25 | 2013-01-16 | 中国科学院地质与地球物理研究所 | 一种微地震信号到时自动拾取的方法及装置 |
CN103412331A (zh) * | 2013-08-30 | 2013-11-27 | 电子科技大学 | 一种三维地震断层自动提取方法 |
CN105467428A (zh) * | 2015-11-17 | 2016-04-06 | 南京航空航天大学 | 一种基于短时能量检测和频谱特征分析的地震波预警方法 |
CN105866835A (zh) * | 2016-03-28 | 2016-08-17 | 中国石油大学(华东) | 一种基于地应力分布的断层三维封闭性定量评价方法 |
Non-Patent Citations (1)
Title |
---|
断层滑动速度的研究;李启成等;《地球物理学进展》;20161231;第31卷(第1期);第442-448页 |
Also Published As
Publication number | Publication date |
---|---|
CN107807388A (zh) | 2018-03-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Ohlmann et al. | Interpretation of coastal HF radar–derived surface currents with high-resolution drifter data | |
CN107807388B (zh) | 一种基于多普勒效应的地震断层滑动速度计算方法 | |
CN106772218B (zh) | 基于移动rfid阅读器仓库货包平面位置分级定位方法 | |
CN103869371A (zh) | 人工场源频率域全梯度电磁测量方法 | |
Garcia | Constraints on upper inner-core structure from waveform inversion of core phases | |
CN105158808A (zh) | 一种浅海瞬变电磁海空探测及其解释方法 | |
Wadhams et al. | Sea ice thickness measurement using episodic infragravity waves from distant storms | |
CN109100722A (zh) | 基于雷达回波图像扇区分量分析的风暴趋势预测方法 | |
CN105676000A (zh) | 透射式ct探地雷达对土壤相对介电常数的测定方法 | |
Hartzell et al. | Seismic site characterization of an urban sedimentary basin, Livermore valley, California: Site response, basin‐edge‐induced surface waves, and 3D simulations | |
Ito et al. | Variation of velocity and volume transport of the Tsugaru Warm Current in the winter of 1999–2000 | |
Karlsson et al. | Accumulation rates during 1311–2011 CE in North-Central Greenland derived from air-borne radar data | |
CN106954190A (zh) | 一种基于指数映射域的wifi室内定位方法 | |
CN105068075B (zh) | 一种近地面大风的计算方法 | |
Wang et al. | Characterizing the capability of mesoscale eddies to carry drifters in the northwest Pacific | |
Sripathi et al. | Characteristics of the equatorial plasma drifts as obtained by using Canadian Doppler ionosonde over southern tip of India | |
Got et al. | Strain localization and fluid migration from earthquake relocation and seismicity analysis in the western Vosges (France) | |
Guo et al. | Ground-penetrating radar survey of subsurface features at the margin of ice sheet, East Antarctica | |
Hejrani et al. | Ambient noise tomography of Australia: application to AusArray deployment | |
Bór et al. | Systematic deviations in source direction estimates of Q‐bursts recorded at Nagycenk, Hungary | |
CN106646638A (zh) | 一种固源瞬变电磁三维隧道超前预报方法 | |
Gradon et al. | Characterization with dense array data of seismic sources in the shallow part of the San Jacinto fault zone | |
CN102939548A (zh) | 确定放置于海底的探测器的位置的方法 | |
CN106646377B (zh) | 基于时间序列相似搜索的震动目标定位方法 | |
CN104777514A (zh) | 一种基于均匀水平层状介质模型的几何扩散补偿方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190607 Termination date: 20201102 |
|
CF01 | Termination of patent right due to non-payment of annual fee |