Nonlinear process control of wave-equation inversion and its application in the detection of gas

Geophysics ◽  
2007 ◽  
Vol 72 (1) ◽  
pp. R9-R18 ◽  
Author(s):  
Yumei Shi ◽  
Wenzhi Zhao ◽  
Hong Cao

The wave equation describes how seismic waves propagate in the subsurface. Inversion methods based on the wave equation naturally take into account the complex behavior of propagating waves and can be used to make accurate estimates of model parameters. However, computational cost and poor convergence have not been overcome, and thus restrict the broad application of this technique. Preconditioned conjugate gradient inversion using back-propagation techniques is a simple, robust implementation of wave-equation inversion in which the step length for correcting the model for each iteration is a crucial factor affecting convergence and hence computational cost. The step length can be calculated by an adaptive controller based on the theory of model reference nonlinear control that ensures that the error energy of the complex system vanishes rapidly. Although the computational cost for each iteration remains the same, the inversion is robust and converges more rapidly than other methods. We tested our method on synthetic data generated from a three-layer fractured model. The inversion converges to the true model after five iterations, and different initial models give similar inversion results. The application to 2D real data from a gas field in western China illustrates that even two iterations yield unambiguous interpretable inversion sections.

Geophysics ◽  
2004 ◽  
Vol 69 (6) ◽  
pp. 1491-1504 ◽  
Author(s):  
W. A. Mulder ◽  
R.‐E. Plessix

Results for wave‐equation migration in the frequency domain using the constant‐density acoustic two‐way wave equation have been compared to images obtained by its one‐way approximation. The two‐way approach produces more accurate reflector amplitudes and provides superior imaging of steep flanks. However, migration with the two‐way wave equation is sensitive to diving waves, leading to low‐frequency artifacts in the images. These can be removed by surgical muting of the input data or iterative migration or high‐pass spatial filtering. The last is the most effective. Iterative migration based on a least‐squares approximation of the seismic data can improve the amplitudes and resolution of the imaged reflectors. Two approaches are considered, one based on the linearized constant‐density acoustic wave equation and one on the full acoustic wave equation with variable density. The first converges quickly. However, with our choice of migration weights and high‐pass spatial filtering for the linearized case, a real‐data migration result shows little improvement after the first iteration. The second, nonlinear iterative migration method is considerably more difficult to apply. A real‐data example shows only marginal improvement over the linearized case. In two dimensions, the computational cost of the two‐way approach has the same order of magnitude as that for the one‐way method. With our implementation, the two‐way method requires about twice the computer time needed for one‐way wave‐equation migration.


Geophysics ◽  
2002 ◽  
Vol 67 (1) ◽  
pp. 38-49 ◽  
Author(s):  
Jorge O. Parra ◽  
Chris L. Hackert ◽  
Anthony W. Gorody ◽  
Valeri Korneev

Guided seismic waves can be used to predict continuity and discontinuity of reservoir structures between wells, with the low‐velocity beds acting as waveguides. We relate guided‐wave signatures to waveguide targets using experimental data acquired at the Stratton gas field in southeast Texas. The observed seismic data indicate the presence of trapped energy in low velocity shale markers between wells 145 and 151. Guided waves in the form of leaky modes are excited, transmitted, and detected in the low‐velocity shale markers at a well separation of 1730 ft (527 m). Dispersion analysis, modeling, frequency–amplitude depth curves, well logs, and lithological information all support the results. Specifically, the characterization of two low‐velocity shale markers, V2 and V5, demonstrates that V2 is more heterogeneous than V5 between the source well 151 and detector well 145. Finally, images of synthetic and real data show the potential applications of the guided‐wave technology as a tool for reservoir characterization.


Author(s):  
Ronan Keane ◽  
H. Oliver Gao

Before a car-following model can be applied in practice, it must first be validated against real data in a process known as calibration. This paper discusses the formulation of calibration as an optimization problem and compares different algorithms for its solution. The optimization consists of an arbitrary car following model, posed as either an ordinary or delay differential equation, being calibrated to an arbitrary source of trajectory data that may include lane changes. Typically, the calibration problem is solved using gradient free optimization. In this work, the gradient of the optimization problem is derived analytically using the adjoint method. The computational cost of the adjoint method does not scale with the number of model parameters, which makes it more efficient than evaluating the gradient numerically using finite differences. Numerical results are presented that show that quasi-Newton algorithms using the adjoint method are significantly faster than a genetic algorithm and also achieve slightly better accuracy of the calibrated model.


2019 ◽  
Vol XVI (2) ◽  
pp. 1-11
Author(s):  
Farrukh Jamal ◽  
Hesham Mohammed Reyad ◽  
Soha Othman Ahmed ◽  
Muhammad Akbar Ali Shah ◽  
Emrah Altun

A new three-parameter continuous model called the exponentiated half-logistic Lomax distribution is introduced in this paper. Basic mathematical properties for the proposed model were investigated which include raw and incomplete moments, skewness, kurtosis, generating functions, Rényi entropy, Lorenz, Bonferroni and Zenga curves, probability weighted moment, stress strength model, order statistics, and record statistics. The model parameters were estimated by using the maximum likelihood criterion and the behaviours of these estimates were examined by conducting a simulation study. The applicability of the new model is illustrated by applying it on a real data set.


Sensors ◽  
2021 ◽  
Vol 21 (8) ◽  
pp. 2704
Author(s):  
Yunhan Lin ◽  
Wenlong Ji ◽  
Haowei He ◽  
Yaojie Chen

In this paper, an intelligent water shooting robot system for situations of carrier shake and target movement is designed, which uses a 2 DOF (degree of freedom) robot as an actuator, a photoelectric camera to detect and track the desired target, and a gyroscope to keep the robot’s body stable when it is mounted on the motion carriers. Particularly, for the accurate shooting of the designed system, an online tuning model of the water jet landing point based on the back-propagation algorithm was proposed. The model has two stages. In the first stage, the polyfit function of Matlab is used to fit a model that satisfies the law of jet motion in ideal conditions without interference. In the second stage, the model uses the back-propagation algorithm to update the parameters online according to the visual feedback of the landing point position. The model established by this method can dynamically eliminate the interference of external factors and realize precise on-target shooting. The simulation results show that the model can dynamically adjust the parameters according to the state relationship between the landing point and the desired target, which keeps the predicted pitch angle error within 0.1°. In the test on the actual platform, when the landing point is 0.5 m away from the position of the desired target, the model only needs 0.3 s to adjust the water jet to hit the target. Compared to the state-of-the-art method, GA-BP (genetic algorithm-back-propagation), the proposed method’s predicted pitch angle error is within 0.1 degree with 1/4 model parameters, while costing 1/7 forward propagation time and 1/200 back-propagation calculation time.


Mathematics ◽  
2021 ◽  
Vol 9 (16) ◽  
pp. 1850
Author(s):  
Rashad A. R. Bantan ◽  
Farrukh Jamal ◽  
Christophe Chesneau ◽  
Mohammed Elgarhy

Unit distributions are commonly used in probability and statistics to describe useful quantities with values between 0 and 1, such as proportions, probabilities, and percentages. Some unit distributions are defined in a natural analytical manner, and the others are derived through the transformation of an existing distribution defined in a greater domain. In this article, we introduce the unit gamma/Gompertz distribution, founded on the inverse-exponential scheme and the gamma/Gompertz distribution. The gamma/Gompertz distribution is known to be a very flexible three-parameter lifetime distribution, and we aim to transpose this flexibility to the unit interval. First, we check this aspect with the analytical behavior of the primary functions. It is shown that the probability density function can be increasing, decreasing, “increasing-decreasing” and “decreasing-increasing”, with pliant asymmetric properties. On the other hand, the hazard rate function has monotonically increasing, decreasing, or constant shapes. We complete the theoretical part with some propositions on stochastic ordering, moments, quantiles, and the reliability coefficient. Practically, to estimate the model parameters from unit data, the maximum likelihood method is used. We present some simulation results to evaluate this method. Two applications using real data sets, one on trade shares and the other on flood levels, demonstrate the importance of the new model when compared to other unit models.


2021 ◽  
Vol 11 (15) ◽  
pp. 6998
Author(s):  
Qiuying Li ◽  
Hoang Pham

Many NHPP software reliability growth models (SRGMs) have been proposed to assess software reliability during the past 40 years, but most of them have focused on modeling the fault detection process (FDP) in two ways: one is to ignore the fault correction process (FCP), i.e., faults are assumed to be instantaneously removed after the failure caused by the faults is detected. However, in real software development, it is not always reliable as fault removal usually needs time, i.e., the faults causing failures cannot always be removed at once and the detected failures will become more and more difficult to correct as testing progresses. Another way to model the fault correction process is to consider the time delay between the fault detection and fault correction. The time delay has been assumed to be constant and function dependent on time or random variables following some kind of distribution. In this paper, some useful approaches to the modeling of dual fault detection and correction processes are discussed. The dependencies between fault amounts of dual processes are considered instead of fault correction time-delay. A model aiming to integrate fault-detection processes and fault-correction processes, along with the incorporation of a fault introduction rate and testing coverage rate into the software reliability evaluation is proposed. The model parameters are estimated using the Least Squares Estimation (LSE) method. The descriptive and predictive performance of this proposed model and other existing NHPP SRGMs are investigated by using three real data-sets based on four criteria, respectively. The results show that the new model can be significantly effective in yielding better reliability estimation and prediction.


2019 ◽  
Vol 7 (2) ◽  
pp. T255-T263 ◽  
Author(s):  
Yanli Liu ◽  
Zhenchun Li ◽  
Guoquan Yang ◽  
Qiang Liu

The quality factor ([Formula: see text]) is an important parameter for measuring the attenuation of seismic waves. Reliable [Formula: see text] estimation and stable inverse [Formula: see text] filtering are expected to improve the resolution of seismic data and deep-layer energy. Many methods of estimating [Formula: see text] are based on an individual wavelet. However, it is difficult to extract the individual wavelet precisely from seismic reflection data. To avoid this problem, we have developed a method of directly estimating [Formula: see text] from reflection data. The core of the methodology is selecting the peak-frequency points to linear fit their logarithmic spectrum and time-frequency product. Then, we calculated [Formula: see text] according to the relationship between [Formula: see text] and the optimized slope. First, to get the peak frequency points at different times, we use the generalized S transform to produce the 2D high-precision time-frequency spectrum. According to the seismic wave attenuation mechanism, the logarithmic spectrum attenuates linearly with the product of frequency and time. Thus, the second step of the method is transforming a 2D spectrum into 1D by variable substitution. In the process of transformation, we only selected the peak frequency points to participate in the fitting process, which can reduce the impact of the interference on the spectrum. Third, we obtain the optimized slope by least-squares fitting. To demonstrate the reliability of our method, we applied it to a constant [Formula: see text] model and the real data of a work area. For the real data, we calculated the [Formula: see text] curve of the seismic trace near a well and we get the high-resolution section by using stable inverse [Formula: see text] filtering. The model and real data indicate that our method is effective and reliable for estimating the [Formula: see text] value.


2020 ◽  
Vol 26 ◽  
pp. 7
Author(s):  
Hui Wei ◽  
Shuguan Ji

This paper is devoted to the study of periodic (in time) solutions to an one-dimensional semilinear wave equation with x-dependent coefficients under various homogeneous boundary conditions. Such a model arises from the forced vibrations of a nonhomogeneous string and propagation of seismic waves in nonisotropic media. By combining variational methods with an approximation argument, we prove that there exist infinitely many periodic solutions whenever the period is a rational multiple of the length of the spatial interval. The proof is essentially based on the spectral properties of the wave operator with x-dependent coefficients.


2013 ◽  
Vol 2013 ◽  
pp. 1-9 ◽  
Author(s):  
Ali Kargarnejad ◽  
Mohsen Taherbaneh ◽  
Amir Hosein Kashefi

Tracking maximum power point of a solar panel is of interest in most of photovoltaic applications. Solar panel modeling is also very interesting exclusively based on manufacturers data. Knowing that the manufacturers generally give the electrical specifications of their products at one operating condition, there are so many cases in which the specifications in other conditions are of interest. In this research, a comprehensive one-diode model for a solar panel with maximum obtainable accuracy is fully developed only based on datasheet values. The model parameters dependencies on environmental conditions are taken into consideration as much as possible. Comparison between real data and simulations results shows that the proposed model has maximum obtainable accuracy. Then a new fuzzy-based controller to track the maximum power point of the solar panel is also proposed which has better response from speed, accuracy and stability point of view respect to the previous common developed one.


Sign in / Sign up

Export Citation Format

Share Document