Reconstruction of borehole radar images by a modified f–k migration

2020 ◽  
Vol 221 (3) ◽  
pp. 1626-1634
Author(s):  
Jianjian Huo ◽  
Binzhong Zhou ◽  
Qing Zhao ◽  
Iain M Mason

SUMMARY Borehole radar (BHR) is an effective imaging tool. It can be used to detect and map faults, fractures, folds, domes, partings and mine workings. Most BHRs have azimuthally omnidirectional radiation patterns. The echoes sensed by such BHRs may come from any direction. Considering a radar in a straight borehole that passes through a stack of flat reflection planes, V-shaped events or crosses appear on the time section. One of the arms of each cross is a real image while the other is an ambiguity of known origin. Directional ambiguities such as these obstruct efforts to interpret the data. In this paper, we address this difficulty by using a modified f–k migration algorithm to translate crosses into lines on the final section that are consistent with a priori information about for example bedding. Compared with conventional strategies, for example migration + f–k dip filter, this approach integrates the two separated processes into one and is straightforward, computationally effective and simple to implement. The method is demonstrated using a synthetic model and a real BHR field data set. It allows the interpreter to use a priori information about fault swarms or plausible bedding planes at an early stage. The reconstructed BHR image helps the search for geological anomalies such as fractures, partings, domes and rolls that could be a hazard for mining.

Geophysics ◽  
2019 ◽  
Vol 84 (3) ◽  
pp. E125-E141 ◽  
Author(s):  
Francesca Pace ◽  
Alessandro Santilano ◽  
Alberto Godio

We implement the particle swarm optimization (PSO) algorithm for the two-dimensional (2D) magnetotelluric (MT) inverse problem. We first validate PSO on two synthetic models of different complexity and then apply it to an MT benchmark for real-field data, the COPROD2 data set (Canada). We pay particular attention to the selection of the PSO input parameters to properly address the complexity of the 2D MT inverse problem. We enhance the stability and convergence of the solution of the geophysical problem by applying the hierarchical PSO with time-varying acceleration coefficients (HPSO-TVAC). Moreover, we parallelize the code to reduce the computation time because PSO is a computationally demanding global search algorithm. The inverse problem was solved for the synthetic data both by giving a priori information at the beginning and by using a random initialization. The a priori information was given to a small number of particles as the initial position within the search space of solutions, so that the swarming behavior was only slightly influenced. We have demonstrated that there is no need for the a priori initialization to obtain robust 2D models because the results are largely comparable with the results from randomly initialized PSO. The optimization of the COPROD2 data set provides a resistivity model of the earth in line with results from previous interpretations. Our results suggest that the 2D MT inverse problem can be successfully addressed by means of computational swarm intelligence.


Author(s):  
EVGENIA DIMITRIADOU ◽  
ANDREAS WEINGESSEL ◽  
KURT HORNIK

In this paper we present a voting scheme for fuzzy cluster algorithms. This voting method allows us to combine several runs of cluster algorithms resulting in a common partition. This helps us to tackle the problem of choosing the appropriate clustering method for a data set where we have no a priori information about it. We mathematically derive the algorithm from theoretical considerations. Experiments show that the voting algorithm finds structurally stable results. Several cluster validity indexes show the improvement of the voting result in comparison to simple fuzzy voting.


2018 ◽  
pp. 36-43
Author(s):  
Володимир Олексійович Габрінець ◽  
Сергій Михайлович Подольчак

The subject matter of the article is an analysis of methods for assessing the reliability of a complex technical system applicable to a rocket engine for solving the problem of estimating the reliability of a rocket engine at an early stage of its development. The goal is to determine the method of estimating the reliability of a rocket engine with insufficient statistical information. The tasks to be solved: to carry out the analysis of existing methods for assessing the reliability of a complex technical system; to justify the expediency of their application for assessing the reliability of the rocket engine at an early stage of its development; choose an effective method for solving the problem. In work the analysis of following methods was carried out: a method of rationing reliability taking into account the complexity of the system and the importance of the elements belonging to the system; method of "weighting factors"; method of accounting for a priori information when using the model "parameter – one-way tolerance field". The following results obtained during the research. As a result of the analysis of existing methods for assessing the reliability of a complex technical system, it argues that they can apply to a rocket engine in terms of the availability of sufficient information When determining the method for estimating the reliability of a rocket engine with insufficient statistical information, a reliability chosen model that allows estimating the level of reliability of a rocket engine at an early stage of its development with insufficient information.  Conclusions. The scientific novelty of the results is following: the proposed method makes it possible to assess the level of reliability of rocket engines of new developments when the designers do not have sufficient statistical information on the results of reliability tests, that is, at the very beginning of the development of the rocket engine; developed methods for assessing reliability of a complex technical system applied to a rocket engine with the use of the "parameter – one-way tolerance field" model was further developed, thanks to the availability of a priori information, which makes it possible to use them for various approaches to assess the reliability of a rocket engine, depending on the amount of information available


2015 ◽  
Vol 8 (3) ◽  
pp. 3399-3422 ◽  
Author(s):  
E. Maillard Barras ◽  
A. Haefele ◽  
R. Stübi ◽  
D. Ruffieux

Abstract. We present a method to derive the site atmospheric state best estimate (SASBE) of the ozone profile combining brightness temperature spectra around the 142 GHz absorption line of ozone measured by the microwave radiometer SOMORA and ozone profiles measured by the radiosonde (RS). The SASBE ozone profile is obtained using the radiosonde ozone profile as a priori information in an optimal estimation retrieval of the SOMORA radiometer. The resulting ozone profile ranges from ground up to 65 km altitude and makes optimal use of the available information at each altitude. The high vertical resolution of the radiosonde profile can be conserved and the uncertainty of the SASBE is well defined at each altitude. A SASBE ozone profile dataset has been generated for Payerne, Switzerland, with a temporal resolution of 3 profiles a week for the time period ranging from 2011 to 2013. The relative difference of the SASBE ozone profiles to the AURA/MLS ozone profiles lies between −3 to 6% over the vertical range of 20–65 km. Above 20 km, the agreement between the SASBE and AURA/MLS ozone profiles is better than the agreement between the operational SOMORA ozone data set and AURA/MLS. Below 20 km the SASBE ozone data are identical to the radiosonde measurements. The same method has been applied to ECWMF-ERA interim ozone profiles and SOMORA data to generate a SASBE dataset with a time resolution of 4 profiles per day. These SASBE ozone profiles agree between −4 and +8% with AURA/MLS. The improved agreement of the SASBE datasets with AURA/MLS above 20 km demonstrates the benefit of better a priori information in the retrieval of ozone from brightness temperature data.


Geophysics ◽  
2019 ◽  
Vol 84 (6) ◽  
pp. E387-E402 ◽  
Author(s):  
Max A. Meju ◽  
Randall L. Mackie ◽  
Federico Miorelli ◽  
Ahmad Shahir Saleh ◽  
Roger V. Miller

Geologic interpretation of 3D anisotropic resistivity models from conventional marine controlled-source electromagnetic (CSEM) data inversion faces difficulties in low-resistivity contrast sediments and structurally complex environments that typify the new frontiers for hydrocarbon exploration. Currently, the typically reconstructed horizontal resistivity [Formula: see text] and vertical resistivity [Formula: see text] models often have conflicting depth structures that are difficult to explain in terms of subsurface geology, and the resulting resistivities may not be close to the true formation resistivities required for estimating reservoir parameters. We have investigated the concept that an objective geologically oriented or structurally tailored inversion can be achieved by requiring that the cross-product of the gradient of horizontal resistivity and the gradient of the vertical resistivity is equal to zero at significant geologic boundaries. We incorporate this boundary-shape criterion in our 3D inverse problem formulations, implemented within nonlinear model-space and conjugate-gradient contexts, for cases in which a priori calibration data from wells and/or seismically derived subsurface boundaries are available and for cases in which these are lacking. The resulting fit-for-purpose solutions serve to better analyze the peculiarity of a given data set. We applied these algorithms to synthetic and field CSEM data sets representing a fold-thrust environment with low-resistivity and low-contrast sediments. The resulting [Formula: see text] and [Formula: see text] models from cross-gradient joint inversion of synthetic data of appropriate frequency bandwidth without a priori information are structurally similar and consistent with the test models, whereas those from the inversions of band-limited field data are consistent with the available seismic and resistivity well-log data. This particular approach will thus be useful for lithologic correlation in frontier regions with limited a priori information using broadband CSEM data. For these band-limited field data, we found that the anisotropic bulk resistivities of the low-contrast sediments are better determined by incorporating a priori calibration data from triaxial resistivity logs and seismic horizons.


1979 ◽  
Vol 69 (2) ◽  
pp. 387-396 ◽  
Author(s):  
James A. Hileman

abstract Phase times for an ensemble of 20 aftershocks of the June 1, 1975, Galway Lake earthquake were inverted to determine hypocenters and local velocity structure simultaneously. A maximum likelihood formulation of linear, least-squares inversion was used so that the data were weighted according to their estimated variances. A trade-off parameter controlled the relative importance of the RMS error and the amount by which the model changed at each iteration. Individual aftershocks must be selected so that a wide variety of travel paths are represented. Initial trials disclosed that careful constraints are necessary for allowed changes in station delays. Fixed velocity boundaries, based on a priori information, were used for each of several starting models, and only layer velocities were allowed to vary. Each of the trials indicated that shallow crustal velocities in the vincinity of Galway Lake are somewhat lower than those of the usual velocity models. The velocities are not strongly constrained by this data set, but the results are consistent with a subsequent, detailed refraction survey by other workers. Mapping the travel-time surface in time-depth-distance space helps clarify the limitations inherent in a given problem.


Geophysics ◽  
1996 ◽  
Vol 61 (2) ◽  
pp. 394-408 ◽  
Author(s):  
Yaoguo Li ◽  
Douglas W. Oldenburg

We present a method for inverting surface magnetic data to recover 3-D susceptibility models. To allow the maximum flexibility for the model to represent geologically realistic structures, we discretize the 3-D model region into a set of rectangular cells, each having a constant susceptibility. The number of cells is generally far greater than the number of the data available, and thus we solve an underdetermined problem. Solutions are obtained by minimizing a global objective function composed of the model objective function and data misfit. The algorithm can incorporate a priori information into the model objective function by using one or more appropriate weighting functions. The model for inversion can be either susceptibility or its logarithm. If susceptibility is chosen, a positivity constraint is imposed to reduce the nonuniqueness and to maintain physical realizability. Our algorithm assumes that there is no remanent magnetization and that the magnetic data are produced by induced magnetization only. All minimizations are carried out with a subspace approach where only a small number of search vectors is used at each iteration. This obviates the need to solve a large system of equations directly, and hence earth models with many cells can be solved on a deskside workstation. The algorithm is tested on synthetic examples and on a field data set.


Sign in / Sign up

Export Citation Format

Share Document