查看更多>>摘要：The attenuation of random noise is important for improving the signal to noise ratio (SNR). However, the precondition for most conventional denoising methods is that the noisy data must be sampled on a uniform grid, making the conventional methods unsuitable for non-uniformly sampled data. In this paper, a denoising method capable of regularizing the noisy data from a non-uniform grid to a specified uniform grid is proposed. Firstly, the denoising method is performed for every time slice extracted from the 3D noisy data along the source and receiver directions, then the 2D non-equispaced fast Fourier transform (NFFT) is introduced in the conventional fast discrete curvelet transform (FDCT). The non-equispaced fast discrete curvelet transform (NFDCT) can be achieved based on the regularized inversion of an operator that links the uniformly sampled curvelet coefficients to the non-uniformly sampled noisy data. The uniform curvelet coefficients can be calculated by using the inversion algorithm of the spectral projected-gradient for l(1)-norm problems. Then local threshold factors are chosen for the uniform curvelet coefficients for each decomposition scale, and effective curvelet coefficients are obtained respectively for each scale. Finally, the conventional inverse FDCT is applied to the effective curvelet coefficients. This completes the proposed 3D denoising method using the non-equispaced curvelet transform in the source receiver domain. The examples for synthetic data and real data reveal the effectiveness of the proposed approach in applications to noise attenuation for non-uniformly sampled data compared with the conventional FDCT method and wavelet transformation. (C) 2018 Elsevier B.V. All rights reserved.
查看更多>>摘要：The location of buried utility pipes is often unknown. We use the time-domain induced polarization method to non-intrusively localize metallic pipes. A new approach, based on injecting a primary electrical current between a pair of electrodes and measuring the time-lapse voltage response on a set of potential electrodes after shutting down this primary current is used. The secondary voltage is measured on all the electrodes with respect to a single electrode used as a reference for the electrical potential, in a way similar to a self-potential time lapse survey. This secondary voltage is due to the formation of a secondary current density in the ground associated with the polarization of the metallic pipes. An algorithm is designed to localize the metallic object using the secondary voltage distribution by performing a tomography of the secondary source current density associated with the polarization of the pipes. This algorithm is first benchmarked on a synthetic case. Then, two laboratory sandbox experiments are performed with buried metallic pipes located in a sandbox filled with some clean sand. In Experiment #1, we use a horizontal copper pipe while in Experiment #2 we use an inclined stainless steel pipe. The result shows that the method is effective in localizing these two pipes. At the opposite, electrical resistivity tomography is not effective in localizing the pipes because they may appear resistive at low frequencies. This is due to the polarization of the metallic pipes which blocks the charge carriers at its external boundaries. (C) 2018 Elsevier B.V. All rights reserved.
查看更多>>摘要：Ground roll is a type of coherent noise in land seismic data that has low frequency, low velocity and high amplitude. It damages reflection events that contain important information about subsurface structures, hence the removal of ground roll is a crucial step in seismic data processing. A suitable transform is needed for removal of ground roll. Curvelet transform is an effective sparse transform that optimally represents seismic events. In addition, the curvelets can provide a multiscale and multidirectional decomposition of the input data in time-frequency and angular domain, which can help distinguish between ground roll and useful signals. In this paper, we apply synchrosqueezed curvelet transform (SSCT) for ground roll attenuation. The synchrosqueezing technique in SSCT is used to precisely reallocate the energy of local wave vectors in order to separate ground roll from the original data with higher resolution and higher fidelity. Examples of synthetic and field seismic data reveal that SSCT performs well in the suppression of aliased and non-aliased ground roll while preserving reflection waves, in comparison with high-pass filtering, wavelet and curvelet methods. (C) 2018 Elsevier B.V. All rights reserved.
Motieyan-Najar, M. E.Sadeghi, J.Zakeri, J. A.Yousefi, B....
11页
查看更多>>摘要：Ballast plays an important role in the stability of railway track systems. The effectiveness of the ballast in maintaining the track stability is very much dependent on its mechanical conditions. The available ballast maintenance approaches are mainly based on only track geometry conditions (such as track profile) which do not sufficiently reflect the ballast mechanical behaviors. That is, the ballast potential of degradation (i.e., ballast long term behaviors) has been omitted. This makes the effectiveness of the current ballast maintenance approach questionable, indicating a need for a more comprehensive and effective ballast conditions assessment technique. In response to this need, two ballast condition indices based on ballast geometry degradation (BGI) and the level of ballast fouling (BFI) as the main indicators of ballast mechanical behavior were developed. The BC! is a function of the standard deviations of track alignment, unevenness and twist. The BFI was developed based on the data obtained from the ground penetration radar (GPR). Making use of the new indices, a more reliable maintenance algorithm was developed. Through illustrations of the applicability of the new maintenance algorithm in a railway line, it was shown that the new algorithm causes a considerable improvement in the maintenance effectiveness and an increase in the life cycle of railway tracks by making more effective allocation of resources and more accurate maintenance planning. (C) 2018 Elsevier B.V. All rights reserved.
查看更多>>摘要：We propose an iterative three-dimensional (3D) deblending scheme using 3D multi-scale shaping operator to separate 3D simultaneous source data. The proposed scheme is based on the property that signal is coherent, whereas interference is incoherent in some domains, e.g., common receiver domain and common midpoint domain. In two-dimensional (2D) blended record, the coherency difference of signal and interference is in only one spatial direction. Compared with 2D deblending, the 3D deblending can take more sparse constraints into consideration to obtain better performance, e.g., in 3D common receiver gather, the coherency difference is in two spatial directions. Furthermore, with different levels of coherency, signal and interference distribute in different scale curvelet domains. In both 2D and 3D blended records, most coherent signal locates in coarse scale curvelet domain, while most incoherent interference distributes in fine scale curvelet domain. The scale difference is larger in 3D deblending, thus, we apply the multi-scale shaping scheme to further improve the 3D deblending performance. We evaluate the performance of 3D and 2D deblending with the multi-scale and global shaping operators, respectively. One synthetic and one field data examples demonstrate the advantage of the 3D deblending with 3D multi-scale shaping operator. (C) 2018 Elsevier B.V. All rights reserved.
查看更多>>摘要：A novel technique is developed to level airborne geophysical data using principal component analysis based on flight line difference. In the paper, flight line difference is introduced to enhance the features of levelling error for airborne electromagnetic (AEM) data and improve the correlation between pseudo tie lines. Thus we conduct levelling to the flight line difference data instead of to the original AEM data directly. Pseudo tie lines are selected distributively cross profile direction, avoiding the anomalous regions. Since the levelling errors of selective pseudo tie lines show high correlations, principal component analysis is applied to extract the local levelling errors by low-order principal components reconstruction. Furthermore, we can obtain the levelling errors of original AEM data through inverse difference after spatial interpolation. This levelling method does not need to fly tie lines and design the levelling fitting function. The effectiveness of this method is demonstrated by the levelling results of survey data, comparing with the results from tie-line levelling and flight-line correlation levelling. (C) 2018 Elsevier B.V. All rights reserved.
Gholamy, Shaban AliBasir, Hadi MandaviJavaherian, AbdolrahimShomali, Zaher Hossein...
11页
查看更多>>摘要：Imaging is a key step in seismic data processing. To date, a myriad of advanced pre-stack depth migration approaches have been developed; however, reverse time migration (RTM) is still considered as the high-end imaging algorithm. The main limitations associated with the performance cost of reverse time migration are the intensive computation of the forward and backward simulations, time consumption, and memory allocation related to imaging condition. Based on the reduced order modeling, we proposed an algorithm, which can be adapted to all the aforementioned factors. Our proposed method benefit from Krylov subspaces method to compute certain mode shapes of the velocity model computed by as an orthogonal base of reduced order modeling. Reverse time migration by reduced order modeling is helpful concerning the highly parallel computation and strongly reduces the memory requirement of reverse time migration. The synthetic model results showed that suggested method can decrease the computational costs of reverse time migration by several orders of magnitudes, compared with reverse time migration by finite element method. (C) 2018 Elsevier B.V. All rights reserved.
查看更多>>摘要：With the development of Logging While Drilling (LWD) technology, dual-induction LWD logging is not only widely applied in deviated wells and horizontal wells, but it is used commonly in vertical wells. Accordingly, it is necessary to simulate the response of LWD tools in vertical wells for logging interpretation. In this paper, the investigation characteristics, the effects of the tool structure, skin effect and drilling environment of a dual induction LWD tool are simulated by the three-dimensional (3D) finite element method (FEM). In order to closely simulate the actual situation, real structure of the tool is taking into account. The results demonstrate that the influence of the background value of the tool structure can be eliminated. The values of deducting the background of a tool structure and analytical solution have a quantitative agreement in homogeneous formations. The effect of measurement frequency could be effectively eliminated by chart of skin effect correction. In addition, the measurement environment, borehole size, mud resistivity, shoulder bed, layer thickness and invasion, have an effect on the true resistivity. To eliminate these effects, borehole correction charts, shoulder bed correction charts and tornado charts are computed based on real tool structure. Based on correction charts, well logging data can be corrected automatically by a suitable interpolation method, which is convenient and fast. Verified with actual logging data in vertical wells, this method could obtain the true resistivity of formation. (C) 2018 Published by Elsevier B.V.
Ranjan, ShashiKambhammettu, B. V. N. P.Peddinti, Srinivasa RaoAdinarayana, J....
10页
查看更多>>摘要：Image reconstruction from discrete electrical responses pose a number of computational and mathematical challenges. Application of smoothness constrained regularized inversion from limited measurements may fail to detect resistivity anomalies and sharp interfaces separated by hydro stratigraphic units. Under favourable conditions, compressed sensing (CS) can be thought of an alternative to reconstruct the image features by finding sparse solutions to highly underdetermined linear systems. This paper deals with the development of a CS assisted, 3-D resistivity inversion algorithm for use with hydrogeologists and groundwater scientists. CS based l(1)-regularized least square algorithm was applied to solve the resistivity inversion problem. Sparseness in the model update vector is introduced through block oriented discrete cosine transformation, with recovery of the signal achieved through convex optimization. The equivalent quadratic program was solved using primal-dual interior point method. Applicability of the proposed algorithm was demonstrated using synthetic and field examples drawn from hydrogeology. The proposed algorithm has outperformed the conventional (smoothness constrained) least square method in recovering the model parameters with much fewer data, yet preserving the sharp resistivity fronts separated by geologic layers. Resistivity anomalies represented by discrete homogeneous blocks embedded in contrasting geologic layers were better imaged using the proposed algorithm. In comparison to conventional algorithm, CS has resulted in an efficient (an increase in R-2 from 0.62 to 0.78; a decrease in RMSE from 125.14 Omega-m to 72.46 Omega-m), reliable, and fast converging (run time decreased by about 25%) solution. (C) 2018 Elsevier B.V. All rights reserved.
查看更多>>摘要：We propose an algorithm for calculating all-time apparent resistivity from transient electromagnetic induction logging. The algorithm is based on the whole-space transient electric field expression of the uniform model and Halley's optimisation. In trial calculations for uniform models, the all-time algorithm is shown to have high accuracy. We use the finite-difference time-domain method to simulate the transient electromagnetic field in radial two-layer models without wall rock and convert the simulation results to apparent resistivity using the all-time algorithm. The time-varying apparent resistivity reflects the radially layered geoelectrical structure of the models and the apparent resistivity of the earliest time channel follows the true resistivity of the inner layer; however, the apparent resistivity at larger times reflects the comprehensive electrical characteristics of the inner and outer layers. To accurately identify the outer layer resistivity based on the series relationship model of the layered resistance, the apparent resistivity and diffusion depth of the different time channels are approximately replaced by related model parameters; that is, we propose an apparent resistivity correction algorithm. By correcting the time-varying apparent resistivity of radial two-layer models, we show that the correction results reflect the radially layered electrical structure and the corrected resistivities of the larger time channels follow the outer layer resistivity. The transient electromagnetic fields of radially layered models with wall rock are simulated to obtain the 2D time-varying profiles of the apparent resistivity and corrections. The results suggest that the time-varying apparent resistivity and correction results reflect the vertical and radial geoelectrical structures. For models with small wall-rock effect, the correction removes the effect of the low resistance inner layer on the apparent resistivity of the larger time channels. (C) 2018 Elsevie
NSTL
Elsevier