Linearized waveform inversion for vertical transversely isotropic elastic media:Methodology and multi-parameter crosstalk analysis

Seismic migration and inversion are closely related techniques to portray subsurface images and identify hydrocarbon reservoirs.Seismic migration aims at obtaining structural images of subsurface geologic discontinuities.More specifically,seismic migration estimates the reflectivity function(stacked average reflectivity or pre-stack angle-dependent reflectivity)from seismic reflection data.On the other hand,seismic inversion quantitatively estimates the intrinsic rock properties of subsurface formulations.Such seismic inversion methods are applicable to detect hydrocarbon reservoirs that may exhibit lateral variations in the inverted parameters.Although there exist many differences,pre-stack seismic migration is similar with the first iteration of the general linearized seismic inversion.Usually,seismic migration and inversion techniques assume an acoustic or isotropic elastic medium.Unconventional reservoirs such as shale and tight sand formation have notable anisotropic property.We present a linearized waveform inversion(LWI)scheme for weakly anisotropic elastic media with vertical transversely isotropic(VTI)symmetry.It is based on two-way anisotropic elastic wave equation and simultaneously inverts for the localized perturbations(ΔVp_(0)/Vp_(0)/Vs_(0)/Vs_(0)/,Δ∈,Δδ)from the long-wavelength reference model.Our proposed VTI-elastic LWI is an iterative method that requires a forward and an adjoint operator acting on vectors in each iteration.We derive the forward Born approximation operator by perturbation theory and adjoint operator via adjoint-state method.The inversion has improved the quality of the images and reduces the multi-parameter crosstalk comparing with the adjoint-based images.We have observed that the multi-parameter crosstalk problem is more prominent in the inversion images for Thomsen anisotropy parameters.Especially,the Thomsen parameter is the most difficult to resolve.We also analyze the multi-parameter crosstalk using scattering radiation patterns.The linearized waveform inversion for VTI-elastic media presented in this article provides quantitative information of the rock properties that has the potential to help identify hydrocarbon reservoirs.

Full waveform inversion based on hybrid gradient

The low-wavenumber components in the gradient of full waveform inversion(FWI)play a vital role in the stability of the inversion.However,when FWI is implemented in some high frequencies and current models are not far away from the real velocity model,an excessive number of low-wavenumber components in the gradient will also reduce the convergence rate and inversion accuracy.To solve this problem,this paper firstly derives a formula of scattering angle weighted gradient in FWI,then proposes a hybrid gradient.The hybrid gradient combines the conventional gradient of FWI with the scattering angle weighted gradient in each inversion frequency band based on an empirical formula derived herein.Using weighted hybrid mode,we can retain some low-wavenumber components in the initial lowfrequency inversion to ensure the stability of the inversion,and use more high-wavenumber components in the high-frequency inversion to improve the convergence rate.The results of synthetic data experiment demonstrate that compared to the conventional FWI,the FWI based on the proposed hybrid gradient can effectively reduce the low-wavenumber components in the gradient under the premise of ensuring inversion stability.It also greatly enhances the convergence rate and inversion accuracy,especially in the deep part of the model.And the field marine seismic data experiment also illustrates that the FWI based on hybrid gradient(HGFWI)has good stability and adaptability.

First-Arrival Picking Method for Active Source Data with Ocean Bottom Seismometers Based on Spatial Waveform Variation Characteristics

The precision and reliability of first-arrival picking are crucial for determining the accuracy of geological structure inversion using active source ocean bottom seismometer(OBS)refraction data.Traditional methods for first-arrival picking based on sample points are characterized by theoretical errors,especially in low-sampling-frequency OBS data because the travel time of seismic waves is not an integer multiple of the sampling interval.In this paper,a first-arrival picking method that utilizes the spatial waveform variation characteristics of active source OBS data is presented.First,the distribution law of theoretical error is examined;adjacent traces exhibit variation characteristics in their waveforms.Second,a label cross-correlation superposition method for extracting highfrequency signals is presented to enhance the first-arrival picking precision.Results from synthetic and field data verify that the proposed approach is robust,successfully overcomes the limitations of low sampling frequency,and achieves precise outcomes that are comparable with those of high-sampling-frequency data.

Co-Sharing Waveform Design for Millimeter-Wave Radar Communication Systems

Millimeter-wave(mmWave)radar communication has emerged as an important technique for future wireless systems.However,the interference between the radar signal and communication data is the main issue that should be considered for the joint radar communication system.In this paper,a co-sharing waveform(CSW)is proposed to achieve communication and radar sensing simultaneously.To eliminate the co-interference between the communication and sensing signal,signal splitting and processing methods for communication data demodulation and radar signal processing are given respectively.Simulation results show that the bit error rate(BER)of CSW is close to that of the pure communication waveform.Moreover,the proposed CSW can achieve better performance than the existing waveforms in terms of range and velocity estimation.

Underwater Pulse Waveform Recognition Based on Hash Aggregate Discriminant Network

Underwater pulse waveform recognition is an important method for underwater object detection.Most existing works focus on the application of traditional pattern recognition methods,which ignore the time-and space-varying characteristics in sound propagation channels and cannot easily extract valuable waveform features.Sound propagation channels in seawater are time-and space-varying convolutional channels.In the extraction of the waveform features of underwater acoustic signals,the effect of high-accuracy underwater acoustic signal recognition is identified by eliminating the influence of time-and space-varying convolutional channels to the greatest extent possible.We propose a hash aggregate discriminative network(HADN),which combines hash learning and deep learning to minimize the time-and space-varying effects on convolutional channels and adaptively learns effective underwater waveform features to achieve high-accuracy underwater pulse waveform recognition.In the extraction of the hash features of acoustic signals,a discrete constraint between clusters within a hash feature class is introduced.This constraint can ensure that the influence of convolutional channels on hash features is minimized.In addition,we design a new loss function called aggregate discriminative loss(AD-loss).The use of AD-loss and softmax-loss can increase the discriminativeness of the learned hash features.Experimental results show that on pool and ocean datasets,which were collected in pools and oceans,respectively,by using acoustic collectors,the proposed HADN performs better than other comparative models in terms of accuracy and mAP.

Joint Design of ISAC Waveform Under PAPR Constraints

In this paper,we formulate the precoding problem of integrated sensing and communication(ISAC)waveform as a non-convex quadratically constrained quadratic programming(QCQP),in which the weighted sum of communication multi-user interference(MUI)and the gap between dual-use waveform and ideal radar waveform is minimized with peak-toaverage power ratio(PAPR)constraints.We propose an efficient algorithm based on alternating direction method of multipliers(ADMM),which is able to decouple multiple variables and provide a closed-form solution for each subproblem.In addition,to improve the sensing performance in both spatial and temporal domains,we propose a new criteria to design the ideal radar waveform,in which the beam pattern is made similar to the ideal one and the integrated sidelobe level of the ambiguity function in each target direction is minimized in the region of interest.The limited memory Broyden-Fletcher-Goldfarb-Shanno(LBFGS)algorithm is applied to the design of the ideal radar waveform which works as a reference in the design of the dual-function waveform.Numerical results indicate that the designed dual-function waveform is capable of offering good communication quality of service(QoS)and sensing performance.

Electroacupuncture with different waveforms for primary dysmenorrhea:A randomized controlled trial

Objective:To observe and compare the clinical effects of different electroacupuncture waveforms on primary dysmenorrhea.Methods: This was a prospective,randomized,three-group,parallel-controlled trial.Participants with primary dysmenorrhea were randomly divided into dense-sparse wave,continuous wave,and discontinuous wave groups in a 1:1:1 ratio.Two lateral Ciliao(BL 32)points were used.All three groups started treatment 3–5 days before menstruation,once a day for six sessions per course of treatment,one course of treatment per menstrual cycle,and three menstrual cycles.The primary outcome measure was the proportion with an average visual analog scale(VAS)score reduction of≥50%from baseline for dysmenorrhea in the third menstrual cycle during treatment.The secondary outcome measures included changes in dysmenorrhea VAS scores,Cox Menstrual Symptom Scale scores and the proportion of patients taking analgesic drugs.Results: The proportion of cases where the average VAS score for dysmenorrhea decreased by≥50%from baseline in the third menstrual cycle was not statistically significant(P>.05).Precisely 30 min after acupuncture and regarding immediate analgesia on the most severe day of dysmenorrhea,there was a statistically significant difference in the dense-sparse wave group compared with the other two groups during the third menstrual cycle(P<.05).Additionally,there was a statistically significant difference between the dense-sparse wave and discontinuous wave groups 24 h after acupuncture(P<.05).Conclusions: Waveform electroacupuncture can alleviate primary dysmenorrhea and its related symptoms in patients.The three groups showed similar results in terms of short-and long-term analgesic efficacy and a reduction in the number of patients taking analgesic drugs.Regarding achieving immediate analgesia,the dense-sparse wave group was slightly better than the other two groups.

Simultaneous Waveform Inverse Modelling for Litho-Fluid Prediction in an Old Marginal, “Agbbo”Field, Onshore Niger Delta, Nigeria

Simultaneous waveform inversion was used to predict lithofacies and fluid type across the field. Very often, characterizing reservoirs in terms of lithology and fluid type using conventional methods is replete with uncertainties, especially in marginal fields. An approach is employed in this study that integrated rock physics and waveform inverse modelling for lithology and fluid-type characterization to appropriately identify potential hydrocarbon saturated zones and their corresponding lithology. Seismic and well-log data were analyzed using Hampson Russel software. The method adopted includes lithofacies and fluid content analysis using rock physics parameters and seismic simultaneous inverse modelling. Rock physics analysis identified 2 broad reservoirs namely: HDZ1 and HDZ2 reservoirs. Results from the inverse modelling showed that low values of acoustic impedance from 19,743 to 20,487 (ft/s)(g/cc) reflect hydrocarbon-bearing reservoirs while medium to high values shows brine and shale respectively, with brine zone ranging from 20,487 to 22,531 (ft/s)(g/cc) and shale above 22,531 (ft/s)(g/cc). Two lithofacies were identified from inversion analysis of Vp/Vs and Mu-Rho, namely: sand and shale with VpVs 1.95 values respectively. Mu-Rho > 12.29 (GPa)(g/cc) and <12.29 (GPa) (g/cc) represent sand and shale respectively. From 3D volume, it was observed that a high accumulation of hydrocarbon was observed to be saturated at the north to the eastern part of the field forming a meandering channel. Sands were mainly distributed around the northeastern to the southwestern part of the field, that tends to be away from Well 029. This was also validated by the volume of rigidity modulus (Mu-Rho) showing high values indicating sands fall within the northeastern part of the field.

Three dimensional shear wave velocity structure of crust and upper mantle in South China Sea and its adjacent regions by surface waveform inversion

We assembled approximately 328 seismic records. The data set was from 4 digitally recording long-period and broadband stations of CDSN. We carried out the inversion based on the partitioned waveform inversion (PWI). It partitions the large-scale optimization problem into a number of independent small-scale problems. We adopted surface waveform inversion with an equal block (2((2() discretization in order to acquire the images of shear velocity structure at different depths (from surface to 430 km) in the crust and upper-mantle. The resolution of all these anomalies has been established with (check-board( resolution tests. These results show significant difference in velocity, lithosphere and asthenosphere structure between South China Sea and its adjacent regions.

Separation of Ions from Volatile Organic Compounds Using High-Field Asymmetric Waveform Ion Mobility Spectrometry-Mass Spectrometer

A combination of high-field asymmetric waveform ion mobility spectrometry （FAIMS） with mass spectrometer （MS） was analyzed. FAIMS separates ions from the volatile organic compounds in the gas-phase as an ion-filter for MS. The sample ions were created at ambient pressure by ion source, which was equipped with a 10.6 eV UV discharge lamp （A=116.5 nm）. The drift tube of FAIMS is composed of two parallel planar electrodes and the dimension is 10 mm×8 mm×0.5 mm. FAIMS was investigated when driven by the high-filed rectangular asymmetric waveform with the peak-to-peak voltage of 1.36 kV at the frequency of 1 MHz and the duty cycle of 30%. The acetone, the butanone, and their mixture were adopted to characterize the FAIMS-MS. The mass spectra obtained from MS illustrate that there are ion-molecular reactions between the ions and the sample neutral molecular. And the proton transfer behavior in the mixture of the acetone and the butanone is also observed. With the compensation voltage tuned from -30 V to 10 V with a step size of 0.1 V, the ion pre-separation before MS is realized.

Multi-scale seismic full waveform inversion in the frequency-domain with a multi-grid method

Although full waveform inversion in the frequency domain can overcome the local minima problem in the time direction, such problem still exists in the space direction because of the media subsurface complexity. Based on the optimal steep descent methods, we present an algorithm which combines the preconditioned bi-conjugated gradient stable method and the multi-grid method to compute the wave propagation and the gradient space. The multiple scale prosperity of the waveform inversion and the multi-grid method can overcome the inverse problems local minima defect and accelerate convergence. The local inhomogeneous three-hole model simulated results and the Marmousi model certify the algorithm effectiveness.

Effect of projectile head style on high gacceleration waveform of Hopkinson bar calibration system

The freestyle Hopkinson bar is a kind of main high g loading equipment utilized widely in calibration of high g accelerometer and other high shock conditions. The calibration experiment of accelerometer was conducted. With one-dimension stress wave theory, ANSYS/LS-DYNA software and experiment, the effect rules of the projectile＇s front-head style and the accelerometer＇s mounted base＇s length on acceleration waveform were analyzed. The results show that the acceleration duration inspired from Hopkinson bar is almost equal to the rising edge time of perfect half sine stress wave, and it is independent to the mounted base＇s length. Moreover, the projectile＇s fronthead style is a main affecting factor, and the projectiles with less Conical degrees will produce the lower amplitude and longer acceleration duration.

INFLUENCE AND CORRECTION OF SIGNAL SOURCE DISTORTION TO EVALUATION OF EFFECTIVE BITS OF WAVEFORM RECORDERS

In most effective bits evaluation of waveform recorders, the prerequisite is that there is no signal source distortion, or the distortion can be neglected. But when the distortion can be neglected or how it affects the evaluation when it can't be neglected it is not determined yet. In this paper, the influence of signal source distortion to the evaluation of the effective bits of waveform recorders is discussed, then, the correction method of the effective bits error caused by the distortion influence is given. Finally , the error limit of the effective bits is given and how to selecte the calibrator is introduced. In the end , some simulation results of the new method in test are described.

Transmission cable fault detector based on waveform reconstruction

In order to improve the accuracy of cable fault position location at a low cost and make the testing results intuitive, a cable fault detector based on wave form reconstruction is designed. In this detector, the cable fault position is located based on the time-domain pulse reflection （TDR） principle. A pulse waveform is injected in the tested cable, and a high-speed comparator with changeable reference voltages is used to binarize the test pulse waveform to a binary sequence on a certain voltage. Through scanning the reference voltage in a full voltage range, multi-sequences are acquired to reconstruct the pulse waveform transmission in the cable, and then the pulse attenuation feature, electrical open circuit fault, electrical short circuit fault, and the fault position of the cable are diagnosed. Experimental results show that the designed cable fault detector can determine the fault type and its position of the cable being tested, and the testing results are intuitive.

AUTOFOCUSING OF SYNTHETIC HRR RADAR USING STEPPED FREQUENCY WAVEFORMS

Most operating radar systems don′t have sufficient frequency bandwidth to produce high range resolution(HRR) profile of a target. But we can use stepped frequency waveform in a narrow band coherent radar to obtain the HRR profile of a target. For moving targets which are of great importance in practical radar usage, autofocusing,i.e. phase correction, is a necessary and critical step of the synthetic HRR processing. The purpose of autofocusing is to remove the radial motion effect of the target from radar echoes, and only reserve the stepped frequency effect which is the basis of synthetic HRR capability. We investigate two autofocusing approaches for synthetic HRR radars using stepped frequency waveform in this paper. The first is motion fitting method. This method depends on a certain parametric model, and is computationally expensive. Then we propose the iterative dominant scatterer method. It is robust, non parametric and simple in computation in comparison with the motion fitting method. Experimental results based on data acquired by using a metallised scale model B 52 in a microwave anechoic chamber reveal the validity and effectiveness of the method.

Improved hybrid iterative optimization method for seismic full waveform inversion

In full waveform inversion （FWI）, Hessian information of the misfit function is of vital importance for accelerating the convergence of the inversion; however, it usually is not feasible to directly calculate the Hessian matrix and its inverse. Although the limited memory Broyden-Fletcher-Goldfarb-Shanno （L-BFGS） or Hessian-free inexact Newton （HFN） methods are able to use approximate Hessian information, the information they collect is limited. The two methods can be interlaced because they are able to provide Hessian information for each other; however, the performance of the hybrid iterative method is dependent on the effective switch between the two methods. We have designed a new scheme to realize the dynamic switch between the two methods based on the decrease ratio （DR） of the misfit function （objective function）, and we propose a modified hybrid iterative optimization method. In the new scheme, we compare the DR of the two methods for a given computational cost, and choose the method with a faster DR. Using these steps, the modified method always implements the most efficient method. The results of Marmousi and overthrust model testings indicate that the convergence with our modified method is significantly faster than that in the L-BFGS method with no loss of inversion quality. Moreover, our modified outperforms the enriched method by a little speedup of the convergence. It also exhibits better efficiency than the HFN method.

Detection optimization for resonance region radar with dense multi-carrier waveform

Unlike the existing resonance region radar systems （RRRS ） that transmit the orthogonal frequency division multiplexing （OFDM）multi-carrier waveform,the dense multi-carrier （DMC）radar waveform which has a narrower frequency interval than the traditional OFDM waveform is proposed.Therefore,in the same frequency bandwidth,the DMC waveform contains more sub-carriers and provides more frequency diversity.Additionally,to further improve detection performance,a novel optimal weight accumulation target detection （OWATD）method is proposed,where the echo electromagnetic waves at different frequencies are accumulated with the optimal weight coefficients.Then,with the signal-to-noise ratio （SNR）of echo waveform approaching infinity,the asymptotic detection performance is analyzed, and the condition that the OWATD method with the DMC outperforms the matched filter with the OFDM is presented.Simulation results show that the DMC outperforms the OFDM in the target detection performance,and the OWATD method can further improve the detection performance of the traditional methods with both the OFDM and DMC radar waveform.

Frequency-domain auto-adapting full waveform inversion with blended source and frequency-group encoding

As a high quality seismic imaging method, full waveform inversion （FWI） can accurately reconstruct the physical parameter model for the subsurface medium. However, application of the FWI in seismic data processing is computationally expensive, especially for the three-dimension complex medium inversion. Introducing blended source technology into the frequency-domain FWI can greatly reduce the computational burden and improve the efficiency of the inversion. However, this method has two issues： first, crosstalk noise is caused by interference between the sources involved in the encoding, resulting in an inversion result with some artifacts; second, it is more sensitive to ambient noise compared to conventional FWI, therefore noisy data results in a poor inversion. This paper introduces a frequency-group encoding method to suppress crosstalk noise, and presents a frequency- domain auto-adapting FWI based on source-encoding technology. The conventional FWI method and source-encoding based FWI method are combined using an auto-adapting mechanism. This improvement can both guarantee the quality of the inversion result and maximize the inversion efficiency.

Determination of focal depth by two waveform-based methods:A case study for the 2008 Panzhihua earthquake

With the 2008 Ms6.1 Panzhihua earthquake as a case study, we demonstrate that the focal depth of the main shock can be well constrained with two approaches： （1） using the depth phase sPL and （2） using full waveform inversion of local and teleseismic data. We also show that focal depths can be well constrained using the depth phase sPL with single broadband seismic station. Our study indicates that the main shock is located at a depth of ii kin, much shallower than those from other studies, confirming that the earthquake occurs in upper crust. Aftershocks are located in the depth range of 11 16 kin, which is consistent with a ruptured near vertical fault whose width is about 10 km, as expected for an Ms6.1 earthquake.

Upper-mantle velocity structures beneath the Tibetan Plateau and surrounding areas inferred from triplicated P waveforms

P-wave waveforms in the distance range between 12°and 30°were analyzed to investigate upper-mantle P velocity structures beneath the Tibetan Plateau and surrounding areas.The waveform data from 504 earthquakes with magnitudes larger than 5.0 between 1990 and 2005 that occurred within 30°from the center of the Plateau were modelled.We divided the study area into 6 regions and modeled upper-mantle-distance P waveforms with turning points beneath each region separately.The results show that the uppermantle P-wave velocity structures beneath India,the Himalayas,and the Lhasa Terrane are similar and contain a high-velocity lid about 250 km thick.The upper-mantle velocities down to 200 km beneath the Qiangtang Terrane,the Tarim Basin,and especially the Songpan-GarzêTerrane are lower than those in the south.The 410-km discontinuity beneath these two terranes is elevated by about 20 km.Highvelocity anomalies are found in the transition zone below 500 km under the Lhasa and Qiangtang Terranes.The results suggest that the Tibetan Plateau was generated by thrusting of the Indian mantle lithosphere under the southern part of Tibet.Portions of the thickened Eurasian mantle lithosphere were delaminated;they are now sitting in the transition zone beneath southern Tibet and atop of the 410-km discontinuity underneath northern Tibet.