查看更多>>摘要：Many problems in applied geophysics can be formulated as a linear inverse problem. The associated problems, however, are large-scale and ill-conditioned. Therefore, regularization techniques are needed to be employed for solving them and generating a stable and acceptable solution. We consider numerical methods for solving such problems in this paper. In order to tackle the ill-conditioning of the problem we use blockiness as a prior information of the subsurface parameters and formulate the problem as a constrained total variation (TV) regularization. The Bregmanized operator splitting (BOS) algorithm as a combination of the Bregman iteration and the proximal forward backward operator splitting method is developed to solve the arranged problem. Two main advantages of this new algorithm are that no matrix inversion is required and that a discrepancy stopping criterion is used to stop the iterations, which allow efficient solution of large-scale problems. The high performance of the proposed TV regularization method is demonstrated using two different experiments: 1) velocity inversion from (synthetic) seismic data which is based on Born approximation, 2) computing interval velocities from RMS velocities via Dix formula. Numerical examples are presented to verify the feasibility of the proposed method for high-resolution velocity inversion. (C) 2018 Elsevier B.V. All rights reserved.
Jeoung, Taek JuKo, Kyung NamKim, Ki SooPark, Kyung-sick...
12页
查看更多>>摘要：Prospective shale plays require a combination of good reservoir and completion qualities. Total organic carbon (TOC) is an important reservoir quality and brittleness is the most critical condition for completion quality. We analyzed seismically-derived brittleness and TOC to investigate the prospectivity of the Horn River Group shale (the Muskwa, Otter Park, Evie shales) of a shale-gas field in the western Horn River Basin, British Columbia, Canada. We used the lambda rho-mu rho brittleness template, constructed from the mineralogy-based brittleness index (MBI) and elastic logs from two wells, to convert the lambda rho and rho rho volumes from prestack seismic inversion to the volume for the brittleness petrotypes (most brittle, intermediate, and least brittle). The probability maps of the most brittle petrotype for the three shales were generated from Bayesian classification, based on the lambda rho-rho rho template. The relationship between TOC and P-wave and S-wave velocity ratio (V-p/V-s) at the wells allowed the conversion of the V-p/V-s volume from prestack inversion to the TOC volume, which in turn was used to construct the TOC maps for the three shales. Increased TOC is correlated with high brittleness, contrasting with the commonly-held understanding. Therefore, the prospectivity of the shales in the study area can be represented by high brittleness and increased TOC. We propose a shale prospectivity index (SPI), computed by the arithmetic average of the normalized probability of the most brittle petrotype and the normalized TOC. The higher SPI corresponds to higher production rates in the Muskwa and Evie shales. The areas of the highest SPI have not been fully tested. The future drilling should be focused on these areas to increase the economic viability of the field. (C) 2018 Elsevier B.V. All rights reserved.
查看更多>>摘要：Multi-frequency measurement of a dispersive electromagnetic (EM) property, such as electrical conductivity, dielectric permittivity, or magnetic permeability, is commonly analyzed for purposes of material characterization. Such an analysis requires inversion of the multi-frequency measurement based on a specific relaxation model, such as Cole-Cole model or Pelton's model. We develop a unified inversion scheme that can be coupled to various type of relaxation models to independently process multi-frequency measurement of varied EM properties for purposes of improved EM-based geomaterial characterization. The proposed inversion scheme is firstly tested in few synthetic cases in which different relaxation models are coupled into the inversion scheme and then applied to multi-frequency complex conductivity, complex resistivity, complex permittivity, and complex impedance measurements. The method estimates up to seven relaxation-model parameters exhibiting convergence and accuracy for random initializations of the relaxation-model parameters within up to 3-orders of magnitude variation around the true parameter values. The proposed inversion method implements a bounded Levenberg algorithm with tuning initial values of damping parameter and its iterative adjustment factor, which are fixed in all the cases shown in this paper and irrespective of the type of measured EM property and the type of relaxation model. Notably, jump-out step and jump-back-in step are implemented as automated methods in the inversion scheme to prevent the inversion from getting trapped around local minima and to honor physical bounds of model parameters. The proposed inversion scheme can be easily used to process various types of EM measurements without major changes to the inversion scheme. (C) 2018 Elsevier B.V. All rights reserved.
查看更多>>摘要：We develop a novel accurate volume integral formula for solving potentials of three dimensional (3D) direct-current resistivity problems with inhomogeneous conductivities in half-space. This new integral formula is composed of the potential, the gradient of Green's function, the gradient of the potential and the anomalous conductivity as the physical variables. First, the unstructured grids are adopted to handle inhomogeneous bodies with complicated shapes. Then, in each anomalous tetrahedron, the potential is represented as its values at vertices and the linear shape functions. Analytical expressions are developed to evaluate singular volume integrals in the final system of linear equations when the observation sites locate in the anomalous body. Finally, two synthetic models are utilized to verify the accuracy and convergence rate of our newly developed volume integral formula and its capability of dealing with complicated models with high conductivity contrasts. (C) 2018 Elsevier B.V. All rights reserved.
查看更多>>摘要：Full-length bonded rock bolts are commonly used in mining, tunneling and slope engineering because of their simple design and resistance to corrosion. However, the length of a rock bolt and grouting quality do not often meet the required design standards in practice because of the concealment and complexity of bolt construction. Non-destructive testing is preferred when testing a rock bolt's quality because of the convenience, low cost and wide detection range. In this paper, a signal analysis method for the non-destructive sound wave testing of full-length bonded rock bolts is presented, which is based on the Hilbert-Huang transform (HHT). First, we introduce the HHT analysis method to calculate the bolt length and identify defect locations based on sound wave reflection test signals, which includes decomposing the test signal via empirical mode decomposition (EMD), selecting the intrinsic mode functions (IMF) using the Pearson Correlation Index (PCI) and calculating the instantaneous phase and frequency via the Hilbert transform (HT). Second, six model tests are conducted using different grouting defects and bolt protruding lengths to verify the effectiveness of the HHT analysis method. Lastly, the influence of the bolt protruding length on the test signal, identification of multiple reflections from defects, bolt end and protruding end, and mode mixing from EMD are discussed. The HHT analysis method can identify the bolt length and grouting defect locations from signals that contain noise at multiple reflected interfaces. The reflection from the long protruding end creates an irregular test signal with many frequency peaks on the spectrum. The reflections from defects barely change the original signal because they are low energy, which cannot be adequately resolved using existing methods. The HHT analysis method can identify reflections from the long protruding end of the bolt and multiple reflections from grouting defects based on mutations in the instantaneous
查看更多>>摘要：We present an electrode layout strategy that allows electrical resistivity profiles to image the third dimension close to the profile plane. This "offset-electrode profile" approach involves laterally displacing electrodes away from the profile line in an alternating fashion and then inverting the resulting data using three-dimensional electrical resistivity tomography software. In our synthetic and field surveys, the offset electrode method succeeds in revealing three-dimensional structures in the vicinity of the profile plane, which we could not achieve using three-dimensional inversions of linear profiles. We confirm and explain the limits of linear electrode profiles through a discussion of the three-dimensional sensitivity patterns: For a homogeneous starting model together with a linear electrode layout, all sensitivities remain symmetric with respect to the profile plane through each inversion step. This limitation can be overcome with offset-electrode layouts by breaking the symmetry pattern among the sensitivities. Thanks to freely available powerful three-dimensional resistivity tomography software and cheap modern computing power, the requirement for full three-dimensional calculations does not create a significant burden and renders the offset-electrode approach a cost-effective method. By offsetting the electrodes in an alternating pattern, as opposed to laying the profile out in a U-shape, we minimize shortening the profile length. (C) 2018 Elsevier B.V. All rights reserved.
Baawain, Mahad S.Al-Futaisi, Ahmed M.Ebrahimi, A.Omidvarborna, Hamid...
9页
查看更多>>摘要：Time Domain Electromagnetic (TDEM) survey as well as drilling investigations were conducted to identify possible contamination of a dumping site in an unsaturated zone located in Barka, Oman. The method was applied to evaluate conductivity of the contaminated plumes in hot and arid/semiarid region, where high temperatures commonly ranged between 35 and 50 degrees C. The drilling investigation was carried out over the survey area to verify the geophysical results. The low-resistivity zone (<80 Omega m), encountered near the subsurface, indicated plume migration caused by liquid waste disposal activities. The combination of TDEM survey results with the lithology of piezometers showed that higher resistivity (>90 Omega m) was correlated with compacted or cemented gravels and cobbles, particularly that of medium dense to very dense gravels and cobbles. Additionally, the TDEM profiles suggested that the plume migration followed a preferential flow path. The resistivity range 40-80 Omega m considered as contaminated areas; however, the drilling results showed the close resistivity domain in the depth >70 m below water table for some profiles (BL1, BL2, BL3, BL4 and BL5). The combined results of drilling wells, piezometers, and TDEM apparent resistivity maps showed a coincidence of the migrated leachate plume and water table. Predicted zone of the probable contamination was located at the depth of around 65 m and horizontal offset ranges 0-280 m, 80-240 m, and 40-85 m in the sounding traverses of BL4, BL6 and BL7, respectively. (C) 2018 Elsevier B.V. All rights reserved.
查看更多>>摘要：We study problems associated with seismic data decomposition and migration imaging. We first represent the seismic data utilizing Gaussian beam basis functions, which have nonzero curvature, and then consider the sparse decomposition technique. The sparse decomposition problem is an l(0)-norm constrained minimization problem. In solving the l(0)-norm minimization, a polynomial Radon transform is performed to achieve sparsity, and a fast gradient descent method is used to calculate the waveform functions. The waveform functions can subsequently be used for sparse Gaussian beam migration. Compared with traditional sparse Gaussian beam methods, the seismic data can be properly reconstructed employing fewer Gaussian beams with nonzero initial curvature. The migration approach described in this paper is more efficient than the traditional sparse Gaussian beam migration. (C) 2018 The Authors. Published by Elsevier B.V.
查看更多>>摘要：Marine gravity anomaly derived from satellite altimetry can be computed using either sea surface height or sea surface slope measurements. Here we consider the slope method and evaluate the errors in the slope of the corrections supplied with the Jason-1 geodetic mission data. The slope corrections are divided into three groups based on whether they are small, comparable, or large with respect to the 1 microradian error in the current sea surface slope models. (1) The small and thus negligible corrections include dry tropospheric correction, inverted barometer correction, solid earth tide and geocentric pole tide. (2) The moderately important corrections include wet tropospheric correction, dual-frequency ionospheric correction and sea state bias. The radiometer measurements are more preferred than model values in the geophysical data records for constraining wet tropospheric effect owing to the highly variable water-vapor structure in atmosphere. The items of dual-frequency ionospheric correction and sea state bias should better not be directly added to range observations for obtaining sea surface slopes since their inherent errors may cause abnormal sea surface slopes and along-track smoothing with uniform distribution weight in certain width is an effective strategy for avoiding introducing extra noises. The slopes calculated from radiometer wet tropospheric corrections, and along-track smoothed dual-frequency ionospheric corrections, sea state bias are generally within +/- 0.5 microradians and no larger than 1 microradians. (3) Ocean tide has the largest influence on obtaining sea surface slopes while most of ocean tide slopes distribute within +/- 3 microradians. Larger ocean tide slopes mostly occur over marginal and island-surrounding seas, and extra tidal models with better precision or with extending process (e.g. Got-e) are strongly recommended for updating corrections in geophysical data records. (C) 2018 Elsevier B.V. All rights reserved.
查看更多>>摘要：A visco-acoustic wave-equation traveltime inversion method is presented that inverts for the shallow subsurface velocity distribution. Similar to the classical wave equation traveltime inversion, this method finds the velocity model that minimizes the squared sum of the traveltime residuals. Even though, wave-equation traveltime inversion can partly avoid the cycle skipping problem, a good initial velocity model is required for the inversion to converge to a reasonable tomogram with different attenuation profiles. When Q model is far away from the real model, the final tomogram is very sensitive to the starting velocity model. Nevertheless, a minor or moderate perturbation of the Q model from the true one does not strongly affect the inversion if the low wavenumber information of the initial velocity model is mostly correct. These claims are validated with numerical tests on both the synthetic and field data sets. (C) 2018 Elsevier B.V. All rights reserved.
NSTL
Elsevier