频率域海洋可控源电磁垂直各向异性三维反演   总被引：1，自引：0，他引：1       下载免费PDF全文

彭荣华  胡祥云  李建慧  刘亚军 《地球物理学报》2019,62(6):2165-2175

地层宏观电性各向异性会对可控源电磁响应产生重要影响.由于海底地层电性结构常表现为电导率各向异性，若仅对海洋可控源电磁（MCSEM）数据进行常规各向同性反演，有可能无法获得准确的反演解释结果，从而削弱MCSEM技术的可靠性.本文实现了电导率垂直各向异性（VTI）条件下频率域海洋可控源电磁数据三维反演算法.其中，三维正演采用基于二次场控制方程的交错网格有限体积法，并利用直接矩阵分解技术来求解离散所得的大型线性方程组，有利于快速计算多场源的响应.反演采用具有近似二次收敛性的高斯牛顿算法对目标函数进行最优化.最后，对具有VTI电性各向异性特征的盐丘构造模型的MCSEM合成数据分别进行了电导率各向同性和垂直各向异性三维反演，结果表明：各向同性三维反演算法无法对受VTI介质影响的MCSEM数据进行正确的反演解释，而垂直各向异性三维反演能够获得更为可靠的地下电阻率结构和异常体分布，展现出对海底电性各向异性结构更为优良的反演解释能力.  相似文献   

GEODETIC AND TELESEISMIC CONSTRAINTS ON SLIP DISTRIBUTION OF 2015 M_W6.4 PISHAN EARTHQUAKE          下载免费PDF全文

CHEN Yun-guo  HE Ping  DING Kai-hua  LI Shui-ping  WANG Qi 《地震地质》2019,41(1):137-149

On July 3^rd, 2015, a M_W6.4 earthquake occurred on Pishan County, Xinjiang, located in the front of western Kunlun thrust belt, which is the largest earthquake(M_W6.0~7.0)in the past 40 years in this region. In this study, we collected both the near-filed geodetic coseismic deformation observations including 4 GPS sites and one high-resolution ALOS-2 InSAR imagery, and far-field teleseismic P waveforms from 25 stations provided by IRIS/USGS, to invert the fault parameters(strike and dip)and coseismic rupture model of 2015 M_W6.4 Pishan earthquake. Using the finite fault theory, a non-linear simulated annealing algorithm was employed to resolve our joint inversion problem. The strike (120°~130°) and dip angle(35°~40°)of optimal models are different from that of some previous studies, and the dip change is strongly constrained by combined data than that of strike. In fixing the geometric parameters of optimal fault model, we also considered data weight(5)(geodetic data/teleseismic P waveforms)and constrained weight from moment and smooth factor(2.5). Clearly, our results indicate that the slip distribution mainly concentrates in the depth range from 9 to 16km and a length range of 20km along the strike direction, which is similar to the spatial distribution of the relocated aftershocks. The maximum slip is~95cm. The seismic moment release is 5.45×10¹⁸N·m, corresponding to M_W6.42. Compared with the single data set, geodetic data or teleseismic waveform, our joint inversion model could simultaneously constrain the seismic moment and slip distribution well, thus avoiding effectively a lower-resolution rupture distribution determined by teleseismic-only inversion and a bias released moment estimated by the geodetic-only inversion. Importantly, we should consider both the near-field geodetic data and far-field teleseismic data in retrieving the rupture model for accurately describing the seismogenic structure of active fault in western Kunlun region.  相似文献   

COSEISMIC DISPLACEMENT FIELD OF THE WENCHUAN EARTHQUAKE DERIVED FROM STRONG MOTION RECORDS AND APPLICATION IN SLIP INVERSION          下载免费PDF全文

LIU Xiao-dong  SHAN Xin-jian  ZHANG Ying-feng  YIN Hao  QU Chun-yan 《地震地质》2019,41(4):1027-1041

The development of high-rate GNSS seismology and seismic observation methods has provided technical support for acquiring the near-field real-time displacement time series during earthquake. But in practice, the limited number of GNSS continuous stations hardly meets the requirement of near-field quasi-real-time coseismic displacement observation, while the macroseismographs could be an important complement. Compared with high-rate GNSS, macroseismograph has better sensitivity, higher resolution(100~200Hz)and larger dynamic range, and the most importantly, lower cost. However, baseline drift exists in strong-motion data, which limits its widespread use. This paper aims to prove the feasibility and reliability of strong motion data in acquiring seismic displacement sequences, as a supplement to high-rate GNSS. In this study, we have analyzed the strong-motion data of Wenchuan M_S8.0 earthquake in Longmenshan fault zone, based on the automatic scheme for empirical baseline correction proposed by Wang et al., which fits the uncorrected displacement by polynomial to obtain the fitting parameters, and then the baseline correction is completed in the velocity sequence. Through correction processing and quadratic integration, the static coseismic displacement field and displacement time series are obtained. Comparison of the displacement time series from the strong motions with the result of high-rate GPS shows a good coincidence. We have worked out the coseismic displacement field in the large area of Wenchuan earthquake using GPS data and strong motion data. The coseismic displacement fields calculated from GPS and strong motions are consistent with each other in terms of magnitude, direction and distribution patterns. High-precision coseismic deformation can provide better data constraint for fault slip inversion. To verify the influence of strong-motion data on slip distribution in Wenchuan earthquake, we used strong motion, GPS and InSAR data to estimate the stress drop, moment magnitude and coseismic slip model, and our results agreed with those of the previous studies. In addition, the inversion results of different data are different and complementary to some extent. The use of strong-motion data supplements the slip of the fault in the 180km segment and the 270~300km segment, thus making the inversion results of fault slip more comprehensive. From this result, we can draw the following conclusions:1)Based on the robust baseline correction method, the use of strong motion data, as an important complement to high-rate GNSS, can obtain reliable surface displacement after the earthquake. 2)The strong motion data provide an effective method to study the coseismic displacement sequence, the surface rupture process and quick seismogenic parameters acquisition. 3)The combination of multiple data can significantly improve the data coverage and give play to the advantages of different data. Therefore, it is suggested to combine multiple data(GPS, strong motion, InSAR, etc.)for joint inversion to improve the stability of fault slip model.  相似文献   

APPLYING 3D INVERSION OF SINGLE-PROFILE MAGNETOTELLURIC DATA TO IDENTIFY THE SHADE AND YUNONGXI FAULTS          下载免费PDF全文

JIANG Feng  CHEN Xiao-bin  DONG Ze-yi  CUI Teng-fa  LIU Zhong-yin  WANG Pei-jie 《地震地质》2019,41(6):1444-1463

Many synthetic model studies suggested that the best way to obtain good 3D interpretation results is to distribute the MT sites at a 2D grid array with regular site spacing over the target area. However, MT 3D inversion was very difficult about 10 years ago. A lot of MT data were collected along one profile and then interpreted with 2D inversion. How to apply the state-of-the-art 3D inversion technique to interpret the accumulated mass MT profiles data is an important topic. Some studies on 3D inversion of measured MT profile data suggested that 2D inversions usually had higher resolution for the subsurface than 3D inversions. Meanwhile, they often made their interpretation based on 2D inversion results, and 3D inversion results were only used to evaluate whether the overall resistivity structures were correct. Some researchers thought that 3D inversions could not resolute the local structure well, while 2D inversion results could agree with the surface geologic features much well and interpret the geologic structures easily. But in the present paper, we find that the result of 3D inversion is better than that of 2D inversion in identifying the location of the two local faults, the Shade Fault(SDF)and the Yunongxi Fault(YNXF), and the deep structures. In this paper, we first studied the electrical structure of SDF and YNXF based on a measured magnetotelluric(MT) profile data. Besides, from the point of identifying active faults, we compared the capacity of identifying deep existing faults between 2D inversion models and 3D models with different inversion parameters. The results show that both 2D and 3D inversion of the single-profile data could obtain reasonable and reliable electrical structures on a regional scale. Combining 2D and 3D models, and according to our present data, we find that both SDF and YNXF probably have cut completely the high resistivity layer in the upper crust and extended to the high conductivity layer in the middle crust. In terms of the deep geometry of the faults, at the profile's location, the SDF dips nearly vertically or dips southeast with high dip angle, and the YNXF dips southeast at depth. In addition, according to the results from our measured MT profile, we find that the 3D inversion of single-profile MT data has the capacity of identifying the location and deep geometry of local faults under present computing ability. Finally, this research suggests that appropriate cell size and reasonable smoothing parameters are important factors for the 3D inversion of single-profile MT data, more specifically, too coarse meshes or too large smoothing parameters on horizontal direction of 3D inversion may result in low resolution of 3D inversions that cannot identify the structure of faults. While, for vertical mesh size and data error thresholds, they have limited effect on identifying shallow tectonics as long as their changes are within a reasonable range. 3D inversion results also indicate that, to some extent, adding tippers to the 3D inversion of a MT profile can improve the model's constraint on the deep geometry of the outcropped faults.  相似文献   

基于三维应变格林函数反演中小震震源机制          下载免费PDF全文

周立  林强  范涛 《地震工程学报》2019,41(6):1553-1559

随着复杂速度结构反演的发展和高性能计算能力的提升,基于高精度3D介质模型计算格林函数反演震源机制更具可行性。中小地震因具有更好的覆盖和近似点源效应,在区域结构成像中有着广泛的应用。基于波形类的反演方法如波动方程层析成像\,全波形反演都需要震源机制解,而传统的震源机制反演方法不能很好地应用于中小地震。采用有限差分法构建应变格林张量(Strain Green Tensor,SGT)数据库,将合成波形和实际波形按震相截窗并滤波到不同的频带范围,先通过最小化互相关走时差来进行震源重定位,再通过最小化波形残差反演震源机制。通过合成数据测试验证方法的正确性,随后将该方法应用于青藏高原东部边缘龙门山地区,反演一系列M_W3.4～5.7的中小地震震源机制。由于应变格林张量数据库可预先构建,该方法可以应用于(近)实时震源机制解反演。  相似文献   

2017年6月16日湖北秭归Ms4.3地震成因初探     

李伟  储日升  王烁帆 《地震》2019,39(3):28-42

北京时间2017年6月16日19时48分, 湖北省秭归县发生M_S4.3地震(震中位置31.06°N, 110.48°E)。 此次地震导致当地160多间建筑物出现不同程度损坏, 100多人生活受到影响, 震区中危岩体出现一定程度破坏。 分析此次地震的震源参数将有助于了解其震源特征和发震机制。 本文首先获取了震中距在3°以内的40个宽频带固定台站的三分量数据, 然后利用FK(Frequency wavenumber)方法计算了在不同深度下Crust2.0模型和改进的1-D模型的格林函数, 分别使用Cut And Paste(CAP)方法反演得到本次地震的震源机制解和深度; 并采用基阶Rayleigh波振幅谱进一步约束了质心深度。 结果表明此次地震的最佳震源机制解的节面Ⅰ: 68°/59°/163°; 节面Ⅱ: 166°/75°/32°, 最佳矩震级为M_W4.3, 最佳深度约为5 km。 对比2013年巴东地震、 2014年秭归地震以及2011年湖北阳新—江西瑞昌地震的震源参数和余震信息, 此次地震有可能是水库诱发地震, 但要判断成因还需进一步研究。  相似文献   

基于MPI并行算法的电阻率法多种装置数据的三维联合反演          下载免费PDF全文

马欢  郭越  吴萍萍  谭捍东 《地球物理学报》2018,61(12):5052-5065

由于地表电阻率法受到浅地表局部异常体的干扰，反演精度受到影响，井中装置数据资料参与反演虽然可以减小浅地表局部异常体的干扰，但是由于钻井位置的局限性，数据量得不到保障，也会导致反演精度降低.为此，本文开发了一套结合地表、地-井、井-地和井-井装置数据的三维联合反演算法.首先，利用有限差分法实现正演模拟，采用非线性共轭梯度法（NLCG）恢复电阻率结构；其次，调用Message Passing Interface（MPI）函数库加速正演模拟和灵敏度矩阵运算，当开辟12个进程时，反演程序获得最大加速比4.51；最后，通过合成数据和实测数据算例证明该反演算法的有效性，也证实了多种装置组合数据体反演结果明显优于单一地表装置数据体反演结果.  相似文献   

2018年9月8日墨江地震及周边地区构造应力场特征分析   总被引：1，自引：0，他引：1       下载免费PDF全文

李泽潇  万永革  崔华伟  李振月  胡晓辉  黄骥超  靳志同 《地球物理学报》2020,63(4):1431-1443

2018年9月8日,云南省墨江县发生MS5.9地震并伴随一系列余震,探究该地震周围的应力场对于理解该地震的发生机制和后续地震的发展趋势具有着重要的参考意义.本研究收集了震源及其邻区中前人研究和Global CMT所给出的震源机制资料,对该地区进行了构造应力场反演,并同时利用反演得到的应力张量模拟墨江地区的震源机制解表现.结果表明:(1)在应力轴整体分布上,自西向东σ1轴(压轴)从NNE-SSW向逐渐转向NNW-SSE向,σ3轴(张轴)从WWN-EES向逐渐转向WWS-EEN向,张轴呈弧形分布,压轴呈放射状分布.(2)在应力轴倾伏角上,研究区域内的压应力轴和张应力轴倾伏角都比较小,即两轴均接近水平.(3)R值分布大体是在东南部相较于西北部大,结合当地地质背景分析得到,物质逃逸自西北向东南呈逐渐变缓的趋势.(4)利用反演得到的应力张量和应力状态计算墨江地震震源区的相对剪应力和相对正应力大小.由此推测,墨江地震恰好发生在相对剪切应力值和相对正应力正值最大的节面上.从而可以确定墨江地震的发震节面的基本参数:走向216.32°,倾角86.91°,滑动角0.27°,相对剪应力值0.9,相对正应力值0.3.本研究为此次墨江地震的发震背景和地震动力学研究提供了基础性资料.  相似文献   

Estimation of the 3D correlation structure of an alluvial aquifer from surface-based multi-frequency ground-penetrating radar reflection data     

Zhiwei Xu  James Irving  Kyle Lindsay  John Bradford  Peimin Zhu  Klaus Holliger 《Geophysical Prospecting》2020,68(2):678-689

Knowledge about the stochastic nature of heterogeneity in subsurface hydraulic properties is critical for aquifer characterization and the corresponding prediction of groundwater flow and contaminant transport. Whereas the vertical correlation structure of the heterogeneity is often well constrained by borehole information, the lateral correlation structure is generally unknown because the spacing between boreholes is too large to allow for its meaningful inference. There is, however, evidence to suggest that information on the lateral correlation structure may be extracted from the correlation statistics of the subsurface reflectivity structure imaged by surface-based ground-penetrating radar measurements. To date, case studies involving this approach have been limited to 2D profiles acquired at a single antenna centre frequency in areas with limited complementary information. As a result, the practical reliability of this methodology has been difficult to assess. Here, we extend previous work to 3D and consider reflection ground-penetrating radar data acquired using two antenna centre frequencies at the extensively explored and well-constrained Boise Hydrogeophysical Research Site. We find that the results obtained using the two ground-penetrating radar frequencies are consistent with each other, as well as with information from a number of other studies at the Boise Hydrogeophysical Research Site. In addition, contrary to previous 2D work, our results indicate that the surface-based reflection ground-penetrating radar data are not only sensitive to the aspect ratio of the underlying heterogeneity, but also, albeit to a lesser extent, to the so-called Hurst number, which is a key parameter characterizing the local variability of the fine-scale structure.  相似文献   

九寨沟Ms7.0地震的同震形变场误差研究及断层滑动分布反演     

庄春晓  王仕祥  王仲想 《测绘与空间地理信息》2020,(3):211-214,224

北京时间2017年8月8日四川省九寨沟发生Ms7.0级地震。为了提高断层滑动分布反演的可靠性,本文利用升降轨InSAR数据联合反演断层滑动分布,再根据模拟形变值二次反演,对比反演结果来探讨升降轨形变场误差和确定断层滑动分布的结果。结果表明,九寨沟升轨同震形变场的质量较好,降轨形变场受余震、震后粘弹性松弛效应影响明显,利用模拟形变值二次反演确定的断层滑动分布可靠性更高。  相似文献