SOLVING THE VELOCITY PARAMETER OF THE 2-D WAVE INVERSE PROBLEMS WITH THE INTEGRATION-CHARACTERISTIC METHOD

For the 2-D wave inverse problems introduced from geophysical exploration, in this paper, the author presents integration-characteristic method to solve the velocity parameter, and then applies it to common shotpoint model data, in noise-free case. The accuracy is quite good.

Fracture network types revealed by well test curves for shale reservoirs in the Sichuan Basin,China

Pressure buildup testing can be used to analyze fracture network characteristics and conduct quantitative interpretation of relevant parameters for shale gas wells,thus providing bases for assessing the well productivity and formulating proper development strategies.This study establishes a new well test interpretation model for fractured horizontal wells based on seepage mechanisms of shale reservoirs and proposes a method for identifying fracturing patterns based on the characteristic slopes of pressure buildup curves and curve combination patterns.The pressure buildup curve patterns are identified to represent three types of shale reservoirs in the Sichuan Basin,namely the moderately deep shale reservoirs with high pressure,deep shale reservoirs with ultra-high pressure,and moderately deep shale reservoirs with normal pressure.Based on this,the relationship between the typical pressure buildup curve patterns and the fracture network types are put forward.Fracturing effects of three types of shale gas reservoir are compared and analyzed.The results show that typical flow patterns of shale reservoirs include bilinear flow in primary and secondary fractures,linear flow in secondary fractures,bilinear flow in secondary fractures and matrix,and linear flow in matrix.The fracture network characteristics can be determined using the characteristic slopes of pressure buildup curves and curve combinations.The linear flow in early secondary fractures is increasingly distinct with an increase in primary fracture conductivity.Moreover,the bilinear flow in secondary fractures and matrix and the subsequent linear flow in the matrix occur as the propping and density of secondary fractures increase.The increase in the burial depth,in-situ stress,and stress difference corresponds to a decrease in the propping of primary fractures that expand along different directions in the shale gas wells in the Sichuan Basin.Four pressure buildup curve patterns exist in the Sichuan Basin and its periphery.The pattern of pressure buildup curves of shale reservoirs in the Yongchuan area can be described as 1/2/→1/4,indicating limited stimulated reservoir volume,poorly propped secondary fractures,and the forming of primary fractures that extend only to certain directions.The pressure buildup curves of shale reservoirs in the main block of the Fuling area show a pattern of 1/4/→1/2 or 1/2,indicating greater stimulated reservoir volume,well propped secondary fractures,and the forming of complex fracture networks.The pattern of pressure buildup curves of shale reservoirs in the Pingqiao area is 1/2/→1/4→/1/2,indicating a fracturing effect somewhere between that of the Fuling and Yongchuan areas.For reservoirs with normal pressure,it is difficult to determine fracture network characteristics from pressure buildup curves due to insufficient formation energy and limited liquid drainage.

The constructing of pore structure factor in carbonate rocks and the inversion of reservoir parameters

With a more complex pore structure system compared with clastic rocks, carbonate rocks have not yet been well described by existing conventional rock physical models concerning the pore structure vagary as well as the influence on elastic rock properties. We start with a discussion and an analysis about carbonate rock pore structure utilizing rock slices. Then, given appropriate assumptions, we introduce a new approach to modeling carbonate rocks and construct a pore structure algorithm to identify pore structure mutation with a basis on the Gassmann equation and the Eshelby-Walsh ellipsoid inclusion crack theory. Finally, we compute a single well＇s porosity using this new approach with full wave log data and make a comparison with the predicted result of traditional method and simultaneously invert for reservoir parameters. The study results reveal that the rock pore structure can significantly influence the rocks＇ elastic properties and the predicted porosity error of the new modeling approach is merely 0.74%. Therefore, the approach we introduce can effectively decrease the predicted error of reservoir parameters.

Parameter optimization of gravity density inversion based on correlation searching and the golden section algorithm

For density inversion of gravity anomaly data, once the inversion method is determined, the main factors affecting the inversion result are the inversion parameters and subdivision scheme. A set of reasonable inversion parameters and subdivision scheme can, not only improve the inversion process efficiency, but also ensure inversion result accuracy. The gravity inversion method based on correlation searching and the golden section algorithm is an effective potential field inversion method. It can be used to invert 2D and 3D physical properties with potential data observed on flat or rough surfaces. In this paper, we introduce in detail the density inversion principles based on correlation searching and the golden section algorithm. Considering that the gold section algorithm is not globally optimized. we present a heuristic method to ensure the inversion result is globally optimized. With a series of model tests, we systematically compare and analyze the inversion result efficiency and accuracy with different parameters. Based on the model test results, we conclude the selection principles for each inversion parameter with which the inversion accuracy can be obviously improved.

Deformation early-warning index for heightened gravity dam during impoundment period

The mechanical parameters of materials in a dam body and dam foundation tend to change when dams are reinforced in aging processes.It is important to use an early-warning index to reflect the safety status of dams,particularly of heightened projects in the impoundment period.Herein,a new method for monitoring the safety status of heightened dams is proposed based on the deformation monitoring data of a dam structure,a statistical model,and finite-element numerical simulation.First,a fast optimization inversion method for estimation of dam mechanical parameters was developed,which used the water pressure component extracted from a statistical model,an improved inversion objective function,and a genetic optimization iterative algorithm.Then,a finite element model of a heightened concrete gravity dam was established,and the deformation behavior of the dam with rising water levels in the impoundment period was simulated.Subsequently,mechanical parameters of aged dam parts were calculated using the fast optimization inversion method with simulated deformation and the water pressure deformation component obtained by the statistical model under the same conditions of water pressure change.Finally,a new earlywarning index of dam deformation was constructed by means of the forward-simulated deformation and other components of the statistical model.The early-warning index is useful for forecasting dam deformation under different water levels,especially high water levels.

Deformation warning index for reinforced concrete dam based on structural health monitoring data and numerical simulation

The material mechanical parameters of the dam body and foundation will change when a dam is reinforced during the aging process.This causes significant changes in the structural state of the project and makes it difficult to ensure its structural safety.In this study,a new deformation warning index for reinforced concrete dams was developed according to the prototype monitoring data,statistical models,three-dimensional finite element model(FEM)numerical simulation,and the critical conditions of the dam structure.A statistical model was established to separate the water pressure component.Then,a three-dimensional FEM of the reinforced concrete dam was constructed to simulate the water pressure component.Furthermore,the deformation components that affected the mechanical parameters of the dam under the same amount of reservoir water level change were separated and quantified accurately.In addition,the method for inversion of comprehensive mechanical parameters after dam reinforcement was used.The influence mechanisms of the deformation behavior of concrete dams under the reservoir water level and temperature changes were investigated.A new deformation warning index was developed by combining the forward-simulated critical water pressure component and temperature component in the period of extreme temperature decrease with the aging component separated by the statistical model.The new deformation warning index considers the structural state of the dam before and after reinforcement and links the structural strength criterion and the deformation evolution mechanisms.It provides a theoretical foundation and decision support for long-term service and operation management of reinforced dams.

Influence of correlation scale errors on aquifer hydraulic conductivity inversion precision

In order to investigate the influence of correlation scale error on the inversion precision of the hydraulic conductivity of the aquifer,the successive linear estimator(SLE)was used to invert the hydraulic conductivity field of a heterogeneous aquifer based on synthetic experiments.By increasing the numbers of observation wells and pumping tests,we analyzed the difference between the estimated and true values of hydraulic conductivity with different correlation scale errors.The relationships between the observation well number and the error in inversion results,and between the pumping test number and the error in inversion results were investigated.The results show that,if the amount of observed head data is insufficient,there will be errors in inversion results with changing correlation scale.Due to the existence of correlation scale error,the improvement of inversion precision gradually slows down with the increase of the amount of observed head data,which indicates that too much observed head data causes data redundancy.Therefore,for the synthetic experiments described in this paper,the observation well number should be less than 41,the pumping test number should be less than 17,and a more suitable method should be selected according to the precision requirements of specific situations in practical engineering.

Genetic Nelder-Mead neural network algorithm for fault parameter inversion using GPS data

The traditional genetic algorithm(GA)has unstable inversion results and is easy to fall into the local optimum when inverting fault parameters.Therefore,this article considers the combination of GA with other non-linear algorithms in order to improve the inversion precision of GA.This paper proposes a genetic Nelder-Mead neural network algorithm(GNMNNA).This algorithm uses a neural network algorithm(NNA)to optimize the global search ability of GA.At the same time,the simplex algorithm is used to optimize the local search capability of the GA.Through numerical examples,the stability of the inversion algorithm under different strategies is explored.The experimental results show that the proposed GNMNNA has stronger inversion stability and higher precision compared with the existing algorithms.The effectiveness of GNMNNA is verified by the BodrumeKos earthquake and Monte Cristo Range earthquake.The experimental results show that GNMNNA is superior to GA and NNA in both inversion precision and computational stability.Therefore,GNMNNA has greater application potential in complex earthquake environment.

Inversion of Oceanic Parameters Represented by CTD Utilizing Seismic Multi-Attributes Based on Convolutional Neural Network

In Recent years,seismic data have been widely used in seismic oceanography for the inversion of oceanic parameters represented by conductivity temperature depth(CTD).Using this technique,researchers can identify the water structure with high horizontal resolution,which compensates for the deficiencies of CTD data.However,conventional inversion methods are modeldriven,such as constrained sparse spike inversion(CSSI)and full waveform inversion(FWI),and typically require prior deterministic mapping operators.In this paper,we propose a novel inversion method based on a convolutional neural network(CNN),which is purely data-driven.To solve the problem of multiple solutions,we use stepwise regression to select the optimal attributes and their combination and take two-dimensional images of the selected attributes as input data.To prevent vanishing gradients,we use the rectified linear unit(ReLU)function as the activation function of the hidden layer.Moreover,the Adam and mini-batch algorithms are combined to improve stability and efficiency.The inversion results of field data indicate that the proposed method is a robust tool for accurately predicting oceanic parameters.

A SCR method for uncertainty estimation in geodesy non-linear error propagation: Comparisons and applications

We review three derivative-free methods developed for uncertainty estimation of non-linear error propagation, namely, MC(Monte Carlo), SUT(scaled unscented transformation), and SI(sterling interpolation). In order to avoid preset parameters like as these three methods need, we introduce a new method to uncertainty estimation for the first time, namely, SCR(spherical cubature rule), which is no need for setting parameters. By theoretical derivation, we prove that the precision of uncertainty obtained by SCR can reach second-order. We conduct four synthetic experiments, for the first two experiments, the results obtained by SCR are consistent with the other three methods with optimal setting parameters, but SCR is easier to operate than other three methods, which verifies the superiority of SCR in calculating the uncertainty. For the third experiment, real-time calculation is required, so the MC is hardly feasible. For the forth experiment, the SCR is applied to the inversion of seismic fault parameter which is a common problem in geophysics, and we study the sensitivity of surface displacements to fault parameters with errors. Our results show that the uncertainty of the surface displacements is the magnitude of ±10 mm when the fault length contains a variance of 0.01 km^(2).

A smart productivity evaluation method for shale gas wells based on 3D fractal fracture network model

The generation method of three-dimensional fractal discrete fracture network(FDFN)based on multiplicative cascade process was developed.The complex multi-scale fracture system in shale after fracturing was characterized by coupling the artificial fracture model and the natural fracture model.Based on an assisted history matching(AHM)using multiple-proxy-based Markov chain Monte Carlo algorithm(MCMC),an embedded discrete fracture modeling(EDFM)incorporated with reservoir simulator was used to predict productivity of shale gas well.When using the natural fracture generation method,the distribution of natural fracture network can be controlled by fractal parameters,and the natural fracture network generated coupling with artificial fractures can characterize the complex system of different-scale fractures in shale after fracturing.The EDFM,with fewer grids and less computation time consumption,can characterize the attributes of natural fractures and artificial fractures flexibly,and simulate the details of mass transfer between matrix cells and fractures while reducing computation significantly.The combination of AMH and EDFM can lower the uncertainty of reservoir and fracture parameters,and realize effective inversion of key reservoir and fracture parameters and the productivity forecast of shale gas wells.Application demonstrates the results from the proposed productivity prediction model integrating FDFN,EDFM and AHM have high credibility.

Gas emission source term estimation with 1-step nonlinear partial swarm optimization-Tikhonov regularization hybrid method

Source term identification is very important for the contaminant gas emission event. Thus, it is necessary to study the source parameter estimation method with high computation efficiency, high estimation accuracy and reasonable confidence interval. Tikhonov regularization method is a potential good tool to identify the source parameters. However, it is invalid for nonlinear inverse problem like gas emission process. 2-step nonlinear and linear PSO （partial swarm optimization）-Tikhonov regularization method proposed previously have estimated the emission source parameters successfully. But there are still some problems in computation efficiency and confidence interval. Hence, a new 1-step nonlinear method combined Tikhonov regularizafion and PSO algorithm with nonlinear forward dispersion model was proposed. First, the method was tested with simulation and experiment cases. The test results showed that 1-step nonlinear hybrid method is able to estimate multiple source parameters with reasonable confidence interval. Then, the estimation performances of different methods were compared with different cases. The estimation values with 1-step nonlinear method were close to that with 2-step nonlinear and linear PSO-Tikhonov regularization method, 1-step nonlinear method even performs better than other two methods in some cases, especially for source strength and downwind distance estimation. Compared with 2-step nonlinear method, 1-step method has higher computation efficiency. On the other hand, the confidence intervals with the method proposed in this paper seem more reasonable than that with other two methods. Finally, single PSO algorithm was compared with 1-step nonlinear PSO-Tikhonov hybrid regularization method. The results showed that the skill scores of 1-step nonlinear hybrid method to estimate source parameters were close to that of single PSO method and even better in some cases. One more important property of 1-step nonlinear PSO-Tikhonov regularization method is its reasonable confidence interval, which is not obtained by single PSO algorithm. Therefore, 1-step nonlinear hybrid regularization method proposed in this paper is a potential good method to estimate contaminant gas emission source term.

A parameter inversion for sea bridge based on high-resolution polarimetric synthetic aperture radar

Each reflection return of a bridge over water is displayed as wide stripe in a high-resolution synthetic aperture radar （SAR） image, which lead to difficulties in a parameter inversion. Therefore, a method of bridge parameter inversion is proposed for high-resolution full polarimetric SAR （PolSAR）. First, the single, double and triple- bounce returns from each component of the bridge are distinguished by the polarization scattering features. Then the reasons which lead to the backscatter echoes oft_he bridge over water being displayed as stripes are analyzed, using a principle of microwave reflection, as well as an extraction method for each reflection return, and a parameter retrieval method is obtained. Finally, the parameters of the bridge, including the height （top and bottom surfaces of the sea bridge）, width, thickness, span, and height of the bridge tower, are retrieved using full polarimetric AIRSAR data. When a comparison of the measured data is completed, the results indicate that the proposed method can invert the parameters with a high accuracy, and that the inversion error of the bridge height （bottom surface） is only 1.3%. Moreover, the results also show that for the high-resolution SAR, the C and L-band images have the same ability in regards to parameter retrieval.

A fast inversion method for ocean parameters based on dispersion curves with a single hydrophone

The dispersion characteristics of shallow water can be described by the dispersion curves,which contain substantial ocean parameter information.A fast ocean parameter inversion method based on dispersion curves with a single hydrophone is presented in this paper.The method is achieved through Bayesian theory.Several sets of dispersion curves extracted from measured data are used as the input function.The inversion is performed by matching a replica calculated with a dispersion formula.The bottom characteristics can be described by the bottom reflection phase shift parameter P.The propagation range and the depth can be inverted quickly when the seabed parameters are represented by on parameter P.The inversion results improve the inversion efficiency of the seabed parameters.Consequently,the inversion efficiency and accuracy are improved while the number of inversion parameters is decreased and the computational speed of replica is increased.The inversion results have lower error than the reference values,and the dispersion curves calculated with inversion parameters are also in good agreement with extracted curves from measured data;thus,the effectiveness of the inversion method is demonstrated.

PARAMETERS INVERSION OF FLUID-SATURATED POROUS MEDIA BASED ON NEURAL NETWORKS

The multi- layers feedforward neural network is used for inversion ofmaterial constants of fluid-saturated porous media. The direct analysis of fluid-saturated porousmedia is carried out with the boundary element method. The dynamic displacement responses obtainedfrom direct analysis for prescribed material parameters constitute the sample sets training neuralnetwork. By virtue of the effective L-M training algorithm and the Tikhonov regularization method aswell as the GCV method for an appropriate selection of regu-larization parameter, the inversemapping from dynamic displacement responses to material constants is performed. Numerical examplesdemonstrate the validity of the neural network method.

DBN-GABP model for estimation of aircraft wake vortex parameters using Lidar data

Aircraft wake turbulence is an inherent outcome of aircraft flight,presenting a substan-tial challenge to air traffic control,aviation safety and operational efficiency.Building upon data obtained from coherent Doppler Lidar detection,and combining Dynamic Bayesian Networks(DBN)with Genetic Algorithm-optimized Backpropagation Neural Networks(GA-BPNN),this paper proposes a model for the inversion of wake vortex parameters.During the wake vortex flow field simulation analysis,the wind and turbulent environment were initially superimposed onto the simulated wake velocity field.Subsequently,Lidar-detected echoes of the velocity field are simulated to obtain a data set similar to the actual situation for model training.In the case study validation,real measured data underwent preprocessing and were then input into the established model.This allowed us to construct the wake vortex characteristic parameter inversion model.The final results demonstrated that our model achieved parameter inversion with only minor errors.In a practical example,our model in this paper significantly reduced the mean square error of the inverted velocity field when compared to the traditional algorithm.This study holds significant promise for real-time monitoring of wake vortices at airports,and is proved a crucial step in developing wake vortex interval standards.

Parameter sensitivity and inversion analysis of a concrete faced rock-fill dam based on HS-BPNN algorithm

Considering the complex nonlinear relationship between the material parameters of a concrete faced rock-fill dam(CFRD) and its displacements, the harmony search(HS) algorithm is used to optimize the back propagation neural network(BPNN), and the HS-BPNN algorithm is formed and applied for the inversion analysis of the parameters of rock-fill materials. The sensitivity of the parameters in the Duncan and Chang's E-B model is analyzed using the orthogonal test design. The case study shows that the parameters φ0, K, Rf, and Kb are sensitive to the deformation of the rock-fill dam and the inversion analysis for these parameters is performed by the HS-BPNN algorithm. Compared with the traditional BPNN, the HS-BPNN algorithm exhibits the advantages of high convergence precision, fast convergence rate, and strong stability.

A topographic parameter inversion method based on laser altimetry

A topographic parameter inversion method based on laser altimetry is developed in this paper, which can be used to deduce the surface vertical profile and retrieve the topographic parameters within the laser footprints by analyzing and simulating return waveforms. This method comprises three steps. The first step is to build the numerical models for the whole measuring procedure of laser altimetry, construct digital elevation models for surfaces with different topographic parameters, and calculate return waveforms. The second step is to analyze the simulated return waveforms to obtain their characteristics parameters, summarize the effects of the topographic parameter variations on the characteristic parameters of simulated return waveforms, and analyze the observed return waveforms of laser altimeters to acquire their characteristic parameters at the same time. The last step is to match the characteristic parameters of the simulated and observed return waveforms, and deduce the topographic parameters within the laser footprint. This method can be used to retrieve the topographic parameters within the laser footprint from the observed return waveforms of spaceborne laser altimeters and to get knowledge about the surface altitude distribution within the laser footprint other than only getting the height of the surface encountered firstly by the laser beam, which extends laser altimeters' function and makes them more like radars.

Wear monitoring method of water-lubricated polymer thrust bearing based on ultrasonic reflection coefficient amplitude spectrum

The water-lubricated thrust bearings of the marine rim-driven thruster(RDT)are usually composed of polymer composites,which are prone to serious wear under harsh working conditions.Ultrasonic is an excellent non-destructive monitoring technology,but polymer materials are characterized by viscoelasticity,heterogeneity,and large acoustic attenuation,making it challenging to extract ultrasonic echo signals.Therefore,this paper proposes a wear monitoring method based on the amplitude spectrum of the ultrasonic reflection coefficient.The effects of bearing parameters,objective function,and algorithm parameters on the identification results are simulated and analyzed.Taking the correlation coefficient and root mean square error as the matching parameters,the thickness,sound velocity,density,and attenuation factor of the bearing are inversed simultaneously by utilizing the differential evolution algorithm(DEA),and the wear measurement system is constructed.In order to verify the identification accuracy of this method,an accelerated wear test under heavy load was executed on a multi-functional vertical water lubrication test rig with poly-ether-etherketone(PEEK)fixed pad and stainless-steel thrust collar as the object.The thickness of pad was measured using the high-precision spiral micrometer and ultrasonic testing system,respectively.Ultimately,the results demonstrate that the thickness identification error of this method is approximately 1%,and in-situ monitoring ability will be realized in the future,which is of great significance to the life prediction of bearings.

Inversion problem for the dimension of fractal rough surface

In the present paper, the fractal rough surface is described by a band-limited Weierstrass-Mandelbrot function. By using the Monte Carlo method and optimal method, a minimal target function method is applied to inverting the fractal dimension of the fractal rough surface. Numerical simulations show that the method can avoid the influence of the fractal characteristic scale, and that the method is of high precision.