用地震波波形拟合方法研究中小地震的震源机制

采用网格搜索方法,利用区域地震三分量全波形资料,对中小地震震源机制的反演方法进行了研究。反演采用初动约束、分段求相关系数,以及幅值加权的混合方法计算误差函数,通过格点搜索法确定震源深度。理论模型试验表明,反演结果具有较好的稳定性和可靠性。云南地区区域数字地震台网的观测资料反演表明,该方法可以合理地拟合观测波形资料。将2001年7月9日(24.93°N,101.47°E)发生的MS 5.3地震的震源机制解与哈佛大学CMT进行了对比,结果表明,Z分量波形拟合略比CMT拟合的好,R分量和T分量的波形拟合明显好于CMT的对应分量。 我们对2001年4月5日在(25.77°N,102.30°E)发生的mb 4.8地震进行了震源机制反演。结果表明,该地震主要表现为走滑断层,震源机制解为(180°,80°,11°)。根据这两个地震震源机制推断的主应力轴,与该地区现代构造应力场特征基本一致。     

Resolution tests of global geodynamic models by travel-time tomography

This study presents results of a 2-D tomographic inversion of synthetic data that examines the ability of seismic tomography to reveal structures created by mantle dynamic processes. Our seismic velocity anomaly model is based on the density heterogeneities obtained from models of thermal and thermo-chemical convection. Both layered and whole-mantle models are employed to produce the synthetic input anomalies. We investigate the resolving power of the inversion of P and pP arrival times, and assess the influence of parameterisation and regularisation (damping). We show that the effect of regularisation is substantial and that the optimum damping depends upon the wavelength of the input structures. The resolution of the inversion decreases considerably at depths greater than 1000 km, therefore the ability of the kinematic inversion to distinguish between whole-mantle and layered flows (coupled via thermal coupling) may be limited.     

2010年1月15日濮阳4．2级地震震源机制反演

使用ISOLA波形反演方法,对2010年1月15日濮阳4．2级地震的震源机制进行研究,结果显示,此次地震的震源参数为：节面A：走向63°,倾角89°滑动角134°;节面B：走向153°,倾角44°,滑动角1°;标量地震矩1．0×10^15N·m,震源深度5km。研究发现,由于使用相对简化的速度模型,波形拟合程度随着震中距的增加而逐渐变差。     

Inversion of source process and related studies of the Taiwan Strait earthquake using genetic algorithm

ＩｎｖｅｒｓｉｏｎｏｆｓｏｕｒｃｅｐｒｏｃｅｓａｎｄｒｅｌａｔｅｄｓｔｕｄｉｅｓｏｆｔｈｅＴａｉｗａｎＳｔｒａｉｔｅａｒｔｈｑｕａｋｅｕｓ－ｉｎｇｇｅｎｅｔｉｃａｌｇｏｒｉｔｈｍＨＡＩ－ＪＵＮＷＡＮＧ１）（王海军）ＢＡＮＧ－ＨＵＩＬＩＮ１）（林邦慧...     

岩石物理驱动的储层裂缝参数与物性参数概率地震反演方法

长波长假设条件下,各向同性背景地层中发育一组平行排列的垂直裂缝可等效为具有水平对称轴的横向各向同性（HTI）介质.基于不同观测方位的岩石地震响应特征变化,宽方位地震数据不仅可实现裂缝岩石弹性参数与各向异性参数的预测,同时也蕴含着丰富的孔隙度等储层物性参数信息.本文结合实际地震资料提出了贝叶斯框架下岩石物理驱动的储层裂缝参数与物性参数概率地震联合反演方法,首先基于AVAZ反演裂缝岩石的弹性参数与各向异性参数,并在此基础上通过统计岩石物理模型表征孔隙度、裂缝密度等各向异性介质储层参数与裂缝岩石参数的相互关联,并采用马尔科夫链蒙特卡洛（MCMC）抽样方法进行大量样本的随机模拟,使用期望最大化（EM）算法估计后验条件概率分布,最终寻找最大后验条件概率对应的孔隙度、裂缝密度等HTI裂缝介质储层参数即为反演结果.测井及实际地震数据处理表明,该方法能够稳定合理地从方位地震资料中获取裂缝岩石弹性参数与各向异性参数,并提供了一种较为可靠的孔隙度、裂缝密度等裂缝介质储层参数概率地震反演方法.

     

纵、横波联合反演方法在海拉尔盆地贝39井区的应用

海拉尔贝39井区构造复杂,砂体边界识别难,常规纵波地震解释存在多解性,引入转换波信息可增加地震解释的可靠性.由测井曲线的岩石物理分析可知,在纵波阻抗-速度比双弹性空间可识别砂泥岩.本文根据海拉尔地区测井资料的特点,利用拟合方法重构横波测井曲线;基于海拉尔贝39井区三维三分量地震数据,采用放大比例尺层拉平的精细层位解释及纵波降颊和转换波提频相结合的方法进行精细层位匹配,在此基础上,由纵、横波联合反演求出速度比属性,定性描述了储层的平面分布,与贝39井钻遇油层特征相符,解决了该地区的油水分布问题.应用属性切片技术得到油层组的砂体平面展布,砂体的展布特征与区域地质认识的物源沉积方向具有较好的吻合程度,并对贝39并区的油藏类型进行了重新认识.     

Inversion of waveform from the data of six moderate strong earthquakes in Chinese mainland to seismic moment tensors and source mechanism

Based on Generalized Seismic Ray Theory (Helmberger, 1968), a new quickly linear inversion method from the data of seismic waveform to seismic moment tensor and source mechanism for domestic earthquake is studied in this paper. Six moderately strong earthquakes which occurred in Chinese mainland in the past few years are studied. The seismic source parameters of these earthquakes, seismic moment tensors, scalar seismic moments, fault plane solutions and source time functionsetc, are obtained. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,15, 261–268, 1993.     

Two attempts of study of seismic source from teleseismic data by simulated annealing non-linear inversion

Two seismic source studies usingteleseismic data are performed by Simulated Annealing(SA), a non-linear inversion method. The Very FastSimulated Annealing (VFSA) algorithm is used and onlyteleseismic data are inverted. We have designed a fastand efficient way of multiple direct problemevaluation, which is based on pre-calculating theelementary Green's function. During the process we setthe values of the inversion control parameters(temperature, number of iterations) and modified thecooling schedule. In the current version, theinversion seeks for the point source mechanism, thedepth of the source, the scalar moment and the sourcetime function (STF). The method is applied to twoearthquakes: 18 Nov. 1992 in Greece, M 5.9 and 14 Sep.1995 in Mexico, M7.3. The calculation is performed ona simple 1D model of the structure. For the firstearthquake the inversion recovered the solution fairlywell; for the second the solution was less acceptable.However, we do not consider this fact to be a failureof the method, but a consequence of an inadequatemodel of the source and of the medium structure. Acasual attempt of reliability determination was alsoperformed; the obtained values of errors arereasonable, except for a few cases when the methodfailed.     

复波场的实部和虚部对同时源波形反演的影响

高精度的成像需要准确的速度场信息,波形反演被认为是目前具有最高分辨率的速度反演方法之一.计算效率是目前全波形反演需要考虑的一个主要问题.为了很好地解决计算效率的问题,本文引入了一种高效的无串扰同时源反演方法,并详细介绍了其原理与计算流程.同时,基于此多震源同时反演方法,本文拓展出实波场反演、虚波场反演及复波场反演的三种反演策略,进而分析复波场的实部和虚部对全波形反演的影响.实验表明,相比于复波场全波形反演,无论是实波场反演还是虚波场反演,反演分辨率有所降低;相比于实波场反演,虚波场反演对初始模型的依赖性较小,目标函数的非线性较弱;最后,通过使用组合反演策略,即初始阶段采用虚波场反演,中后期阶段采用复波场反演,不仅可以降低反演的非线性,而且能够保证高精度建模.

     

Inversion of source process and related studies of the Taiwan Strait earthquake using genetic algorithm

Applying genetic algorithm to inversion of seismic moment tensor solution and using the data of P waveform from digital network and initial motion directions of P waves of Taiwan network stations, we studied the moment tensor solutions and focal parameters of the earthquake of M=7.3 on 16 September of 1994 in Taiwan Strait and other four quakes of M _L≥5.8 in the near region (21°–26°N, 115°–120°E). Among the five earthquakes, the quake of M=7.3 on September 16, 1994 in Taiwan Strait is the strongest one in the southeastern coast area since Nan’ao earthquake of M=7.3 in 1918. The results show that moment tensor solution of M=7.3 earthquake is mainly double-couple component, and is normal fault whose fault plane is near NW. The strike of the fault plane resembles that of the distributive bands of earthquakes before the main event and fracture pattern shown by aftershocks. The tension stress axis of focal mechanism is about horizontal, near in NE strike, the compressive stress axis is approximately vertical, near in NWW strike. It seems that this quake is controlled by the force of Philippine plate’s pressing Eurasian plate in NW direction. But from the viewpoint of P axis of near vertical and T axis of near horizontal, it is a normal fault of strong tensibility. There are relatively big difference between focal mechanism solution of this quake and those of the four other strong quakes. The complexity of source mechanism solution of these quakes represents the complexity of the process of the strait earthquake sequences. Contribution No. 98A01001, Institute of Geophysics, State Seismological Bureau, China. The subject is supported and helped by Academician Yun-Tai CHEN, Profs. Qing-Yao HONG, Zhen-Xing YAO, Tian-Yu ZHENG, Yao-Lin SHI, Ji-An XU, Bo-Shou HUANG and colleague Mei-Jian AN, Xue-Reng DING, Rui-Feng LIU. De-Chong ZHANG and Ming Li provided the digital data warm-heartedly. Lin-Ying WANG offered us the catalogue of earthquakes in southeastern coastal area in China. Xi-Li WANG and Tong-Xia BAI provided us the issued annual reports data. The authors would like to express their gratitude to all of these people. This paper is sponsored by the National Natural Science Foundation of China and Scientific and Technological Commission of Shantou, Guangdong Province.     

VSP正反演综述

垂直地震剖面法是一种发展很迅速的地球物理方法,由于垂直地震剖面法是在介质内部点上直接观测,因而能够避开或减弱剖面上部低速带的干扰及外界噪声干扰,可以更直接、更有效地研究波的运动学和动力学特征,解决其地质问题,因此在地球物理勘探领域中最为活跃.本文从垂直地震剖面正反演的角度,并结合国内外实际的例子来介绍垂直地震剖面,同时讨论了垂直地震剖面未来发展的方向.     

反演磁性地质界面的意义与方法

磁性地质界面的反演是地球物理勘探的重要课题之一,对于提高磁测资料的地质解释效果有重要意义.通过磁测资料反演磁性地质界面主要分两步：第一步,从总异常中分离有磁性界面起伏引起的磁异常;第二步,由分离的磁异常反演磁性界面的深度.本文分别从这两个方面回顾了提取异常和反演界面的较有代表性的技术.     

基于可变临界孔隙度模型的致密砂岩储层参数地震反演方法

储层参数地震反演是通过地震观测数据定量评估储层物性及含油气性的重要途径.中国西部致密砂岩油气分布范围广泛, 具有巨大的勘探开发潜力, 但其孔隙结构复杂, 储层参数的定量预测难度较大.本文提出了一种基于可变临界孔隙度模型的储层参数地震反演方法.为有效描述储层的复杂孔隙结构, 本方法将临界孔隙度作为变量, 使其具有随深度变化的特征, 以表征复杂孔隙结构, 进而提高致密砂岩储层岩石物理建模的精度.此外, 为直接从观测地震数据反演储层参数, 结合岩石物理模型和地震反射系数方程建立正演算子; 基于贝叶斯框架构建反演问题, 采用混合高斯模型构建储层参数先验约束, 用以描述储层的岩性(相)差异特征; 基于临界孔隙度、泥质含量与孔隙度的统计关系, 结合快速模拟退火算法, 在迭代优化中使临界孔隙度与储层参数交替更新, 用以提高反演结果稳定性.选用中国西部四川盆地某工区实际资料, 对本方法进行应用测试与验证, 结果表明本方法对于含气储层具有良好的指示性, 相比常规固定临界孔隙度方法, 本方法可以有效提升反演结果准确性.

     

叠前地震反演技术的进展及其在岩性油气藏勘探中的应用

与叠后地震反演相比,叠前地震反演可得到更丰富的储层信息,可提高储层的描述精度.目前,叠前地震反演主要包括弹性阻抗反演、叠前P波阻抗和S波阻抗联合反演、叠前地震波形反演.文中概述了叠前地震反演各项技术的主要进展及其典型应用实例,提出了叠前地震反演技术在岩性油气藏勘探开发中的应用策略.叠前地震反演技术在未来的油气勘探开发中将有广阔的应用前景.     

Moment tensor inversion and source rupture process of the September 27, 2003 M_S=7.9 earthquake occurred in the border area of China, Russia and Mongolia

IntroductionOn September 27, 2003, an earthquake of MS=7.9 struck the border area of China, Russia and Mongolia. According to the field investigation from the Earthquake Administration of XinjiangAutonomous Region, the whole northern Tianshan region felt the hit. Buildings and structures within six counties and one city in Altay region, which is total about 0.11×106 km2 area, were damaged to different extent and caused certain economic losses. The epicenter determined by China National …     

南北地震带南段地壳厚度重震联合最优化反演

重力反演方法是研究地壳结构和物性界面起伏的有效地球物理手段之一.本文收集了南北地震带南段67个已有的固定台站接收函数反演的Moho面深度结果,并使用基于EGM2008重力异常模型计算的布格重力异常,验证了本文提出的重震联合密度界面反演方法的有效性.利用接收函数对台站下方Moho面深度估计作为先验约束,定义了一类评价函数,通过对重力反演算法中尺度因子,平移因子和稳定性因子的最优选择,最小化重力反演结果与接收函数模型之间的差异.结果表明,本文提出的方法,可以有效地同化不同地球物理方法获得的反演模型,且通过重震联合反演可以改进由于对空间分布不均匀的接收函数结果插值可能而引起的误差.本文还通过引入Crust1.0的Moho面深度为初值,同时考虑地壳密度的横向不均匀分布,通过模型之间的联合反演有效改善了地球物理反演模型间的不一致性问题.本文反演得到的最优化Moho面深度模型与已知67个台站位置接收函数模型之间的标准差约1.9 km,小于Crust1.0与接收函数结果模型之间标准差为3.73 km的统计结果.本文研究结果对于同化重震反演结果、精化地壳密度界面模型,都具有十分重要的参考意义.     

支持向量机在重震联合反演中的应用研究

本文对利用支持向量机进行重震联合反演问题做了深入研究,特别是对支持向量机用于重震联合反演时的重力资料的处理、参数选取、特征量的提取等具体实现问题进行了讨论,最后用所设计的支持向量机重震联合反演模型对东营北部区域结晶基底岩做了预测反演,取得了满意的效果.     

矩张量反演法在江苏及附近海域的应用

利用江苏区域数字地震台网宽频带地震计记录的Pnl波波形资料,采用TDMT_INV时间域地震矩张量反演方法,尝试反演2009-2011年江苏及附近海域7次中小地震(ML3.5以上,最大震级ML5.3)的地震矩张量解,并通过记录清晰的P波初动在震源球的下半球投影对反演结果进行评价.结果表明:反演结果与初动方向一致较好.由此说明在研究区域用矩张量反演法测定中小地震震源机制解具有一定的可行性.     

Moment tensor inversion and source rupture process of the September 27, 2003 M S=7.9 earthquake occurred in the border area of China, Russia and Mongolia

We conducted moment tensor inversion and studied source rupture process for M _S=7.9 earthquake occurred in the border area of China, Russia and Mongolia on September 27 2003, by using digital teleseismic P-wave seismograms recorded by long-period seismograph stations of the global seismic network. Considering the aftershock distribution and the tectonic settings around the epicentral area, we propose that the M _S=7.9 earthquake occurred on a fault plane with the strike of 127°, the dip of 79° and the rake of 171°. The rupture process inversion result of M _S=7.9 earthquake shows that the total rupture duration is about 37 s, the scalar moment tensor is M ₀=0.97×10²⁰ N·m. Rupture mainly occurred on the shallow area with 110 km long and 30 km wide, the location in which the rupture initiated is not where the main rupture took place, and the area with slip greater than 0.5 m basically lies within 35 km deep middle-crust under the earth surface. The maximum static slip is 3.6 m. There are two distinct areas with slip larger than 2.0 m. We noticed that when the rupture propagated towards northwest and closed to the area around the M _S=7.3 hypocenter, the slip decreased rapidly, which may indicate that the rupture process was stopped by barriers. The consistence of spatial distribution of slip on the fault plane with the distribution of aftershocks also supports that the rupture is a heterogeneous process owing to the presence of barriers.