三河-平谷8.0级地震区浅层三维密度结构反演研究

为研究三河-平谷地震区浅层构造背景、地震孕育机理及地震与构造的关系,利用高精度重力异常数据,采用基于块体生长模式的重力三维反演算法对地震区浅层三维密度结构进行反演,并通过模拟试算验证基于块体生长模式反演方法的有效性和稳定性。高精度布格重力异常显示,三河-平谷8.0级地震位于大兴重力局部高、三河-马坊重力局部高与大厂重力低之间的交汇过渡低值区域。研究区浅层三维密度结构反演结果表明,1679年三河-平谷M8.0地震明显受NE向夏垫断裂控制,断裂两侧密度差异明显且向下延伸约10 km,推测发震部位深约10 km。     

Multi-component joint inversion of gravity gradient based on fast forward calculation

With the development of gravity gradient full tensor measurement technique , three-dimensional (3D) inversion based on gravity gradient tensor can provide more accurate information .But the forward calcula-tion of 3D full tensor sensitivity matrix is very time-consuming, which restricts its development and application . According to the symmetry of the kernel function , the authors reconstruct the underground source of geological body to avoid repeat computation of the same value , and work out the corresponding relationship between the re-sponse of geological body to the observation point and the response of reconstructed geological body to the obser -vation point .According to the relationship , rapid calculation of full tensor gravity sensitivity matrix can be achieved .The model calculation shows that this method can increase the speed of 30-45 times compared with the traditional calculation method .The sensitivity matrix is applied to the multi-component inversion of gravity gradient .The application of this method on the measured data provides the basis for the promotion of the meth -od.     

基于重力数据反演南海及邻区莫霍面深度分布特征

采用重力异常数据,在频率域内通过Parker-Oldenburg反演算法计算南海及邻区的莫霍面深度。通过EGM2008重力场模型联合GOCE重力场模型,计算研究区的自由空气重力异常;通过重力场分离,得到莫霍面起伏引起的重力异常,用于反演莫霍面深度。计算结果与地震剖面计算的莫霍面深度结果平均偏差1.59 km,标准偏差1.9 km。结果表明,南海及邻区的莫霍面深度分布具有明显的区域化构造特征。     

基于压缩感知原理的L1范数AVO反演算法

基于压缩感知稀疏信号采样与重构理论,利用AVO反演方法将传统的L2范数改变为L1范数,反演地下地层在L1范数下的稀疏脉冲反射系数。反演得到的稀疏尖峰将局部地下结构通过有限数量的层状结构的叠加来表示,能够提高纵向精度,与传统的AVO反演算法相比提高了薄层的反演效果且具有一定的抗噪性。数值模型及实际数据的结果表明,基于压缩感知原理的L1范数AVO反演方法更加准确、分辨率更高。     

����7.0������ǰ��ؿ��ܶȶ�̬�仯��ά����

?????????鳡1994??1996??6??????????????????????????????????????仯????????????????????????????????????????????????仯??????????????????????????仯???????????????7.0??????????????????????????????????????????????????????????仯???????С-?????С-?????С-????-????-??????????????????????????仯????????????-????-??С-?????С-????-????????????????????????????     

基于重力测点和段差变化的前郭5.8级震群研究与三维反演

利用东北地区近几年流动重力观测数据,分析前郭5.8级震群周边重力点值的时序变化特征,并基于重力段差变化对东北地区整体和局部重力变化引入能够反映异常显著性程度的指标量G和C值,最后对前郭5.8级震群前的重力场动态变化作三维密度反演。结果表明：1）前郭5.8级震群发震构造两侧的局部重力变化具有较明显的差异性,其变化趋势可作为识别构造活化或解耦运动的标志;2）显著性指标量G和C值能够为地震重力前兆的定量描述提供新思路;3）震前重力变化的反演结果显示地震发生在质量运移的过渡区域,垂向反演结果对本次地震的震源深度有一定的揭示意义。     

重力-地震联合反演的改进及在遂宁-阿坝剖面的应用

针对传统顺序联合反演的缺陷,通过改变地下空间的网格剖分方式和约束条件的约束方式来进行相应的改进,并将拟BP神经网络算法和小波多尺度分析应用到反演中,以提高反演结果的分辨率和收敛速度。龙门山地区遂宁-阿坝剖面的反演结果表明,本文方法具有更好的反演效果。     

青藏东缘三维Moho界面的位场遗传算法反演

针对波数域中Parker-Oldenburg位场迭代反演方法中存在向下延拓的问题,提出了利用改进的遗传算法进行三维Moho界面反演的位场遗传反演方法。通过模型试算证明了该方法的有效性和稳定性,并利用该方法反演了青藏高原东缘三维Moho界面的分布。     

滇西试验场区三维重力反演研究

本文依最新获得的滇西试验场区布格重力异常成果和近几年该区已有的人工地震测深资料，利用Ｐａｒｋｅｒ－Ｏｌｄｅｎｂｕｒｇ快速位场反演方法和文献［１］设计的反演计算程序，实施了多个密度界面的三维重力反演，首次获得了该区地壳上、中、下三界面的整体展布图像，分析了它们的细化结构，并结合本区重力异常特征讨论了界面图像的构造物理意义，结果表明：区内杨子板块地壳界面的去耦滑动远强于三江微板块，洱源—永胜一线的地壳半细颈结构是区内现代构造运动最活跃的地段。     

3 D inversion of borehole gravity data using cokriging

Borehole gravity has been used in mineral exPloration recently with the advent of slim-hole gravime-ters. It is logical to Perform inversion to utilize the information in the newly acquired data. The i...     

�����ܶȱ仯��ά���ݵ��Ŵ��㷨�о�

?????????????????????????????????????仯?????????????????÷???????Ч???????????????????????????????????????????仯????????????????????????仯????????????????4???????????????????2003??4??7??????????????????????????????仯???????????仯????????Щ???????????     

A new method to determine the start depth in gravity anomaly inversion

The start depth in gravity anomaly inversion is discussed in this paper, and a new method called “test comparing method” is presented to determine the start depth when both boring data and seismic data are not available. Anticline and syncline density interface models are designed to test the method; the calculated result agreed well with the designed value, which proved that the new method can be used to determine the start depth in the gravity anomaly inversion when boring data and seismic data are not available. Contribution No. 3370 from the Institute of Oceanology, Chinese Academy of Sciences. Project 49625609 supported by NSFC, and Young Scientist Lab. of Submarine Hydrothermal Activity Research.     

活动断层运动的重力场反演

阐述了断层错动引起地表重力变化的基本原理 ,给出了弹性半无限空间点源位错到矩形位错引起的重力变化的基本解析式 ,并将结论推广到多个断层。介绍了利用重力变化反演断层运动的基本思路。以滇西实验场作为应用实例 ,对 1996年 2月丽江 7.1级地震发生前后场区内的活动断层运动进行了演算 ,给出了断层位移解 ,并与已有资料相比较 ,结果基本一致 ,但误差分析结果表明本文位错模型还有待于改善。     

重磁力联合反演在研究复杂盆地基底构造中的应用

本文采用重磁力联合反演方法,计算了民和盆地基底深度和视磁化强度。针对该盆地特点,设计了等效和非均质重力反演模型,在约束条件下分块反演得到准确可靠的基底深度;用基底深度和航磁异常反演基底视磁化强度,根据岩石磁化率参数,判定基底岩性横向分布。该方法不仅能提供盆地基底构造形态,而且能提供基底岩性信息。     

姚安Ms6.0地震典型重力前兆孕震模型量化研究

基于闭锁剪力模式与地震位错理论,提出构造“等效位错”来描述地震的震前“闭锁剪力”,并建立震前地表重力变化与“等效位错”之间的定量化模型。考虑到粗差干扰,采用稳健-贝叶斯最小二乘方法对姚安M_S6.0地震的“等效位错”进行反演。结果表明,其孕震源为高倾角走滑型,“孕震矩”与USGS给出的矩震级接近,证明该方法可依据震前重力四象限图像对未来震源特性和孕震强度进行预测。     

2017年四川九寨沟MS7.0地震震前地壳密度变化特征

基于紧凑约束最小二乘优化反演算法,利用2014~2017年地表重力变化观测成果,反演得到2017年九寨沟MS7.0地震前震中及周边地区的地壳三维密度变化。结果表明：1）在震中及周边地区出现大范围的重力变化,并沿张掖-西宁-九寨沟一带形成最大值达9 μGal的重力变化梯度带;2） 地壳密度变化量随深度增加而增大,其中在45 km附近达到最大,为0.2 g/cm3;3） 震中及邻区范围内形成密度变化的正负四象限分布（2014-04~2017-04）,其表现的物质运移方向呈左旋走滑式分布,与震源机制解基本一致。据此推断,在震前产生了与地震孕育相关的深部构造运动,而强震易发生在物质运移汇聚地方。     

论多种数据联合反演的模式及算法

针对日益丰富的数据资料，综述了大地测量反演的发展进程，论述了多种数据联合反演的意义和必要性，给出了兼顾权比确定的大地测量、地震和地质三类数据的联合反演模式及采用水准测量、GPS和重力三种数据联合反演断层参数的具体公式，特别介绍了求解优化问题所有全局最优解的区间算法，并给出了基本算例，结果体现了区间算法无可比拟的优越性及其在大地测量反演领域中的应用前景。     

���ڶ�߶ȱ�Ե�������쳣����

????λ??С????????????????????????????????????????????????Щ???????????????????????????????????????????????????????÷???????????????κμ??輴??????????β?????????????????н????????????????????????????????????????÷??????????????     

北京大灰厂试验场重力变化特征及其解释

介绍了北京大灰厂试验场（简称京西试验场）的概况，对１９９３年３月以来重力重复观测资料进行了处理。计算了各期重力观测精度，结果显示１９９３年３月～１９９５年５月京西试验场存在一定重力变化。采用给定的模型反演了八宝山断层介质的密度变化，认为断层破碎带介质密度的变化可能是引起地表重力变化的主要因素     

NOTEA NEW METHOD TO DETERMINE THE START DEPTH IN GRAVITY ANOMALY INVERSION

INTRODUCTIONThevalueofthestartdepthDshouldbegivenfirstinthegraVityanoInalinversion.DisthecmstthicknesswherethecorresPondinggraVityanomaliszerooragivenknownvalue.TheburieddepthofthedensityinterfacecalculatedfrOmthegraVityanomalinver-siondependsonthestaItdePthD-Whnthegivenvalue0fthestartdepthissInall(large)lthewholeburieddepthofthecalculateddensityinteifaceissInall(larg),sointhegravityan0Inalyinversion,thedeterminati0nofthestartdepthDisveryimPortant.Usually,Disgivenbasedonseismic,ear…