Reflection full-waveform inversion using a modified phase misfit function

Reflection full-waveform inversion （RFWI） updates the low- and high- wavenumber components, and yields more accurate initial models compared with conventional full-waveform inversion （FWI）. However, there is strong nonlinearity in conventional RFWI because of the lack of low-frequency data and the complexity of the amplitude. The separation of phase and amplitude information makes RFWI more linear. Traditional phase-calculation methods face severe phase wrapping. To solve this problem, we propose a modified phase-calculation method that uses the phase-envelope data to obtain the pseudo phase information. Then, we establish a pseudophase-information-based objective function for RFWI, with the corresponding source and gradient terms. Numerical tests verify that the proposed calculation method using the phase-envelope data guarantees the stability and accuracy of the phase information and the convergence of the objective function. The application on a portion of the Sigsbee2A model and comparison with inversion results of the improved RFWI and conventional FWI methods verify that the pseudophase-based RFWI produces a highly accurate and efficient velocity model. Moreover, the proposed method is robust to noise and high frequency.

Reflection-based traveltime and waveform inversion with second-order optimization

Reflection-based inversion that aims to reconstruct the low-to-intermediate wavenumbers of the subsurface model, can be a complementary to refraction-data-driven full-waveform inversion(FWI), especially for the deep target area where diving waves cannot be acquired at the surface. Nevertheless, as a typical nonlinear inverse problem, reflection waveform inversion may easily suffer from the cycleskipping issue and have a slow convergence rate, if gradient-based first-order optimization methods are used. To improve the accuracy and convergence rate, we introduce the Hessian operator into reflection traveltime inversion(RTI) and reflection waveform inversion(RWI) in the framework of second-order optimization. A practical two-stage workflow is proposed to build the velocity model, in which Gauss-Newton RTI is first applied to mitigate the cycle-skipping problem and then Gauss-Newton RWI is employed to enhance the model resolution. To make the Gauss-Newton iterations more efficiently and robustly for large-scale applications, we introduce proper preconditioning for the Hessian matrix and design appropriate strategies to reduce the computational costs. The example of a real dataset from East China Sea demonstrates that the cascaded Hessian-based RTI and RWI have good potential to improve velocity model building and seismic imaging, especially for the deep targets.

Research on Thin Layer Structure Identification and Sedimentary Facies of Middle and Deep Layers Based on Reflection Coefficient Inversion—By Taking Dongying Formation of CFD Oilfield in Bohai Offshore as an Example

The sand layer B of Dongying Formation of CFD oilfield in Bohai offshore belongs to the middle deep layer of buried hill overlap deposit. Its reservoir distribution has the characteristics of large burial depth, thin thickness and rapidly lateral change. Because of low resolution of seismic data and overlying sand layer. It is difficult to identify and interpret the structure of sand layer accurately. The uncertainty of structure and reservoir restricts the fine development of B sand layer. In order to identify the top surface of reservoir effectively. The seismic data are processed by using the reflection coefficient inversion method. The results show that the inversion resolution of reflection coefficient is significantly higher than that of original data. The top surface of sand layer B and its overlying sand layer can be well identified and traced. Carrying out structural interpretation of B sand layer based on reflection coefficient inversion data and the microstructure and the formation tip extinction point are implemented. Based on the constraint of new interpretation level, the sedimentary facies plane distribution of B sand layer is described and make prediction of dominant reservoir development area in detail combining with sedimentary paleogeomorphology, along layer attribute section and limited drilling data. The research shows that the study area is mainly from the northwest material sources, the slope belt in the northwest is close to the lake shoreline with a gentle slope and shallow water depositional environment, which is located on the main transport and deposition channels. The shallow water gentle slope landform is suitable for forming large-area sand bar deposition, mainly composed of underwater distributary channel and debouch bars facies, which is the dominant reservoir development area. The research conclusion guides the deployment and implementation of the development well location effectively.

Reflection full waveform inversion

Because of the combination of optimization algorithms and full wave equations, full-waveform inversion(FWI) has become the frontier of the study of seismic exploration and is gradually becoming one of the essential tools for obtaining the Earth interior information. However, the application of conventional FWI to pure reflection data in the absence of a highly accurate starting velocity model is difficult. Compared to other types of seismic waves, reflections carry the information of the deep part of the subsurface. Reflection FWI, therefore, is able to improve the accuracy of imaging the Earth interior further. Here, we demonstrate a means of achieving this successfully by interleaving least-squares RTM with a version of reflection FWI in which the tomographic gradient that is required to update the background macro-model is separated from the reflectivity gradient using the Born approximation during forward modeling. This provides a good update to the macro-model. This approach is then followed by conventional FWI to obtain a final high-fidelity high-resolution result from a poor starting model using only reflection data.Further analysis reveals the high-resolution result is achieved due to a deconvolution imaging condition implicitly used by FWI.

Geoacoustic inversion based on reflection model of effective density fluid approximation

In order to obtain the physical and geoacoustic properties of marine sediments,an inverse method using reflection loss of different grazing angles is presented.The reflection loss is calculated according to the reflection model of effective density fluid approximation.A two-step hybrid optimization algorithm combining differential evolution and particle swarm optimization along with Bayesian inversion is employed in estimation of porosity,mean grain size,mass density and bulk modulus of grains.Based on the above physical parameters,geoacoustic parameters,including sound speed and attenuation,are further calculated.According to the numerical simulations,we can draw a conclusion that all the parameters can be well estimated with the exception of bulk modulus of grains.In particular,this indirect inverse method for bottom geoacoustic parameters performs high accuracy and strong robustness.The relative errors are 0.092%and 17%,respectively.Finally,measured reflection loss data of sandy sediments at the bottom of a water tank is analyzed,and the estimation value,uncertainty and correlation of each parameter are presented.The availability of this inverse method is verified through comparison between inverse results and part of measured parameters.

AN INVERSION METHOD FOR OBTAINING BOTTOM REFLECTION LOSS

A method for the inversion of the transmission losses for the bottom reflection loss is proposed on the basis of the theory of the smooth- averaged sound field. The procedure of the inversion is based on the criterion of the least square error in the transmission losses between calculation and measurement. By using the Gauss - Newton iterative approach, the non - linear least square aloqrithm is equivalent to solving a sequence of lineared least square problems. The physical causes of the instability of the inversion problem are discussed and the stability is improved by means of the Levenberg- Marquardt method. Both numerical simulations with noise and experimental results show that the inversion for the bottom reflection loss of small grazing angle has high precision and the certain perturbation in the measured transmission losses does not lead to serious deviation in the inversion result of the bottom reflection loss.

The absorption coefficient measurement of underwater materials in a water tube by using the wide-band pulse method

The existing absorption coefficient measurement in a water tube is limited in low frequency due to the limitation of the pulse tube length.A wide-band pulse tube method based on the post-processing inverse filter is proposed.After obtaining the response of the entire system,the transducer is excited with wide-band,short-duration signals.Then,the received signals are processed with inverse filters and the ＂blurred＂ reflective signals of standard reflector and test sample are recovered respectively,finally,the reflection coefficient and absorption coefficient are calculated.Simulation work is implemented to verify the advantage of postprocessing inverse filter over pre-processing inverse filter in terms of low frequency measurement.To verify the proposed method,two different wide-band pulses in the separate frequency ranges are adopted to measure the same rubber material.The experimental results indicate that the wide-band pulse method and CW（Continuous Wave） pulse method have good agreement in both low frequency range and high frequency range.Wide-band pulse method is of high efficiency,and the low frequency reaches 350 Hz,which is able to extend the low frequency range.