Experimental Research on Quantitative Inversion Models of Suspended Sediment Concentration Using Remote Sensing Technology

Research on quantitative models of suspended sediment concentration （SSC） using remote sensing technology is very important to understand the scouting and siltation variation in harbors and water channels. Based onlaboratory study of the relationship between different suspended sediment concentrations and reflectance spectra measured synchronously, quantitative inversion models of SSC based on single factor, band ratio and sediment parameter were developed, which provides an effective method to retrieve the SSC from satellite images. Results show that the bl （430-500nm） and b3 （670-735nm） are the optimal wavelengths for the estimation of lower SSC and the b4 （780-835nm） is the optimal wavelength to estimate the higher SSC. Furthermore the band ratio B2/B3 can be used to simulate the variation of lower SSC better and the B4/B1 to estimate the higher SSC accurately. Also the inversion models developed by sediment parameters of higher and lower SSCs can get a relatively higher accuracy than the single factor and band ratio models.

Metaheuristic algorithms for groundwater model parameter inversion:Advances and prospects

Groundwater inverse modeling is a vital technique for estimating unmeasurable model parameters and enhancing numerical simulation accuracy.This paper comprehensively reviews the current advances and future prospects of metaheuristic algorithm-based groundwater model parameter inversion.Initially,the simulation-optimization parameter estimation framework is introduced,which involves the integration of simulation models with metaheuristic algorithms.The subsequent sections explore the fundamental principles of four widely employed metaheuristic algorithms-genetic algorithm(GA),particle swarm optimization(PSO),simulated annealing(SA),and differential evolution(DE)-highlighting their recent applications in water resources research and related areas.Then,a solute transport model is designed to illustrate how to apply and evaluate these four optimization algorithms in addressing challenges related to model parameter inversion.Finally,three noteworthy directions are presented to address the common challenges among current studies,including balancing the diverse exploration and centralized exploitation within metaheuristic algorithms,local approxi-mate error of the surrogate model,and the curse of dimensionality in spatial variational heterogeneous pa-rameters.In summary,this review paper provides theoretical insights and practical guidance for further advancements in groundwater inverse modeling studies.

The sensitivity of the array resistivity log to mud filtrate invasion and its primary five-parameter inversion for improved oil water recognition

In order to improve reservoir fluid recognition, the sensitivity of array resistivity response to the difference of the invasion properties in both oil-bearing layers and water layers is analyzed. Then the primary inversion is carried out based on the array resistivity log. The mud invasion process is numerically simulated based on the oil-water flow equation and water convection diffusion equation. The results show that the radial resistivity of a fresh mud-invaded oil-bearing layer presents complex distribution characteristics, such as nonlinear increase, increasing to decreasing and low resistivity annulus, and the resistive invasion profile of a water layer is monotonic. Under specific conditions, array resistivity log can reflect these changes and the array induction log is more sensitive. Nevertheless, due to the effect of factors like large invasion depth, reservoir physical and oil-bearing properties, the measured apparent resistivity may differ greatly from the actual mud filtrate invasion profile in an oil-bearing layer. We proposed a five-parameter formation model to simulate the complex resistivity distribution of fresh mud-invaded formation. Then, based on the principle of non-linear least squares, the measured array resistivity log is used for inversion with the Marquardt method. It is demonstrated that the inverted resistivity is typically non-monotonic in oil-bearing layers and is monotonic in water layers. Processing of some field data shows that this is helpful in achieving efficient reservoir fluid recognition.

Seismic traveltime inversion of 3D velocity model with triangulated interfaces

Seismic traveltime tomographic inversion has played an important role in detecting the internal structure of the solid earth. We use a set of blocks to approximate geologically complex media that cannot be well described by layered models or cells. The geological body is described as an aggregate of arbitrarily shaped blocks,which are separated by triangulated interfaces. We can describe the media as homogenous or heterogeneous in each block. We define the velocities at the given rectangle grid points for each block,and the heterogeneous velocities in each block can be calculated by a linear interpolation algorithm. The parameters of the velocity grid positions are independent of the model parameterization,which is advantageous in the joint inversion of the velocities and the node depths of an interface. We implement a segmentally iterative ray tracer to calculate traveltimes in the 3D heterogeneous block models.The damped least squares method is employed in seismic traveltime inversion,which includes the partial derivatives of traveltime with respect to the depths of nodes in the triangulated interfaces and velocities defined in rectangular grids. The numerical tests indicate that the node depths of a triangulated interface and homogeneous velocity distributions can be well inverted in a stratified model.

3-D multi-parameter type traveltime tomography in a spherical coordinate frame: comparison of double and triple class simultaneous inversions

It is now common practice to perform simultaneous traveltime inversion for the velocity field and the reflector geometry in reflection/refraction tomography, or the velocity field and the hypocenter locations in regional earthquake tomography, but seldom are all three classes of model parameters updated simultaneously. This is mainly due to the trade-off between the different types of model parameters and the lack of different seismic phases to constrain the model parameters. Using a spherical-coordinate ray tracing algorithm for first and later（primary reflected） arrival tracing algorithm in combination with a popular linearized inversion solver, it is possible to simultaneously recover the three classes of model parameters in regional or global tomographic studies. In this paper we incorporate the multistage irregular shortest-path ray tracing algorithm（in a spherical coordinate system） with a subspace inversion solver to formulate a simultaneous inversion algorithm for triple model parameters updating using direct and later arrival time information.Comparison tests for two sets of data（noise free and added noise） indicate that the new triple-class parameter inversion algorithm is capable of obtaining nearly the same results as the double-class parameter inversion scheme. Furthermore,the proposed multi-parameter type inversion method is not sensitive to a modest level of picking error in the traveltime data, and also performs well with a relatively large uncertainty in earthquake hypocentral locations. This shows it to be a feasible and promising approach in regional or global tomographic applications.

Improving the Accuracy of Vegetation Index Retrieval for Biomass by Combining Ground-UAV Hyperspectral Data-A New Method for Inner Mongolia Typical Grasslands

Grassland biomass is an important parameter of grassland ecosystems.The complexity of the grassland canopy vegetation spectrum makes the long-term assessment of grassland growth a challenge.Few studies have explored the original spectral information of typical grasslands in Inner Mongolia and examined the influence of spectral information on aboveground biomass(AGB)estimation.In order to improve the accuracy of vegetation index inversion of grassland AGB,this study combined ground and Unmanned Aerial Vehicle(UAV)remote sensing technology and screened sensitive bands through ground hyperspectral data transformation and correlation analysis.The narrow band vegetation indices were calculated,and ground and airborne hyperspectral inversion models were established.Finally,the accuracy of the model was verified.The results showed that:(1)The vegetation indices constructed based on the ASD FieldSpec 4 and the UAV were significantly correlated with the dry and fresh weight of AGB.(2)The comparison between measured R^(2) with the prediction R^(2) indicated that the accuracy of the model was the best when using the Soil-Adjusted Vegetation Index(SAVI)as the independent variable in the analysis of AGB(fresh weight/dry weight)and four narrow-band vegetation indices.The SAVI vegetation index showed better applicability for biomass monitoring in typical grassland areas of Inner Mongolia.(3)The obtained ground and airborne hyperspectral data with the optimal vegetation index suggested that the dry weight of AGB has the best fitting effect with airborne hyperspectral data,where y=17.962e^(4.672x),the fitting R^(2) was 0.542,the prediction R^(2)was 0.424,and RMSE and REE were 57.03 and 0.65,respectively.Therefore,established vegetation indices by screening sensitive bands through hyperspectral feature analysis can significantly improve the inversion accuracy of typical grassland biomass in Inner Mongolia.Compared with ground monitoring,airborne hyperspectral monitoring better reflects the inversion of actual surface biomass.It provides a reliable modeling framework for grassland AGB monitoring and scientific and technological support for grazing management.

A Hybrid Compensation Scheme for the Input Rate-Dependent Hysteresis of the Piezoelectric Ceramic Actuators

A hybrid compensation scheme for piezoelectric ceramic actuators(PEAs)is proposed.In the hybrid compensation scheme,the input rate-dependent hysteresis characteristics of the PEAs are compensated.The feedforward controller is a novel input rate-dependent neural network hysteresis inverse model,while the feedback controller is a proportion integration differentiation(PID)controller.In the proposed inverse model,an input ratedependent auxiliary inverse operator(RAIO)and output of the hysteresis construct the expanded input space(EIS)of the inverse model which transforms the hysteresis inverse with multi-valued mapping into single-valued mapping,and the wiping-out,rate-dependent and continuous properties of the RAIO are analyzed in theories.Based on the EIS method,a hysteresis neural network inverse model,namely the dynamic back propagation neural network(DBPNN)model,is established.Moreover,a hybrid compensation scheme for the PEAs is designed to compensate for the hysteresis.Finally,the proposed method,the conventional PID controller and the hybrid controller with the modified input rate-dependent Prandtl-Ishlinskii(MRPI)model are all applied in the experimental platform.Experimental results show that the proposed method has obvious superiorities in the performance of the system.

Spectrum Characteristics of Cotton Canopy Infected with Verticillium Wilt and Applications

Hyper spectrum remote sensing with fine spectrum information is an efficient method to estimate the verticillium wilt of cotton. The research was conducted in Xinjiang, the largest cotton plant region of China, by using the data which were collected both by canopy spectrum infected with verticillium wilt and severity level （SL） in the year 2005-2006. The quantitative correlation was analyzed between SL and canopy of reflectance spectrum or derivative spectrum reflectance. The results indicated that spectrum characteristics of cotton canopy infected with verticillium wilt changed regularly with the increase of SL in different periods and varieties, Spectrum reflectance increased in the visible light region （620-700 nm） with the increase of the SL, which inverted in near-infrared region and was extremely significant in the region of （780-1 300 nm）. When SL attained b2 （DI = 25）, cotton canopy infected with verticillium wilt was used as a watershed and diagnosed index in the beginning stages of the disease. The results also indicated that there were marked different characteristics of the first derivative spectrum in these SL, it changed significantly in the red edge ranges （680-760 nm） with different SL, i.e., red edge swing decreased, and red edge position equally moved to the blue. In this study 1 001-1 110 nm and 1 205- 1 320 nm were selected out as sensitive bands for SL of canopy. Inversion models established for estimating cotton canopy infected with verticillium wilt reached the most significant level. Finally, the different spectrum characteristics of cotton canopy infected with verticillium wilt were marked, some inversion models were established, which could estimate SL of canopy infected with verticillium wilt. The best recognized model was the first derivative spectra at （FD 731 nm- FD 1317 nm）, and it might be used to forecast the position of cotton canopy infected with verticillium wilt quantitatively.

Genetic Algorithm-Based Estimation of Nonlinear Transducer

This paper describes an innovative, genetic algorithm based inverse model of nonlinear transducer. In the inverse modeling, using a genetic algorithm, the unknown coefficients of the model are estimated accurately. The simulation results indicate that this technique provides greater flexibility and suitability than the existing methods. It is very easy to modify the nonlinear transducer on line. Thus the method improves the transducer's accuracy. With the help of genetic algorithm (GA), the model coefficients' training are less likely to be trapped in local minima than traditional gradient based search algorithms.

Similarity measure of sedimentary successions and its application in inverse stratigraphic modeling

This paper presents a unique and formal method of quantifying the similarity or distance between sedimentary facies successions from measured sections in outcrop or drilled wells and demonstrates its first application in inverse stratigraphic modeling. A sedimentary facies succession is represented with a string of symbols, or facies codes in its natural vertical order, in which each symbol brings with it one attribute such as thickness for the facies. These strings are called attributed strings. A similarity measure is defined between the attributed strings based on a syntactic pattern-recognition technique. A dynamic programming algorithm is used to calculate the similarity. Inverse stratigraphic modeling aims to generate quantitative 3D facies models based on forward stratigraphic modeling that honors observed datasets. One of the key techniques in inverse stratigraphic modeling is how to quantify the similarity or distance between simulated and observed sedimentary facies successions at data locations in order for the forward model to condition the simulation results to the observed dataset such as measured sections or drilled wells. This quantification technique comparing sedimentary successions is demonstrated in the form of a cost function based on the defined distance in our inverse stratigraphic modeling implemented with forward modeling optimization.

Retrieval of leaf biochemical properties by inversed PROSPECT model and hyperspectral indices:an application to Populus euphratica polymorphic leaves

Leaf biochemical properties have been widely assessed using hyperspectral reflectance information by inversion of PROSPECT model or by using hyperspectral indices, but few studies have focused on arid ecosystems. As a dominant species of riparian ecosystems in arid lands, Populus euphratica Oliv. is an unusual tree species with polymorphic leaves along the vertical profile of canopy corresponding to different growth stages. In this study, we evaluated both the inversed PROSPECT model and hyperspectral indices for estimating biochemical properties of P. euphratica leaves. Both the shapes and biochemical properties of P. euphratica leaves were found to change with the heights from ground surface. The results indicated that the model inversion calibrated for each leaf shape performed much better than the model calibrated for all leaf shapes, and also better than hyperspectral indices. Similar results were obtained for estimations of equivalent water thickness （EWT） and leaf mass per area （LMA）. Hyperspectral indices identified in this study for estimating these leaf properties had root mean square error （RMSE） and R2 values between those obtained with the two calibration strategies using the inversed PROSPECT model. Hence, the inversed PROSPECT model can be applied to estimate leaf biochemical properties in arid ecosystems, but the calibration to the model requires special attention.

Algebraic Method‑Based Point‑to‑Point Trajectory Planning of an Under‑Constrained Cable‑Suspended Parallel Robot with Variable Angle and Height Cable Mast

To avoid impacts and vibrations during the processes of acceleration and deceleration while possessing flexible working ways for cable-suspended parallel robots(CSPRs),point-to-point trajectory planning demands an under-constrained cable-suspended parallel robot(UCPR)with variable angle and height cable mast as described in this paper.The end-effector of the UCPR with three cables can achieve three translational degrees of freedom(DOFs).The inverse kinematic and dynamic modeling of the UCPR considering the angle and height of cable mast are completed.The motion trajectory of the end-effector comprising six segments is given.The connection points of the trajectory segments(except for point P3 in the X direction)are devised to have zero instantaneous velocities,which ensure that the acceleration has continuity and the planned acceleration curve achieves smooth transition.The trajectory is respectively planned using three algebraic methods,including fifth degree polynomial,cycloid trajectory,and double-S velocity curve.The results indicate that the trajectory planned by fifth degree polynomial method is much closer to the given trajectory of the end-effector.Numerical simulation and experiments are accomplished for the given trajectory based on fifth degree polynomial planning.At the points where the velocity suddenly changes,the length and tension variation curves of the planned and unplanned three cables are compared and analyzed.The OptiTrack motion capture system is adopted to track the end-effector of the UCPR during the experiment.The effectiveness and feasibility of fifth degree polynomial planning are validated.

Coseismic and postseismic slip ruptures for 2015Mw 6.4 Pishan earthquake constrained by static GPS solutions

On 3 July 2015, a Mw 6.4 earthquake occurred on a blind fault struck Pishan, Xinjiang,China. By combining Crustal Movement Observation Network of China（CMONOC） and other Static Global Positioning System（GPS） sites surrounding Pishan region, it provides a rare chance for us to constrain the slip rupture for such a moderate event. The maximum displacement is up to 12 cm, 2 cm for coseismic and postseismic deformation, respectively,and both the deformation patterns show a same direction moving northeastward. With rectangular dislocation model, a magnitude of Mw6.48, Mw6.3 is calculated based on coseismic, postseismic deformation respectively. Our result indicates the western Kunlun range is still moving toward Tarim Basin followed by an obvious postseismic slip associated with this earthquake. To determine a more reasonable model for postseismic deformation, a longer GPS dataset will be needed.

Direct shear tests on cemented paste backfill-rock wall and cemented paste backfill-backfill interfaces

Even though a large number of large-scale arch dams with height larger than 200 m have been built in the world, the transient groundwater flow behaviors and the seepage control effects in the dam foundations under difficult geological conditions are rarely reported. This paper presents a case study on the transient groundwater flow behaviors in the rock foundation of Jinping I double-curvature arch dam, the world＇s highest dam of this type to date that has been completed. Taking into account the geological settings at the site, an inverse modeling technique utilizing the time series measurements of both hydraulic head and discharge was adopted to back-calculate the permeability of the foundation rocks,which effectively improves the uniqueness and reliability of the inverse modeling results. The transient seepage flow in the dam foundation during the reservoir impounding was then modeled with a parabolic variational inequality（PVI） method. The distribution of pore water pressure, the amount of leakage, and the performance of the seepage control system in the dam foundation during the entire impounding process were finally illustrated with the numerical results.

A new method for multi-point pollution source identification

Rapid and accurate identification of the characteristics(source location,number,and intensity)of pollution sources is essential for emergency assessment of contamination events.Compared with single-point source iden-tification,the reconstruction of multiple sources is more challenging.In this study,a two-step inversion method is proposed for multi-point pollution source reconstruction from limited measurements with the number of sources unknown.The applicability of the proposed method is validated with a set of synthetic experiments correspond-ing to one-,two-,and three-point pollution sources.The results show that the number and locations of pollution sources are retrieved exactly the same as prescribed,and the source intensities are estimated with negligible errors.The algorithm exhibits good performance in single-and multi-point pollution source identification,and its accuracy and efficiency of identification do not deteriorate with the increase in the number of sources.Some limitations of the algorithm,together with its capabilities,are also discussed in this paper.

Quantitative description of infrared radiation characteristics of preflawed sandstone during fracturing process

Previous studies show that infrared radiation temperature(IRT)abnormalities are always accompanied by the crack development in rocks under external loads.In this context,experiments were conducted on preflawed sandstone to investigate the infrared radiation characteristics during failure process.Two indicators were defined herein,i.e.coefficient of variation of IRT(CVIRT)and skewness of IRT(SIRT).The regression analysis shows that the IRT probability distributions during loading process fit the Gaussian model.The variations in the CVIRT are characterized by four stages:primary stage,steady stage,accelerating stage and post-peak stage.Besides,the variations in the SIRT are divided into three stages:primary stage,steady stage and failure and post-peak stage.The precursor point for preflawed rock failure is identified based on the CVIRTetime curve,with average precursor point of 83%of the peak stress.Compared with other IRT indicators,the proposed two IRT indicators have higher sensitivity to IRT abnormalities during failure process.Furthermore,the connection between the IRT indicators and the rock fracturing was investigated to interpret the IRT indicator abnormalities.Based on the Verhulst inverse function,a new quantitative model was presented to describe the primary stage,steady stage and accelerating stage of the CVIRTetime curve.The results obtained in this study can provide early-warning information for rock failure prediction.

Geophysical and Production Data History Matching Based on Ensemble Smoother with Multiple Data Assimilation

The Ensemble Kalman Filter(EnKF),as the most popular sequential data assimilation algorithm for history matching,has the intrinsic problem of high computational cost and the potential inconsistency of state variables updated at each loop of data assimilation and its corresponding reservoir simulated result.This problem forbids the reservoir engineers to make the best use of the 4D seismic data,which provides valuable information about the fluid change inside the reservoir.Moreover,only matching the production data in the past is not enough to accurately forecast the future,and the development plan based on the false forecast is very likely to be suboptimal.To solve this problem,we developed a workflow for geophysical and production data history matching by modifying ensemble smoother with multiple data assimilation(ESMDA).In this work,we derived the mathematical expressions of ESMDA and discussed its scope of applications.The geophysical data we used is P-wave impedance,which is typically included in a basic seismic interpretation,and it directly reflects the saturation change in the reservoir.Full resolution of the seismic data is not necessary,we subsampled the P-wave impedance data to further reduce the computational cost.With our case studies on a benchmark synthetic reservoir model,we also showed the supremacy of matching both geophysical and production data,than the traditional reservoir history matching merely on the production data:the overall percentage error of the observed data is halved,and the variances of the updated forecasts are reduced by two orders of the magnitude.

Study of Spectral Response Characteristics of Oilseed Rape (Brassica napus) to Particulate Matters Based on Hyper-Spectral Technique

Haze is mainly caused by the suspended particulate matters in the air,of which the particulate matters pollution harms leaf vegetables.In this paper,oilseed rapes at four different growing periods were investigated in a simulated particulate pollution environment.In combination of hyper-spectral technology and micro examination,the response of hyper-spectral characteristics of the leaf to particulate matters was investigated in-depth.The hyperspectral,chlorophyll content,net photosynthetic rate and stomatal conductance of leaf were obtained.The deposition and adsorption of particulate matters on the leaf were observed by Environmental Scanning Electron Microscope(ESEM).Normalized difference vegetation index(NDVI),modified red edge normalized(mNDVI705)and modified red edge simple ratio index(mSR705)were selected as characteristic parameters and the range of 510 nm~620 nm as the sensitive band.16 methods were used to establish the physiological information inversion model.The main results were as follows:Under the influence of particulate matters,the spectral reflectance decreased as a whole.With the increase of leaf age,the phenomenon of blue shift aggravated.The amplitude of yellow and blue edge decreased with overall decreasing vegetation indices.The furrows and irregular band protrusions in leaves were favorable for keeping particulate matters.With longer affecting time and more deposition of particle matters on the leaf,the stomatal opening became smaller.After comparing,principal component regression(PCR)+multiple scatter correction(MSC)+second derivative(SD)+Savitzky-Golay smooth(SG),and partial least square(PLS)+multiple scatter correction(MSC)+first derivative(FD)+Savitzky-Golay smooth(SG)were determined the best method to establish the inversion model of chlorophyll content and net photosynthetic rate respectively.This study may bring novel ideas for the diagnosis and analysis of the physiological response of leaf vegetables under particulate matters pollution using hyper-spectral technology.

Inverse Calculation of Wave-Absorbing Structure Dimensions Based on Extended ANFIS Model

A new wave energy dissipation structure is proposed, aiming to optimize the dimensions of the structure and make the reflection of the structure maintain a low level within the scope of the known frequency band. An optimal extended ANFIS model combined with the wave reflection coefficient analysis for the estimation of the structure dimensions is established. In the premise of lower wave reflection coefficient, the specific sizes of the structure are obtained inversely, and the contribution of each related parameter on the structural reflection performance is analyzed. The main influencing factors are determined. It is found that the optimal dimensions of the proposed structure exist, which make the wave absorbing performance of the structure reach a perfect status under a wide wave frequency band.

Calculation and analysis of losses of magnetic-valve controllable reactor

Magnetic-valve controllable reactor(MCR)has characteristics of DC bias and different types of magnetic flux density in the magnetic circuit and winding current distortion.These characteristics not only lead to loss calculation method of MCR different from that of power transformer,but also make it more difficult to calculate the core loss and wingding loss of MCR accurately.Our study combines core partition method with dynamic inverse J-A model to calculate the core loss of MCR.The winding loss coefficient of MCR is proposed,which takes into account the influence of harmonics and magnetic flux leakage on the winding loss of MCR.The result shows that the proposed core loss calculation method and winding loss coefficient are effective and correct for the loss calculation of MCR.