天然气水合物似海底反射层的全波形反演   总被引：14，自引：4，他引：14       下载免费PDF全文

宋海斌 《地球物理学报》2003,46(1):42-46

建立了天然气水合物似海底反射层(BSR)研究的全波形反演方法. 这是一种将 水平层状弹性介质的反射共中心点道集转换为截距时间-水平慢度域的反演方法. 反演过程 中采用了全局搜索方法与非线性局部搜索方法. 分两步进行. 第一步是根据走时数据应用非 常快速模拟算法求得速度结构的长波长分量. 第二步，利用波形资料用共轭梯度法求得速度 的短波长扰动分量. 这样，最后反演得到的速度结构模型包含了长波长与短波长分量. 反演 中利用了多网格参数化技术. 日本东南海海槽双BSR的速度结构的反演表明，全波形反演是 天然气水合物BSR研究的重要手段之一.  相似文献   

A finite element approach for predicting the full resistance profile of a spudcan deeply penetrating in dense sand overlying clay     

《Applied Ocean Research》2019

This paper presents a finite element approach to calculate the full resistance profile of a spudcan deeply penetrating in dense sand overlying clay, in which a potential for an installing spudcan to experience a sudden uncontrolled punch-through failure exists. A modified Mohr-Coulomb model characterized by incorporating a four-phase variation of the mobilized strength and dilation parameters with an equivalent accumulated plastic strain is developed and tested for the overlying dense silica sand. An extended Tresca model is used for the strain softening of the underlying clay. A series of large deformation finite-element (LDFE) analyses are carried out, varying the strength and dilation parameters as well as the spudcan geometries. A fairly good performance of the present approach is verified by validating against groups of centrifuge tests data, allowing the numerical study to be extended parametrically. The four-phase variation of the mobilized strength and dilation parameters involved in the progressive failure of the upper dense sand is parametrically studied and extended to cover the range of sand relative densities that are of practical interest. Additionally, comparisons with the typical existing LDFE analyses using both simple and sophisticated constitutive models are carried out. It shows that the present approach performs fairly well to calculate the full resistance profile of a spudcan deeply penetration in both thin and thick dense sand overlying clay, especially the peak and post-peak resistance, within around 5% of the corresponding centrifuge tests results.  相似文献   

Theoretical background for the inversion of seismic waveforms including elasticity and attenuation   总被引：3，自引：0，他引：3  

Albert Tarantola 《Pure and Applied Geophysics》1988,128(1-2):365-399

To account for elastic and attenuating effects in the elastic wave equation, the stress-strain relationship can be defined through a general, anisotropic, causal relaxation function ^ijkl (x, ). Then, the wave equation operator is not necessarily symmetric (self-adjoint), but the reciprocity property is still satisfied. The representation theorem contains a term proportional to the history of strain. The dual problem consists of solving the wave equation withfinal time conditions and an anti-causal relaxation function. The problem of interpretation of seismic waveforms can be set as the nonlinear inverse problem of estimating the matter density (x) and all the functions ^ijkl (x, ). This inverse problem can be solved using iterative gradient methods, each iteration consisting of the propagation of the actual source in the current medium, with causal attenuation, the propagation of the residuals—acting as if they were sources—backwards in time, with anti-causal attenuation, and the correlation of the two wavefields thus obtained.  相似文献   

Source mechanism of the 1992 Roermond earthquake from surface-wave inversion of regional data     

Jochen Braunmiller  Torsten Dahm  Klaus-Peter Bonjer 《Geophysical Journal International》1994,116(3):663-672

  相似文献   

The Wavelength-Smoothing Method For Approximating Broad-Band Wave Propagation Through Complicated Velocity Structures     

Anthony Lomax 《Geophysical Journal International》1994,117(2):313-334

  相似文献   

Effects of Sedimentary Structure and Bathymetry Near the Source On Teleseismic P Waveforms From Shallow Subduction Zone Earthquakes     

Taro Okamoto 《Geophysical Journal International》1993,112(3):471-480

  相似文献   

全波震相分析的应用   总被引：10，自引：0，他引：10  

赵鸿儒  孙进忠等 《地球物理学报》1990,33(1):54-63

全波震相分析法是以弹性波传播理论为基础,以现代数字技术为手段,对地震波进行全面综合分析的理论和方法.它既包括对单一分量震波观测记录中的各种震相的识别、确认和分析,又包括对多分量震波记录的合成、分析和图示.全波震相分析法涉及震源、射线路径、初动、走时、振幅、波形、时域和空域的瞬态谱、质点振动矢量等诸多方面,因而能更全面地揭示地震波场与地下介质之间的关系.本文以实例展示在不同领域的研究中全波震相分析的应用情况.全波震相分析法为实现多种波型的联合应用奠定了理论基础,并提供了新的方法.  相似文献   

Interpretation of Gravity Anomalies with the Normalized Full Gradient (NFG) Method and an Example     

Ali Aydin 《Pure and Applied Geophysics》2007,164(11):2329-2344

The Normalized Full Gradient (NFG) method which was put forward about 50 years ago has been used for downward continuation of gravity potential data, especially in the former Union of Soviet Socialist Republics. This method nullifies perturbations due to the passage of mass depth during downward continuation. The method depends on the downwards analytical continuation of normalized full gradient values of gravity data. Analytical continuation discriminates certain structural anomalies which cannot be distinguished in the observed gravity field. This method has been used in various petroleum and tectonic studies. The Trapeze method was used for the determination of Fourier coefficients during the application of this method. No other techniques for calculating these coefficients have been used. However, the Filon method was used for the determination of Fourier coefficients during the application of the NFG method in this work. This method, rather than the Trapeze method, should be preferred for indicating abnormal mass resources at the lower harmonics. In this study, the NFG method using the Filon method has been applied the first time to theoretical models of gravity profiles as example field at the Hasankale-Horasan petroleum exploration province where successful results were achieved. Hydrocarbon presence was shown on the NFG sections by the application of NFG downward continuation operations on theoretical models. Important signs of hydrocarbon structure on the NFG section for field and model data at low harmonics are obtained more effectively using this method.  相似文献   

Preliminary study of crust-upper mantle structure of the Tibetan Plateau by using broadband teleseismic body waveforms   总被引：2，自引：0，他引：2  

Lu-Pei Zhu  Rong-Sheng Zeng  Francis T. Wu  Thomas J. Owens  George E. Randall 《地震学报(英文版)》1993,6(2):305-316

As part of a joint Sino-U.S. research project to study the deep structure of the Tibetan Plateau, 11 broadband digital seismic recorders were deployed on the Plateau for one year of passive seismic recording. In this report we use teleseimic P waveforms to study the seismic velocity structure of crust and upper mantle under three stations by receiver function inversion. The receiver function is obtained by first rotating two horizontal components of seismic records into radial and tangential components and then deconvolving the vertical component from them. The receiver function depends only on the structure near the station because the source and path effects have been removed by the deconvolution. To suppress noise, receiver functions calculated from events clustered in a small range of back-azimuths and epicentral distances are stacked. Using a matrix formalism describing the propagation of elastic waves in laterally homogeneous stratified medium, a synthetic receiver function and differential receiver functions for the parameters in each layer can be calculated to establish a linearized inversion for one-dimensional velocity structure. Preliminary results of three stations, Wen-quan, Golmud and Xigatze (Coded as WNDO, TUNL and XIGA), located in central, northern and southern Plateau are given in this paper. The receiver functions of all three stations show clear P-S converted phases. The time delays of these converted phases relative to direct P arrivals are: WNDO 7.9s (for NE direction) and 8.3s (for SE direction), TUNL 8.2s, XIGA 9.0s. Such long time delays indicate the great thickness of crust under the Plateau. The differences between receiver function of these three station shows the tectonic difference between southern and north-central Plateau. The waveforms of the receiver functions for WNDO and TUNL are very simple, while the receiver function of XIGA has an additional midcrustal converted phase. The S wave velocity structures at these three stations are estimated from inversions of the receiver function. The crustal shear wave velocities at WNDO and TUNL are vertically homogeneous, with value between 3.5–3.6 km/s down to Moho. This value in the lower crust is lower than the normal value for the lower crust of continents, which is consistent with the observed strong Sn attenuation in this region. The velocity structure at XIGA shows a velocity discontinuity at depth of 20 km and high velocity value of 4.0 km/s in the midcrust between 20–30 km depth. Similar results are obtained from a DSS profile in southern Tibet. The velocity under XIGA decreases below a depth of 30 km, reaching the lowest value of 3.2 km/s between 50–55 km. depth. This may imply that the Indian crust underthrusts the low part of Tibetan crust in the southern Plateau, forming a “double crust”. The crustal thickness at each of these sites is: WNDO, 68 km; TUNL, 70 km; XI-GA, 80 km. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,14, Supp., 581–592, 1992.  相似文献   

基于改进全变差正则化的GPR多尺度全波形双参数同步反演          下载免费PDF全文

王珣  冯德山  王向宇 《地球物理学报》1954,63(12):4485-4501

针对探地雷达（GPR）双参数全波形反演中电导率反演精度差、双参数存在串扰现象、反演计算量大、易陷入局部极值等问题.作者将具有多参数调节功能的L-BFGS算法引入到GPR时间域全波形反演中,它避免了对Hessian矩阵的直接存储与精确求解,减小了存储量和计算量.结合参数调节因子的选取,有效减小了同步反演时介电常数与电导率的串扰影响,在不降低介电常数反演精度的前提下,提高电导率参数的反演精度.通过在反演目标函数中加载改进全变差正则化方法,提高了反演的稳定性,使目标体边缘轮廓更加清晰.首先以简单模型为例,对比了单尺度反演与多尺度串行反演策略的优劣,说明多尺度串行反演有利于逐步搜索全局最优解;而开展参数调节因子的选取实验,说明合适的参数调节因子可以有效改善介质电导率的反演精度;测试了不同正则化的反演效果,表明改进全变差正则化能提高反演稳定性,显著降低模型重构误差.最后,分别对含噪合成数据和实测数据进行了反演测试,说明本文提出的多尺度、双参数反演具有较强的鲁棒性,能提供更丰富的信息约束,重构图像界面清晰、反演效果好.  相似文献