查看更多>>摘要:We present a goal-oriented adaptive finite element (FE) modelling algorithm for 3-D magnetotelluric fields in generally anisotropic conductivity media. The model consists of a background layered structure, containing an isotropic blocks. Each block and layer might be anisotropic by assigning to them 3 x 3 conductivity tensors. The second-order partial differential equations are solved using the adaptive finite element method (FEM). The computational domain is subdivided into unstructured tetrahedral elements, which allow for complex geometries including bathymetry and dipping interfaces. The grid refinement process is guided by a global posteriori error estimator and is performed iteratively.
查看更多>>摘要:In towed marine seismic data acquisition, a gap between the source and the nearest recording channel is typical. Therefore, extrapolation of the missing near-offset traces is often required to avoid unwanted effects in subsequent data processing steps. However, most existing interpolation methods perform poorly when extrapolating traces. Interferometric interpolation methods are one particular method that have been developed for filling in trace gaps in shot gathers. Interferometry-type interpolation methods differ from conventional interpolation methods as they utilize information from several adjacent shot records to fill in the missing traces. In this study, we aim to improve upon the results generated by conventional time-space domain interferometric interpolation by performing interferometric interpolation in the Radon domain, in order to overcome the effects of irregular data sampling and limited source-receiver aperture. We apply both time-space and Radon-domain interferometric interpolation methods to the Sigsbee2B synthetic dataset and a real towed marine dataset from the Baltic Sea with the primary aim to improve the image of the seabed through extrapolation into the near-offset gap. Radon-domain interferometric interpolation performs better at interpolating the missing near offset traces than conventional interferometric interpolation when applied to data with irregular geometry and limited source-receiver aperture. We also compare the interferometric interpolated results with those obtained using solely Radon transform (RT) based interpolation and show that interferometry-type interpolation performs better than solely RT-based interpolation when extrapolating the missing near-offset traces. After data processing, we show that the image of the seabed is improved by performing interferometry-type interpolation, especially when Radon-domain interferometric interpolation is applied. (C) 2018 Elsevier B.V. All rights reserved.
查看更多>>摘要:This paper focuses on an evaluation of rotary-percussion drilling (RPD) as a seismic source. Two field experiments were conducted to characterize seismic sources from different rocks with different strengths, i.e. weak shale and hard arkose. Characterization of RPD sources consist of spectral analysis and mean power measurements, along with field measurements of the source radiation patterns. Spectral analysis shows that increase of rock strength increases peak frequency and widens bandwidth, which makes harder rock more viable for seismic-while-drilling purposes. Mean power analysis infers higher magnitude of body waves in RPD than in conventional drillings. Within the horizontal plane, the observed P-wave energy radiation pattern partially confirms the theoretical radiation pattern under a single vertical bit vibration. However a horizontal lobe of energy is observed close to orthogonal to the axial bit vibration. From analysis, this lobe is attributed to lateral bit vibration, which is not documented elsewhere during RPD. Within the horizontal plane, the observed radiation pattern of P-waves is generally consistent with a spherically-symmetric distribution of energy. In addition, polarization analysis is conducted on P-waves recorded at surface geophones for understanding the particle motions. P-wave particle motions are predominantly in the vertical direction showing the interference of the free-surface. (C) 2018 Elsevier B.V. All rights reserved.
查看更多>>摘要:Explicit depth frequency-space (f - x) prestack imaging is an attractive mechanism for seismic imaging. To date, the main focus of this method was data migration assuming an acoustic medium, but until now very little work assumed visco-acoustic media. Real seismic data usually suffer from attenuation and dispersion effects. To compensate for attenuation in a visco-acoustic medium, new operators are required. We propose using the L-1-norm minimization technique to design visco-acousticf-x extrapolators. To show the accuracy and compensation of the operators, prestack depth migration is performed on the challenging Marmousi model for both acoustic and visco-acoustic datasets. The final migrated images show that the proposed L-1-norm extrapolation results in practically stable and improved resolution of the images. (C) 2018 Published by Elsevier B.V.
查看更多>>摘要:Petrophysical description of reservoirs requires proper knowledge of elastic parameters like P-and S wave velocities (V-p and V-s) and density (rho), which can be retrieved from pre-stack seismic data using the concept of elastic impedance (El). We propose an inversion algorithm which recovers elastic parameters from pre-stack seismic data in two sequential steps. In the first step, using the multichannel blind seismic inversion method (exploited recently for recovering acoustic impedance from post-stack seismic data), high-resolution blocky El models are obtained directly from partial angle-stacks. Using an efficient total variation (TV) regularization, each angle-stack is inverted independently in a multichannel form without prior knowledge of the corresponding wavelet. The second step involves inversion of the resulting El models for elastic parameters. Mathematically, under some assumptions, the El's are linearly described by the elastic parameters in the logarithm domain. Thus a linear weighted least squares inversion is employed to perform this step. Accuracy of the concept of elastic impedance in predicting reflection coefficients at low and high angles of incidence is compared with that of exact Zoeppritz elastic impedance and the role of low frequency content in the problem is discussed. The performance of the proposed inversion method is tested using synthetic 2D data sets obtained from the Marmousi model and also 2D field data sets. The results confirm the efficiency and accuracy of the proposed method for inversion of pre-stack seismic data. (C) 2018 Published by Elsevier B.V.
查看更多>>摘要:The pore structure in heterogeneous carbonate rock is usually very complex. This complex pore system makes the relationship between the velocity and porosity of the rock highly scattered, so that for the classical two-dimensional rock physics template (2D RPT) it is not enough to accurately describe the quantitative relationship between the rock elastic parameters of this kind of reservoir and its porosity and water saturation. Therefore it is possible to attribute the effect of pore type to that of the porosity or water saturation, and leads to great deviations when applying such a 2D RPT to predict the porosity and water saturation in seismic reservoir prediction and hydrocarbon detection. This paper first presents a method to establish a new three-dimensional rock physics template (3D RPT) by integrating the Gassmann equations and the porous rock physics model, and use it to characterize the quantitative relation between rock elastic properties and the reservoir parameters including the pore aspect ratio, porosity and water saturation, and to predict these parameters from the known elastic properties. The test results on the real logging and seismic inversion data show that the 3D RPT can accurately describe the variations of elastic properties with the porosity, water saturation and pore -structure parameters, and effectively improve the accuracy of reservoir parameters prediction. (C) 2018 Elsevier B.V. All rights reserved.
查看更多>>摘要:Full waveform inversion (FWI) of ground penetrating radar (GPR) is a promising technique to quantitatively evaluate the permittivity and conductivity of near subsurface. However, these two parameters are simultaneously inverted in the GPR FWI, increasing the difficulty to obtain accurate inversion results for both parameters. In this study, I present a structural constrained GPR FWI procedure to jointly invert the two parameters, aiming to force a structural relationship between permittivity and conductivity in the process of model reconstruction. The structural constraint is enforced by a cross-gradient function. In this procedure, the permittivity and conductivity models are inverted alternately at each iteration and updated with hierarchical frequency components in the frequency domain. The joint inverse problem is solved by the truncated Newton method which considering the effect of Hessian operator and using the approximated solution of Newton equation to be the perturbation model in the updating process. The joint inversion procedure is tested by three synthetic examples. The results show that jointly inverting permittivity and conductivity in GPR FWI effectively increases the structural similarities between the two parameters, corrects the structures of parameter models, and significantly improves the accuracy of conductivity model, resulting in a better inversion result than the individual inversion. (C) 2018 Elsevier B.V. All rights reserved.
查看更多>>摘要:This paper reports on the fundamental role played by Ground Penetrating Radar (GPR), alongside advanced processing and presentation methods, during the tunnel boring project at a Dam and Hydro -Electric Power Station. It identifies from collected GPR data such issues as incomplete grouting and the presence of karst conduits and voids and provides full details of the procedures adopted. In particular, the application of collected GPR data to the Neural Network (NN) method is discussed. (C) 2018 Elsevier B.V. All rights reserved.
Leite, David NakamuraBortolozo, Cassiano AntonioPorsani, Jorge LuisCouto, Marco Antonio, Jr....
16页
查看更多>>摘要:Although Brazil is well known by the large rivers and the Amazon Rain Forest most cities do not have access to sufficient quantities of surface water to supply the population. Because of this 61% of Brazilian population (IBGE, 2003) depends on groundwater resources. In order to help the conscious exploration of this resource in Urupes city (Sao Paulo State) which is characterized by problems of lack of water, this research applied the transient electromagnetic method (TDEM) and Vertical Electrical Sounding (VES) for the geoelectrical characterization of the interest region. So, the objective of this work was increase the hydrogeological basis for groundwater exploitation of Bauru sedimentary aquifer and Serra Geral fractured aquifer (Parana Basin). A total of 23 TDEM and 15 VES soundings were conducted during the years of 2009, 2011 and 2012. In addition, 10 pairs of VES/IDEM soundings were acquired with coincident centers to be able to perform the joint inversion. The joint inversion technique is a promising tool, which enables to get the best of both methods, where the VES add the shallow information and TDEM the deeper one. In this work, the individual and joint inversions were performed using the "Curupira" software. After data process and inversion, the results were interpreted based on geological well information provided by the Department of Water and Electrical Power (DAEE) and the Brazilian Geological Survey (CPRM) which enabled to estimate favorable places to exploitation of water in Bauru and Serra Geral aquifers. For the Baum aquifer, the results suggest areas where thickness exceeds 100 m. In these areas, the resistivity calculated was about 20 Omega center dot m. Therefore, the sediments have been interpreted as saturated sandy clay. In the basalt layer of Serra Geral Formation, the suggested locations present resistivity values <100 Omega center dot m at similar to 200 m depth. The indicated places in sedimentary aquifer and the locations in the fr