Time-Lapse Full-Waveform Inversion Using Cross-Correlation-Based Dynamic Time Warping

Offshore carbon capture, utilization, and storage(OCCUS) is regarded as a crucial technology for mitigating greenhouse gas emissions.Quantitative monitoring maps of sealed carbon dioxide are necessary in a comprehensive OCCUS project. A potential high-resolution method for the aforementioned purpose lies in the full-waveform inversion(FWI) of time-lapse seismic data. However, practical applications of FWI are severely restricted by the well-known cycle-skipping problem. A new time-lapse FWI method using cross-correlation-based dynamic time warping(CDTW) is proposed to detect changes in the subsurface property due to carbon dioxide(CO_(2)) injection and address the aforementioned issue. The proposed method, namely CDTW, which combines the advantages of cross-correlation and dynamic time warping, is employed in the automatic estimation of the discrepancy between the seismic signals simulated using the baseline/initial model and those acquired. The proposed FWI method can then back-project the estimated discrepancy to the subsurface space domain, thereby facilitating retrieval of the induced subsurface property change by taking the difference between the inverted baseline and monitor models. Numerical results on pairs of signals prove that CDTW can obtain reliable shifts under amplitude modulation and noise contamination conditions. The performance of CDTW substantially outperforms that of the conventional dynamic time warping method. The proposed time-lapse fullwaveform inversion(FWI) method is applied to the Frio-2 CO_(2) storage model. The baseline and monitor models are inverted from the corresponding time-lapse seismic data. The changes in velocity due to CO_(2) injection are reconstructed by the difference between the baseline and the monitor models.

Inversion of Seabed Geotechnical Properties in the Arctic Chukchi Deep Sea Basin Based on Time Domain Adaptive Search Matching Algorithm

The chirp sub-bottom profiler,for its high resolution,easy accessibility and cost-effectiveness,has been widely used in acoustic detection.In this paper,the acoustic impedance and grain size compositions were obtained based on the chirp sub-bottom profiler data collected in the Chukchi Plateau area during the 11th Arctic Expedition of China.The time-domain adaptive search matching algorithm was used and validated on our established theoretical model.The misfit between the inversion result and the theoretical model is less than 0.067%.The grain size was calculated according to the empirical relationship between the acoustic impedance and the grain size of the sediment.The average acoustic impedance of sub-seafloor strata is 2.5026×10^(6) kg(s m^(2))^(-1)and the average grain size(θvalue)of the seafloor surface sediment is 7.1498,indicating the predominant occurrence of very fine silt sediment in the study area.Comparison of the inversion results and the laboratory measurements of nearby borehole samples shows that they are in general agreement.

A New Proposal of Inverse-Time Characteristics for Overcurrent Relay Based on Microcomputer

The conventional time function of electromechanical relays is hard to coordinate with other relays. In order to promote the application of inverse-time overcurrent relays, a new time function for microprocessor-type relay is proposed. The setting of the trip time for this relay is performed by determining the shortest trip time and the longest trip time, respectively. The results of analysis show that with the new time function, the inverse-time overcurrent relay is easy to coordinate with other relays and has a comparatively shorter trip time, and that the fault happens in the protective zone.

Inverse problem of quadratic time-dependent Hamiltonians

Using an algebraic approach, it is possible to obtain the temporal evolution wave function for a Gaussian wavepacket obeying the quadratic time-dependent Hamiltonian(QTDH). However, in general, most of the practical cases are not exactly solvable, for we need general solutions of the Riccatti equations which are not generally known. We therefore bypass directly solving for the temporal evolution wave function, and study its inverse problem. We start with a particular evolution of the wave-packet, and get the required Hamiltonian by using the inverse method. The inverse approach opens up a new way to find new exact solutions to the QTDH. Some typical examples are studied in detail. For a specific timedependent periodic harmonic oscillator, the Berry phase is obtained exactly.

Clustering-Inverse: A Generalized Model for Pattern-Based Time Series Segmentation

Patterned-based time series segmentation (PTSS) is an important task for many time series data mining applications. In this paper, according to the characteristics of PTSS, a generalized model is proposed for PTSS. First, a new inter-pretation for PTSS is given by comparing this problem with the prototype-based clustering (PC). Then, a novel model, called clustering-inverse model (CI-model), is presented. Finally, two algorithms are presented to implement this model. Our experimental results on artificial and real-world time series demonstrate that the proposed algorithms are quite effective.

Shifted first arrival point travel time NMO inversion

Serious stretch appears in shallow long offsset signals after NMO correction. In this article we study the generation mechanism of NMO stretch, demonstrate that the conventional travel time equation cannot accurately describe the travel time of the samples within the same reflection wavelet. As a result, conventional NMO inversion based on the travel time of the wavelet＇s central point occurs with errors. In this article, a travel time equation for the samples within the same wavelet is reconstructed through our theoretical derivation （the shifted first arrival point travel time equation）, a new NMO inversion method based on the wavelet＇s first arrival point is proposed. While dealing with synthetic data, the semblance coefficient algorithm equation is modified so that wavelet first arrival points can be extracted. After that, NMO inversion based on the new velocity analysis is adopted on shot offset records. The precision of the results is significantly improved compared with the traditional method. Finally, the block move NMO correction based on the first arrival points travel times is adopted on long offset records and non-stretched results are achieved, which verify the proposed new equation.

A time domain induced polarization relaxation time spectrum inversion method based on a damping factor and residual correction

Relaxation time spectra （RTS） derived from time domain induced polarization data （TDIP） are helpful to assess oil reservoir pore structures. However, due to the sensitivity to the signal-to-noise ratio （SNR）, the inversion accuracy of the traditional singular value decomposition （SVD） inversion method reduces with a decrease of SNR. In order to enhance the inversion accuracy and improve robustness of the inversion method to the SNR, an improved inversion method, based on damping factor and spectrum component residual correction, is proposed in this study. The numerical inversion results show that the oscillation of the RTS derived from the SVD method increased with a decrease of SNR, which makes it impossible to get accurate inversion components. However, the SNR has little influence on inversion components of the improved method, and the RTS has high inversion accuracy and robustness. Moreover, RTS derived from core sample data is basically in accord with the pore-size distribution curve, and the RTS derived from the actual induced polarization logging data is smooth and continuous, which indicates that the improved method is practicable.

Joint AVO inversion in the time and frequency domain with Bayesian interference

Amplitude variations with offset or incident angle （AVO/AVA） inversion are typically combined with statistical methods, such as Bayesian inference or deterministic inversion. We propose a joint elastic inversion method in the time and frequency domain based on Bayesian inversion theory to improve the resolution of the estimated P- and S-wave velocities and density. We initially construct the objective function using Bayesian inference by combining seismic data in the time and frequency domain. We use Cauchy and Gaussian probability distribution density functions to obtain the prior information for the model parameters and the likelihood function, respectively. We estimate the elastic parameters by solving the initial objective function with added model constraints to improve the inversion robustness. The results of the synthetic data suggest that the frequency spectra of the estimated parameters are wider than those obtained with conventional elastic inversion in the time domain. In addition, the proposed inversion approach offers stronger antinoising compared to the inversion approach in the frequency domain. Furthermore, results from synthetic examples with added Gaussian noise demonstrate the robustness of the proposed approach. From the real data, we infer that more model parameter details can be reproduced with the proposed joint elastic inversion.

3D parallel inversion of time-domain airborne EM data

To improve the inversion accuracy of time-domain airborne electromagnetic data, we propose a parallel 3D inversion algorithm for airborne EM data based on the direct Gauss-Newton optimization. Forward modeling is performed in the frequency domain based on the scattered secondary electrical field. Then, the inverse Fourier transform and convolution of the transmitting waveform are used to calculate the EM responses and the sensitivity matrix in the time domain for arbitrary transmitting waves. To optimize the computational time and memory requirements, we use the EM ＂footprint＂ concept to reduce the model size and obtain the sparse sensitivity matrix. To improve the 3D inversion, we use the OpenMP library and parallel computing. We test the proposed 3D parallel inversion code using two synthetic datasets and a field dataset. The time-domain airborne EM inversion results suggest that the proposed algorithm is effective, efficient, and practical.

Quasi-real time inversion method of three-dimensional epicenter coordinate, trigger time, and magnitude based on CORS

This study explores the quasi-real time inversion principle and precision estimation of three-dimensional coordinates of the epicenter, trigger time and magnitude of earthquakes with the aim to improve traditional methods, which are flawed due to missing information or distortion in the seismograph records. The epicenter, trigger time and magnitude from the Lushan earthquake are inverted and analyzed based on high-frequency GNSS data. The inversion results achieved a high precision, which are consistent with the data published by the China Earthquake Administration. Moreover, it has been proven that the inversion method has good theoretical value and excellent application prospects.

Frequency Time Analysis (FTAN) and Moment Tensor Inversion Solutions from Short Period Surface Waves in Cameroon (Central Africa)

Short period surface waves generated by a local earthquake recorded by broadband seismometers at distances of about 186 to 778 km from the earthquake’s epicenter located in Cameroon (Central Africa) were processed for group velocity maps and dispersion waveforms using the frequency time analysis (FTAN) method. The resulting group velocity fundamental modes of the extracted Rayleigh and Love waves were used for a joint amplitude spectral and P polarity inversion using moment tensor inversion. The corresponding group velocity dispersion curves, the residual as a function of depth, the amplitude spectra and the moment tensor solutions of the regions from the epicenter to the different stations up to a depth of about 10 km were obtained.

Geoacoustic Inversion Using Time Reversal of Ocean Noise

We present a passive geoacoustic inversion method using two hydrophones, which combines noise interferometry and time reversal mirror （TRM） techniques. Numerical simulations are firstly performed, in which strong fo- cusing occurs in the vicinity of one hydrophone when Green＇s function （GF） is back-propagated from the other hydrophone, with the position and strength of the focus being sensitive to sound speed and density in the bottom. We next extract the GF from the noise cross-correlation function measured by two hydrophones with 8025-m distance in the Shallow Water ＇06 experiment. After realizing the TRM process, sound speed and density in the bottom are inverted by optimizing focusing of the back-propagated GF. The passive inversion method is inherently environmentally friendly and low-cost.

Diagonally loaded SMI algorithm based on inverse matrix recursion

The derivation of a diagonally loaded sample-matrix inversion （LSMI） algorithm on the busis of inverse matrix recursion （i.e.LSMI-IMR algorithm） is conducted by reconstructing the recursive formulation of covariance matrix. For the new algorithm, diagonal loading is by setting initial inverse matrix without any addition of computation. In addition, a corresponding improved recursive algorithm is presented, which is low computational complexity. This eliminates the complex multiplications of the scalar coefficient and updating matrix, resulting in significant computational savings. Simulations show that the LSMI-IMR algorithm is valid.

Relevance vector machine technique for the inverse scattering problem

A novel method based on the relevance vector machine（RVM） for the inverse scattering problem is presented in this paper.The nonlinearity and the ill-posedness inherent in this problem are simultaneously considered.The nonlinearity is embodied in the relation between the scattered field and the target property,which can be obtained through the RVM training process.Besides,rather than utilizing regularization,the ill-posed nature of the inversion is naturally accounted for because the RVM can produce a probabilistic output.Simulation results reveal that the proposed RVM-based approach can provide comparative performances in terms of accuracy,convergence,robustness,generalization,and improved performance in terms of sparse property in comparison with the support vector machine（SVM） based approach.

2-D crustal Poisson′s ratio from seismic travel time inversion in Changbaishan Tianchi volcanic region

Based on the inversion method of 2D velocity structure and interface, the crustal velocity structures of P-wave and S-wave along the profile L1 are determined simultaneously with deep seismic sounding data in Changbaishan Tianchi volcanic region, and then its Poisson's ratio is obtained. Calculated results show that this technique overcomes some defects of traditional forward calculation method, and it is also very effective to determine Poisson's ratio distribution of deep seismic sounding profile, especially useful for study on volcanic magma and crustal fault zone. Study result indicates that there is an abnormally high Poisson's ratio body that is about 30 km wide and 12 km high in the low velocity region under Tianchi crater. Its value of Poisson's ratio is 8% higher than that of surrounding medium and it should be the magma chamber formed from melted rock with high temperature. There is a high Poisson's ratio zone ranging from magma chamber to the top of crust, which may be the uprise passage of hot substance. The lower part with high Poisson's ratio, which stretches downward to Moho, is possibly the extrusion way of hot substance from the uppermost mantle. The conclusions above are consistent with the study results of both tomographic determination of 3D crustal structure and magnetotelluric survey in this region.

Some Equations for the External Space and Their Implications on the Inverse Cone Concept and the Geometry of the Physical Universe

In previous works, they were proposed a photonic model of the Big Bang [1] and several parameters derived from the Hubble-Lemaitre equation [2]. Since these parameters result higher than the classical ones and, otherwise, the General Theory of Relativity does not apply far away the Physical Universe, in this paper, it will be revised the adequacy of such parameters in the external Space and their influence on the relativistic concept of the cone of time. As well, it will be intended to define the Physical Universe geometry accordingly to a thermo-dynamical analysis of the Big Bang.

A NEW INVERSION METHOD OF TIME-LAPSE SEISMIC

Time-Lapse Seismic improves oil recovery ratio by dynamic reservoir monitoring. Because of the large number of seismic explorations in the process of time-lapse seismic inversion, traditional methods need plenty of inversion calculations which cost high computational works. The method is therefore inefficient. In this paper, in order to reduce the repeating computations in traditional, a new time-lapse seismic inversion method is put forward. Firstly a homotopy-regularization method is proposed for the first time inversion. Secondly, with the first time inversion results as the initial value of following model, a model of the second time inversion is rebuilt by analyzing the characters of time-lapse seismic and localized inversion method is designed by using the model. Finally, through simulation, the comparison between traditional method and the new scheme is given. Our simulation results show that the new scheme could save the algorithm computations greatly.

Range-instantaneous Doppler imaging of inverse synthetic aperture sonar

Because the existing range-Doppler algorithm in inverse synthetic aperture sonar (ISAS) is based on target model of uniform motion, it may be invalidated for maneuvering targets due to the time-varying changes of both individual scatter′s Doppler and imaging projection plane. To resolve the problem, a new range-instantaneous Doppler imaging method is proposed for imaging maneuvering targets based on time-frequency analysis. The proposed approach is verified using real underwater acoustic data.

More Peripheral Visualization of Hepatic Arteries by Using Respiratory-Triggered 3D True Steady-State Free-Precession Projection Magnetic Resonance Angiographic Sequences with Time-Spatial Labeling Inversion Pulse

Purpose: To evaluate respiratory-triggered three-dimensional (3D) true steady-state free-precession (SSFP) projection magnetic resonance angiographic sequences with time-spatial labeling inversion pulse (Time-SLIP) for visualizing the hepatic arteries and to optimize the image acquisition protocol. Materials and Methods: A 1.5-T clinical magnetic resonance imager was used to perform abdominal magnetic resonance angiography (MRA) in 25 consecutive patients before transcatheter arterial chemoembolization or surgery. We compared two selective space-labeling inversion pulse (tag pulse) patterns (Patterns I and II, oblique and parallel tag pulses, respectively). Two experienced radiologists evaluated the number of hepatic arterial branches visible on the acquired MRA images, and the results were referenced with those on images from intra-arterial digital subtraction angiography. Results: Images were acquired from all patients. The two radiologists clearly visualized branches of the left and right hepatic arteries. More peripheral hepatic arterial branches were identified in MRA images captured by using tag pulse Pattern I than in those acquired by using Pattern II (P P > 0.05). Conclusion: Non-contrast-enhanced Time-SLIP hepatic MRA with true SSFP allowed selective visualization of peripheral hepatic vessels.

基岩中入射波三维时域反演及自由场构建

基岩入射波一维和二维时域反演不能反映入射方位角和斜入射角的空间任意性,导致构建的自由场可能与实际偏离较远,因此有必要开展基岩入射波三维时域反演研究。基于波函数组合法,利用基岩中P波、SV波和SH波的入射波时程表达地表控制点自由场,建立基于设计地震动的基岩入射波三维时域反演方法,并进一步构建了多波组合斜入射下的空间非一致自由场,分析了基于一维和二维反演方法构建的自由场偏差,揭示了不同反演和自由场构建方法对于沥青混凝土心墙土石坝地震响应的影响机制。结果表明:基于地震动三维反演方法构建的自由场具有空间非一致性,能够更全面地反映场地的地震动场情况。当入射方位角在0°~60°、斜入射角在40°~90°范围内,与三维反演方法相比,一维反演方法构建的特征点x向自由场偏差较大,位移峰值最大偏差为55.8%;当入射方向与某个水平坐标轴平行时,二维反演方法构建的特征点另一个水平向自由场偏差较大,位移峰值最大偏差为100%。介质泊松比增大,一维反演下x向自由场偏差减小,二维反演下x向自由场偏差增大,特征点位移峰值最大偏差分别为46%和36%。与基于地震动三维反演方法获得的心墙拉应力最大值相比,二维和一维反演下拉应力最大值分别减小了83.3%和20.0%。