基于全波形采样的激电多信息提取方法研究与应用     

李建华  林品荣  何畏  郑采君  丁卫忠  刘昕卓 《地球物理学进展》2020,(1):132-138

激电法是地质勘查中的一种重要方法,时、频域测量各具优势,但传统的时、频域测量实现方法相对独立.根据频域和时域可相互转换的理论,本文提出一种基于全波形采样的时频激电多参数提取的地球物理数据处理方法.采取时间域激电法的观测模式,接收机高精度同步记录整周期电压-电流全部采样点的波形数据,通过该文提出的处理方法,仅一次供电和测量即可提取多种时、频域激电参数:时域激电的视电阻率、视极化率和频域激电的多频视电阻率、视相位、视频散率、去耦后的视相位等参数.在云南保山某典型矿区开展方法的试验应用,提取了地质体大量的相关激电信息,总结了目标体的激电参数响应规律.应用结果表明,该方法便捷高效,抗干扰能力强,相较于传统的时域激电测量,多种参数组合可提升对目标地质体的认识,增强激电法的应用效果.  相似文献   

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.  相似文献   

全张量磁梯度数据的斜导数特征值边界识别方法研究   总被引：1，自引：0，他引：1       下载免费PDF全文

郭华  韩松  郑强  梁韧  段然  刘浩军  李西子 《地球物理学报》2019,62(9):3580-3590

本文提出了一种全张量磁梯度数据的斜导数边界识别方法，该方法利用全张量磁梯度数据定义了北向斜导数、东向斜导数和垂向斜导数.模型试验表明，北向和东向斜导数无法有效识别磁性异常体边界，异常形态复杂虚假异常较多，而垂向斜导数在高纬度斜磁化或者垂直磁化条件下可以清晰而准确得识别地质体的边界，具有一定的深度识别能力.与现有磁张量梯度的边界识别方法相比，识别效果好，可以有效均衡不同深度地质体的响应.将上述方法应用于大兴安岭地区实测航磁三分量数据转换得到的全张量数据，获得了研究区浅部磁性体的分布，并与三分量数据进行联合解释，取得良好的地质效果.  相似文献   

CO2 saturation estimates at Sleipner (North Sea) from seismic tomography and rock physics inversion     

Hong Yan  Bastien Dupuy  Anouar Romdhane  Børge Arntsen 《Geophysical Prospecting》2019,67(4):1055-1071

CO₂ saturations are estimated at Sleipner using a two-step imaging workflow. The workflow combines seismic tomography (full-waveform inversion) and rock physics inversion and is applied to a two-dimensional seismic line located near the injection point at Sleipner. We use baseline data (1994 vintage, before CO₂ injection) and monitor data that was acquired after 12 years of CO₂ injection (2008 vintage). P-wave velocity models are generated using the Full waveform inversion technology and then, we invert selected rock physics parameters using an rock physics inversion methodology. Full waveform inversion provides high-resolution P-wave velocity models both for baseline and monitor data. The physical relations between rock physics properties and acoustic wave velocities in the Utsira unconsolidated sandstone (reservoir formation) are defined using a dynamic rock physics model based on well-known Biot–Gassmann theories. For data prior to injection, rock frame properties (porosity, bulk and shear dry moduli) are estimated using rock physics inversion that allows deriving physically consistent properties with related uncertainty. We show that the uncertainty related to limited input data (only P-wave velocity) is not an issue because the mean values of parameters are correct. These rock frame properties are then used as a priori constraint in the monitor case. For monitor data, the Full waveform inversion results show nicely resolved thin layers of CO₂–brine saturated sandstones under intra-reservoir shale layers. The CO₂ saturation estimation is carried out by plugging an effective fluid phase in the rock physics model. Calculating the effective fluid bulk modulus of the brine–CO₂ mixture (using Brie equation in our study) is shown to be the key factor to link P-wave velocity to CO₂ saturation. The inversion tests are done with several values of Brie/patchiness exponent and show that the CO₂ saturation estimates are varying between 0.30 and 0.90 depending on the rock physics model and the location in the reservoir. The uncertainty in CO₂ saturation estimation is usually lower than 0.20. When the patchiness exponent is considered as unknown, the inversion is less constrained and we end up with values of exponent varying between 5 and 20 and up to 33 in specific reservoir areas. These estimations tend to show that the CO₂–brine mixing is between uniform and patchy mixing and variable throughout the reservoir.  相似文献   

Frequency‐domain waveform modelling and inversion for coupled media using a symmetric impedance matrix          下载免费PDF全文

Jungkyun Shin  Hyunggu Jun  Dong‐Joo Min  Changsoo Shin 《Geophysical Prospecting》2017,65(1):170-183

To simulate the seismic signals that are obtained in a marine environment, a coupled system of both acoustic and elastic wave equations is solved. The acoustic wave equation for the fluid region simulates the pressure field while minimizing the number of degrees of freedom of the impedance matrix, and the elastic wave equation for the solid region simulates several elastic events, such as shear waves and surface waves. Moreover, by combining this coupled approach with the waveform inversion technique, the elastic properties of the earth can be inverted using the pressure data obtained from the acoustic region. However, in contrast to the pure acoustic and elastic cases, the complex impedance matrix for the coupled media does not have a symmetric form because of the boundary (continuity) condition at the interface between the acoustic and elastic elements. In this study, we propose a manipulation scheme that makes the complex impedance matrix for acoustic–elastic coupled media to take a symmetric form. Using the proposed symmetric matrix, forward and backward wavefields are identical to those generated by the conventional approach; thus, we do not lose any accuracy in the waveform inversion results. However, to solve the modified symmetric matrix, LDLT factorization is used instead of LU factorization for a matrix of the same size; this method can mitigate issues related to severe memory insufficiency and long computation times, particularly for large‐scale problems.  相似文献   

Full waveform inversion of SH‐ and Love‐wave data in near‐surface prospecting          下载免费PDF全文

E. Dokter  D. Köhn  D. Wilken  D. De Nil  W. Rabbel 《Geophysical Prospecting》2017,65(Z1):216-236

We develop a two‐dimensional full waveform inversion approach for the simultaneous determination of S‐wave velocity and density models from SH ‐ and Love‐wave data. We illustrate the advantages of the SH/Love full waveform inversion with a simple synthetic example and demonstrate the method's applicability to a near‐surface dataset, recorded in the village ?achtice in Northwestern Slovakia. Goal of the survey was to map remains of historical building foundations in a highly heterogeneous subsurface. The seismic survey comprises two parallel SH‐profiles with maximum offsets of 24 m and covers a frequency range from 5 Hz to 80 Hz with high signal‐to‐noise ratio well suited for full waveform inversion. Using the Wiechert–Herglotz method, we determined a one‐dimensional gradient velocity model as a starting model for full waveform inversion. The two‐dimensional waveform inversion approach uses the global correlation norm as objective function in combination with a sequential inversion of low‐pass filtered field data. This mitigates the non‐linearity of the multi‐parameter inverse problem. Test computations show that the influence of visco‐elastic effects on the waveform inversion result is rather small. Further tests using a mono‐parameter shear modulus inversion reveal that the inversion of the density model has no significant impact on the final data fit. The final full waveform inversion S‐wave velocity and density models show a prominent low‐velocity weathering layer. Below this layer, the subsurface is highly heterogeneous. Minimum anomaly sizes correspond to approximately half of the dominant Love‐wavelength. The results demonstrate the ability of two‐dimensional SH waveform inversion to image shallow small‐scale soil structure. However, they do not show any evidence of foundation walls.  相似文献   

1D elastic full‐waveform inversion and uncertainty estimation by means of a hybrid genetic algorithm–Gibbs sampler approach          下载免费PDF全文

Mattia Aleardi  Alfredo Mazzotti 《Geophysical Prospecting》2017,65(1):64-85

Stochastic optimization methods, such as genetic algorithms, search for the global minimum of the misfit function within a given parameter range and do not require any calculation of the gradients of the misfit surfaces. More importantly, these methods collect a series of models and associated likelihoods that can be used to estimate the posterior probability distribution. However, because genetic algorithms are not a Markov chain Monte Carlo method, the direct use of the genetic‐algorithm‐sampled models and their associated likelihoods produce a biased estimation of the posterior probability distribution. In contrast, Markov chain Monte Carlo methods, such as the Metropolis–Hastings and Gibbs sampler, provide accurate posterior probability distributions but at considerable computational cost. In this paper, we use a hybrid method that combines the speed of a genetic algorithm to find an optimal solution and the accuracy of a Gibbs sampler to obtain a reliable estimation of the posterior probability distributions. First, we test this method on an analytical function and show that the genetic algorithm method cannot recover the true probability distributions and that it tends to underestimate the true uncertainties. Conversely, combining the genetic algorithm optimization with a Gibbs sampler step enables us to recover the true posterior probability distributions. Then, we demonstrate the applicability of this hybrid method by performing one‐dimensional elastic full‐waveform inversions on synthetic and field data. We also discuss how an appropriate genetic algorithm implementation is essential to attenuate the “genetic drift” effect and to maximize the exploration of the model space. In fact, a wide and efficient exploration of the model space is important not only to avoid entrapment in local minima during the genetic algorithm optimization but also to ensure a reliable estimation of the posterior probability distributions in the subsequent Gibbs sampler step.  相似文献   

Acoustic full waveform inversion of synthetic land and marine data in the Laplace domain     

Wansoo Ha  Sukjoon Pyun  Jewoo Yoo  Changsoo Shin 《Geophysical Prospecting》2010,58(6):1033-1048

Elastic waves, such as Rayleigh and mode‐converted waves, together with amplitude versus offset variations, serve as noise in full waveform inversion using the acoustic approximation. Heavy preprocessing must be applied to remove elastic effects to invert land or marine data using the acoustic inversion method in the time or frequency domains. Full waveform inversion using the elastic wave equation should be one alternative; however, multi‐parameter inversion is expensive and sensitive to the starting velocity model. We implement full acoustic waveform inversion of synthetic land and marine data in the Laplace domain with minimum preprocessing (i.e., muting) to remove elastic effects. The damping in the Laplace transform can be thought of as an automatic time windowing. Numerical examples show that Laplace‐domain acoustic inversion can yield correct smooth velocity models even with the noise originating from elastic waves. This offers the opportunity to develop an accurate smooth starting model for subsequent inversion in the frequency domain.  相似文献   

Full scale decay test of a tanker: field data and theoretical analysis     

K. Nishimoto  F. Kaster  I.Q. Masetti  J. Matsuura  J.A.P. Aranha 《Ocean Engineering》1998,26(2):125-145

This paper presents results from a full scale decay test made with a tanker in a relatively protected area in the Brazilian coast. In at least two tests the environmental loads (wind, waves and current) were very small and the time history of the surge motion was well behaved, making it possible to check some proposed models for the damping in the hull and mooring lines. Field data seem to confirm that the damping is indeed of the fluid viscosity type and the theoretical models are able to recover roughly 75% of the observed damping, the energy dissipation in the mooring lines being, by far, the major contribution. The remaining 25% are likely due to non modeled effects, such as the environment influence, which although small and not measured certainly exists, and to the friction between the mooring lines and the seabed.  相似文献   

Lie group variational integrators for the full body problem in orbital mechanics     

Taeyoung Lee  Melvin Leok  N. Harris McClamroch 《Celestial Mechanics and Dynamical Astronomy》2007,98(2):121-144

Equations of motion, referred to as full body models, are developed to describe the dynamics of rigid bodies acting under their mutual gravitational potential. Continuous equations of motion and discrete equations of motion are derived using Hamilton’s principle. These equations are expressed in an inertial frame and in relative coordinates. The discrete equations of motion, referred to as a Lie group variational integrator, provide a geometrically exact and numerically efficient computational method for simulating full body dynamics in orbital mechanics; they are symplectic and momentum preserving, and they exhibit good energy behavior for exponentially long time periods. They are also efficient in only requiring a single evaluation of the gravity forces and moments per time step. The Lie group variational integrator also preserves the group structure without the use of local charts, reprojection, or constraints. Computational results are given for the dynamics of two rigid dumbbell bodies acting under their mutual gravity; these computational results demonstrate the superiority of the Lie group variational integrator compared with integrators that are not symplectic or do not preserve the Lie group structure.  相似文献