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.

Robust inversion analysis of local gravity anomalies caused by geological dislocation model of faults

Theoretical analysis and practical observations show that fault dislocations can change the gravity field around the fault. Gravity changes which were caused by the repeated dislocations over a long period of time were superimposed on the Bougeur gravity anomalies. These anomalies became the evidence of historical movement of fault as well as provide a way for the study of paleo earthquakes. This paper investigates inversion methods for the geological dislocation modeling of faults using the local Bouguer's gravity anomalies. To remove the effects of the irrelevant part of gravity anomalies to fault movements, we propose the robust nonlinear inversion method and set up the corresponding algorithm. Modeling examples indicate that the Marquardt's and Baye's least squares solutions depart from the true solution due to the attraction of gross errors in the data. The more seriously the data is contaminated, the more seriously the solutions are biased. In contrast, the proposed robust Marquardt's and Baye's inversion solutions can still maintain consistency with the solution without gross errors, even though 50 percent of the data is contaminated. This indicates that the proposed robust methods are effective. Using the proposed methods, we invert the geological dislocation models of the faults around the Erhai Lake in West Yunnan. The results show that the Northern Cangdong fault and the Erhai fault are normal dip slip faults with about 4 to 5 km dislocations; and that the Southern Cangdong fault has a less dip slip compared with the former two. A satisfactory fitting between the theoretical values of the inversion solution and the actual local gravity field is achievable.

Model fitting of the seasonal population dynamics of the soybean aphid, Aphis glycines Matsumura, in the field

The soybean aphid, Aphis glycines Matsumura(Hemiptera: Aphididae), is one of the greatest threats to soybean production, and both trend analysis and periodic analysis of its population dynamics are important for integrated pest management(IPM). Based on systematically investigating soybean aphid populations in the field from 2018 to 2020, this study adopted the inverse logistic model for the first time, and combined it with the classical logistic model to describe the changes in seasonal population abundance from colonization to extinction in the field. Then, the increasing and decreasing phases of the population fluctuation were divided by calculating the inflection points of the models, which exhibited distinct seasonal trends of the soybean aphid populations in each year. In addition, multifactor logistic models were then established for the first time, in which the abundance of soybean aphids in the field changed with time and relevant environmental conditions. This model enabled the prediction of instantaneous aphid abundance at a given time based on relevant meteorological data. Taken as a whole, the successful approaches implemented in this study could be used to build a theoretical framework for practical IPM strategies for controlling soybean aphids.

Machine Learning-based Inverse Model for Few-Mode Fiber Designs

The medium for next-generation communication is considered as fiber for fast,secure communication and switching capability.Mode division and space division multiplexing provide an excellent switching capability with high data transmission rate.In this work,the authors have approached an inverse modeling technique using regression-based machine learning to design a weakly coupled few-mode fiber for facilitating mode division multiplexing.The technique is adapted to predict the accurate profile parameters for the proposed few-mode fiber to obtain the maximum number of modes.It is for a three-ring-core few-mode fiber for guiding five,ten,fifteen,and twenty modes.Three types of regression models namely ordinary least-square linear multi-output regression,k-nearest neighbors of multi-output regression,and ID3 algorithm-based decision trees for multi-output regression are used for predicting the multiple profile parameters.It is observed that the ID3-based decision tree for multioutput regression is the robust,highly-accurate machine learning model for fast modeling of FMFs.The proposed fiber claims to be an efficient candidate for the next-generation 5G and 6G backhaul networks using mode division multiplexing.

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.

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.

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.

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.

Nonlinear inverse modeling of sensor based on back-propagation fuzzy logical system

Objective To correct the nonlinear error of sensor output,a new approach to sensor inverse modeling based on Back-Propagation Fuzzy Logical System(BP FS) is presented.Methods The BP FS is a computationally efficient nonlinear universal approximator,which is capable of implementing complex nonlinear mapping from its input pattern space to the output with fast convergence speed.Results The neuro-fuzzy hybrid system,i.e.BP FS,is then applied to construct nonlinear inverse model of pressure sensor.The experimental results show that the proposed inverse modeling method automatically compensates the associated nonlinear error in pressure estimation,and thus the performance of pressure sensor is significantly improved.Conclusion The proposed method can be widely used in nonlinearity correction of various kinds of sensors to compensate the effects of nonlinearity and temperature on sensor output.

Inversion Methods of Optical Constants of Semitransparent Solid Materials from Transmittance Spectrograms

The direct calculation models of spectral transmittance of single and double slabs consisted of semitransparent solid materials were developed based on ray trace method,and a new inversion method of optical constants ( k is extinction coefficient and n is refractive index) of materials was proposed based on transmittance spectrograms of double slabs. Differences between the new method and two others currently used methods were studied,and application range of methods was also investigated. Optical constants of selenide glass attained in references were selected as true values,and spectral transmittances of glass simulated based on direct calculation model were regarded as experimental values. Optical constants of selenide glass were achieved by inverse models. Influences of measurement error on inverse results were also determined. The results showed that: (1) based on transmittance spectrograms of double slabs in which thickness of single slab is the same,the new proposed method can attain optical constants of materials; (2) the effect of optical constants n and k on three inversion methods are urgent larger,but inversed calculation precision of optical constants are higher in most application ranges; (3) the influence of measurement errors existed in experimental datum on the inverse precision of three methods are urgent

Model predictive inverse method for recovering boundary conditions of two-dimensional ablation

A model predictive inverse method (MPIM) is presented to estimate the time- and space-dependent heat flux onthe ablated boundary and the ablation velocity of the two-dimensional ablation system. For the method, first of all, therelationship between the heat flux and the temperatures of the measurement points inside the ablation material is establishedby the predictive model based on an influence relationship matrix. Meanwhile, the estimation task is formulated as aninverse heat transfer problem (IHTP) with consideration of ablation, which is described by an objective function of thetemperatures at the measurement point. Then, the rolling optimization is used to solve the IHTP to online estimate theunknown heat flux on the ablated boundary. Furthermore, the movement law of the ablated boundary is reconstructedaccording to the estimation of the boundary heat flux. The effects of the temperature measurement errors, the numberof future time steps, and the arrangement of the measurement points on the estimation results are analyzed in numericalexperiments. On the basis of the numerical results, the effectiveness of the presented method is clarified.

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.

Geochemical Modeling and Statistical Analysis for Groundwater Evolution Assessment in Wadi Qasab, Sohag, Eastern Desert, Egypt

Qasab basin is one of the most promising areas for the sustainable development in the Eastern Desert fringes of the Nile Valley, Egypt. The integration between statistical analysis, stable isotopes as well as geochemical modeling tools delineated the geochemical possesses affecting groundwater quality and detected the main recharge source in Qasab basin. The most of groundwater samples are brackish (88%), while the minority (12%) of the samples are fresh. The electrical conductivity of groundwater ranged from 1135 to 10,030 μS/cm. The statistical analysis and hydrochemical diagrams suggest that the groundwater quality is mainly controlled by several intermixed processes (rock weathering and agricultural activities). The mineralization of the Pleistocene groundwater is regulated by the rock weathering source, evaporation processes and reverse cation exchange. The isotopic signatures (δ2H and δ18O) represent two groundwater groups. The first group, is enriched with the isotopic signature of δ18O, which ranges from 0.9‰ to 5.5‰. This group is mostly affected by the recent meteoric recharge from the surface water leakage. The second group, is relatively depleted with the isotopic signature of δ18O, reflecting a palaeo recharge source of colder climate. The δ18O‰ varies from -10.1‰ to -6.4‰, indicating upward leakage of the Nubian sandstone aquifer through deep seated faults. The inverse geochemical model reflects that the salinity source of the groundwater samples is due to the leaching and dissolution processes of carbonate, sulphate and chloride minerals from the aquifer matrix. This study can demonstrate the hydrochemistry assessment guide to support sustainable development in Qasab basin to ensure that adequate groundwater management can play to reduce poverty and support socioeconomic development.

Estimation of the water–oil–gas relative permeability curve from immiscible WAG coreflood experiments using the cubic B-spline model

Immiscible water-alternating-gas（WAG） flooding is an EOR technique that has proven successful for water drive reservoirs due to its ability to improve displacement and sweep efficiency.Nevertheless,considering the complicated phase behavior and various multiphase flow characteristics,gas tends to break through early in production wells in heterogeneous formations because of overriding,fingering,and channeling,which may result in unfavorable recovery performance.On the basis of phase behavior studies,minimum miscibility pressure measurements,and immiscible WAG coreflood experiments,the cubic B-spline model（CBM） was employed to describe the three-phase relative permeability curve.Using the Levenberg-Marquardt algorithm to adjust the vector of unknown model parameters of the CBM sequentially,optimization of production performance including pressure drop,water cut,and the cumulative gas-oil ratio was performed.A novel numerical inversion method was established for estimation of the water-oil-gas relative permeability curve during the immiscible WAG process.Based on the quantitative characterization of major recovery mechanisms,the proposed method was validated by interpreting coreflood data of the immiscible WAG experiment.The proposed method is reliable and can meet engineering requirements.It provides a basic calculation theory for implicit estimation of oil-water-gas relative permeability curve.

Nonlinear direct data-driven control for UAV formation flight system

This paper proposes the nonlinear direct data-driven control from theoretical analysis and practical engineering,i.e.,unmanned aerial vehicle(UAV)formation flight system.Firstly,from the theoretical point of view,consider one nonlinear closedloop system with a nonlinear plant and nonlinear feed-forward controller simultaneously.To avoid the complex identification process for that nonlinear plant,a nonlinear direct data-driven control strategy is proposed to design that nonlinear feed-forward controller only through the input-output measured data sequence directly,whose detailed explicit forms are model inverse method and approximated analysis method.Secondly,from the practical point of view,after reviewing the UAV formation flight system,nonlinear direct data-driven control is applied in designing the formation controller,so that the followers can track the leader’s desired trajectory during one small time instant only through solving one data fitting problem.Since most natural phenomena have nonlinear properties,the direct method must be the better one.Corresponding system identification and control algorithms are required to be proposed for those nonlinear systems,and the direct nonlinear controller design is the purpose of this paper.

Dynamic ocean inverse modeling based on differentiable rendering

Learning and inferring underlying motion patterns of captured 2D scenes and then re-creating dynamic evolution consistent with the real-world natural phenomena have high appeal for graphics and animation.To bridge the technical gap between virtual and real environments,we focus on the inverse modeling and reconstruction of visually consistent and property-verifiable oceans,taking advantage of deep learning and differentiable physics to learn geometry and constitute waves in a self-supervised manner.First,we infer hierarchical geometry using two networks,which are optimized via the differentiable renderer.We extract wave components from the sequence of inferred geometry through a network equipped with a differentiable ocean model.Then,ocean dynamics can be evolved using the reconstructed wave components.Through extensive experiments,we verify that our new method yields satisfactory results for both geometry reconstruction and wave estimation.Moreover,the new framework has the inverse modeling potential to facilitate a host of graphics applications,such as the rapid production of physically accurate scene animation and editing guided by real ocean scenes.

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.

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.

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.