Evaluation of unconstrained and constrained mathematical functions to model girth growth of rubber trees (Hevea brasiliensis) using young agemeasurements

No attempt has been made to date to model growth in girth of rubber tree （Hevea brasiliansis）. We evaluated the few widely used growth functions to identify the most parsimonious and biologically reasonable model for describing the girth growth of young rubber trees based on an incomplete set of young age measurements. Monthly data for girth of immature trees （age 2 to 12 yearsi from two locations were sub- jected to modelling. Re-parameterized, unconstrained and constrained growth functions,of Richards （RM）, Gompertz （GM） and the monomo- lecular ＇model ^（MM） were fitted to data. Duration of growth was the firsf constraint introduced. In the stagel We attempted a population aver- age （PA） model to capture the trend in growth. The best PA model was fitted as a subject specific （SS） model. We used appropriate error vari- ance-covariance structure to account for correlation due to repeated measurements over time. Unconstrainecl functions underestimated the asymptotic maximum that did not reflective carrying capacity of the locations. Underestimafions were attributed to the partial set＇ of meas- urements made during the early growth phase of the trees. MM proved superior to RM and GM. In the randomcoefficient models, both Gf and Go appeared to be influenced by tree level effects. Inclusion of diagonal definite positive matrix removed the correlation between random effects. The results were similar at both locations. In the overall assessment MM appeared as the candidate model for studying the girth-age relationships in Hevea trees. Based on the fitted model we conclude that, in Hevea trees, growth rate is maintained at maximum value at to, then decreases until the final state at dG/dt 〉 0, resulting in yield curve with no period of accelerating growth. One physiological explanation is that photosynthetic activity in Hevea trees decreases as girth increases and constructive metabolism is larger than destructive metabolism.

The distribution of deep source rocks in the GS Sag:joint MT-gravity modeling and constrained inversion

The coal-bearing strata of the deep Upper Paleozoic in the GS Sag have high hydrocarbon potential. Because of the absence of seismic data, we use electromagnetic （MT） and gravity data jointly to delineate the distribution of deep targets based on well logging and geological data. First, a preliminary geological model is established by using three-dimensional （3D） MT inversion results. Second, using the formation density and gravity anomalies, the preliminary geological model is modified by interactive inversion of the gravity data. Then, we conduct MT-constrained inversion based on the modified model to obtain an optimal geological model until the deviations at all stations are minimized. Finally, the geological model and a seismic profile in the middle of the sag is analysed. We determine that the deep reflections of the seismic profile correspond to the Upper Paleozoic that reaches thickness up to 800 m. The processing of field data suggests that the joint MT-gravity modeling and constrained inversion can reduce the multiple solutions for single geophysical data and thus improve the recognition of deep formations. The MT-constrained inversion is consistent with the geological features in the seismic section. This suggests that the joint MT and gravity modeling and constrained inversion can be used to delineate deep targets in similar basins.

A TRUST REGION METHOD WITH A CONIC MODEL FOR NONLINEARLY CONSTRAINED OPTIMIZATION

Trust region methods are powerful and effective optimization methods. The conic model method is a new type of method with more information available at each iteration than standard quadratic-based methods. The advantages of the above two methods can be combined to form a more powerful method for constrained optimization. The trust region subproblem of our method is to minimize a conic function subject to the linearized constraints and trust region bound. At the same time, the new algorithm still possesses robust global properties. The global convergence of the new algorithm under standard conditions is established.

Prediction of constrained modulus for granular soil using 3D discrete element method and convolutional neural networks

To efficiently predict the mechanical parameters of granular soil based on its random micro-structure,this study proposed a novel approach combining numerical simulation and machine learning algorithms.Initially,3500 simulations of one-dimensional compression tests on coarse-grained sand using the three-dimensional(3D)discrete element method(DEM)were conducted to construct a database.In this process,the positions of the particles were randomly altered,and the particle assemblages changed.Interestingly,besides confirming the influence of particle size distribution parameters,the stress-strain curves differed despite an identical gradation size statistic when the particle position varied.Subsequently,the obtained data were partitioned into training,validation,and testing datasets at a 7:2:1 ratio.To convert the DEM model into a multi-dimensional matrix that computers can recognize,the 3D DEM models were first sliced to extract multi-layer two-dimensional(2D)cross-sectional data.Redundant information was then eliminated via gray processing,and the data were stacked to form a new 3D matrix representing the granular soil’s fabric.Subsequently,utilizing the Python language and Pytorch framework,a 3D convolutional neural networks(CNNs)model was developed to establish the relationship between the constrained modulus obtained from DEM simulations and the soil’s fabric.The mean squared error(MSE)function was utilized to assess the loss value during the training process.When the learning rate(LR)fell within the range of 10-5e10-1,and the batch sizes(BSs)were 4,8,16,32,and 64,the loss value stabilized after 100 training epochs in the training and validation dataset.For BS?32 and LR?10-3,the loss reached a minimum.In the testing set,a comparative evaluation of the predicted constrained modulus from the 3D CNNs versus the simulated modulus obtained via DEM reveals a minimum mean absolute percentage error(MAPE)of 4.43%under the optimized condition,demonstrating the accuracy of this approach.Thus,by combining DEM and CNNs,the variation of soil’s mechanical characteristics related to its random fabric would be efficiently evaluated by directly tracking the particle assemblages.

Synchrotron-Based Data-Constrained Modeling Analysis of Microscopic Mineral Distributions in Limestone

Three dimensional (3D) microscopic distributions of dolomite and calcite in a limestone sample have been analyzed with a data-constrained modeling (DCM) technique using synchrotron radiation-based multi-energy X-ray computed tomography (CT) data as constraints. In order to optimize the experimental parameters, X-ray CT simulations and DCM analysis of a numerical phantom consisting of calcite (CaCO3) and dolomite (CaMg(CO3)2) have been used to investigate the effects on the predicted results in relation to noise, X-ray energy and sample-to-detector distance (SDD). The simulation results indicate that the optimal X-ray energies are 25 and 35 keVs, and the SDD is 10 mm. The high resolution 3D distributions of mineral phases of a natural limestone have been obtained. The results are useful for quantitative understanding of mineral, porosity, and physical property distributions in relation to oil and gas reservoirs hosted in carbonate rocks, which account for more than half of the world’s conventional hydrocarbon resources. The case studied is also instructive for the applicability of the DCM methods for other types of composite materials with modest atomic number contrasts between the mineral phases.

Fourth-Order Conservative Transport on Overset Grids Using Multi-Moment Constrained Finite Volume Scheme for Oceanic Modeling

With an increase in model resolution,compact high-order numerical advection scheme can improve its effectiveness and competitiveness in oceanic modeling due to its high accuracy and scalability on massive-processor computers.To provide high-quality numerical ocean simulation on overset grids,we tried a novel formulation of the fourth-order multi-moment constrained finite volume scheme to simulate continuous and discontinuous problems in the Cartesian coordinate.Utilizing some degrees of freedom over each cell and derivatives at the cell center,we obtained a two-dimensional(2D)cubic polynomial from which point values on the extended overlap can achieve fourth-order accuracy.However,this interpolation causes a lack of conservation because the flux between the regions are no longer equal;thus,a flux correction is implemented to ensure conservation.A couple of numerical experiments are presented to evaluate the numerical scheme,which confirms its approximately fourth-order accuracy in conservative transportation on overset grid.The test cases reveal that the scheme is effective to suppress numerical oscillation in discontinuous problems,which may be powerful for salinity advection computing with a sharp gradient.

Improved Differential Evolution with Shrinking Space Technique for Constrained Optimization

Most of the current evolutionary algorithms for constrained optimization algorithm are low computational efficiency. In order to improve efficiency, an improved differential evolution with shrinking space technique and adaptive trade-off model, named ATMDE, is proposed to solve constrained optimization problems. The proposed ATMDE algorithm employs an improved differential evolution as the search optimizer to generate new offspring individuals into evolutionary population. For the con- straints, the adaptive trade-off model as one of the most important constraint-handling techniques is employed to select better individuals to retain into the next population, which could effectively handle multiple constraints. Then the shrinking space technique is designed to shrink the search region according to feedback information in order to improve computational efficiency without losing accuracy. The improved DE algorithm introduces three different mutant strategies to generate different offspring into evo- lutionary population. Moreover, a new mutant strategy called ＂DE/rand/best/l＂ is constructed to generate new individuals according to the feasibility proportion ofcurrent population. Finally, the effectiveness of the pro- posed method is verified by a suite of benchmark functions and practical engineering problems. This research presents a constrained evolutionary algorithm with high efficiency and accuracy for constrained optimization problems.

Identification and nonlinear model predictive control of MIMO Hammerstein system with constraints

This work is concerned with identification and nonlinear predictive control method for MIMO Hammerstein systems with constraints. Firstly, an identification method based on steady-state responses and sub-model method is introduced to MIMO Hammerstein system. A modified version of artificial bee colony algorithm is proposed to improve the prediction ability of Hammerstein model. Next, a computationally efficient nonlinear model predictive control algorithm(MGPC) is developed to deal with constrained problem of MIMO system. The identification process and performance of MGPC are shown. Numerical results about a polymerization reactor validate the effectiveness of the proposed method and the comparisons show that MGPC has a better performance than QDMC and basic GPC.

Incorporating shape constraints in generalized additive modelling of the height-diameter relationship for Norway spruce

Background： Measurements of tree heights and diameters are essential in forest assessment and modelling. Tree heights are used for estimating timber volume, site index and other important variables related to forest growth and yield, succession and carbon budget models. However, the diameter at breast height （dbh） can be more accurately obtained and at lower cost, than total tree height. Hence, generalized height-diameter （h-d） models that predict tree height from dbh, age and other covariates are needed. For a more flexible but biologically plausible estimation of covariate effects we use shape constrained generalized additive models as an extension of existing h-d model approaches. We use causal site parameters such as index of aridity to enhance the generality and causality of the models and to enable predictions under projected changeable climatic conditions. Methods： We develop unconstrained generalized additive models （GAM） and shape constrained generalized additive models （SCAM） for investigating the possible effects of tree-specific parameters such as tree age, relative diameter at breast height, and site-specific parameters such as index of aridity and sum of daily mean temperature during vegetation period, on the h-d relationship of forests in Lower Saxony, Germany. Results： Some of the derived effects, e.g. effects of age, index of aridity and sum of daily mean temperature have significantly non-linear pattern. The need for using SCAM results from the fact that some of the model effects show partially implausible patterns especially at the boundaries of data ranges. The derived model predicts monotonically increasing levels of tree height with increasing age and temperature sum and decreasing aridity and social rank of a tree within a stand, The definition of constraints leads only to marginal or minor decline in the model statistics like AIC An observed structured spatial trend in tree height is modelled via 2-dimensional surface fitting. Conclusions： We demonstrate that the SCAM approach allows optimal regression modelling flexibility similar to the standard GAM but with the additional possibility of defining specific constraints for the model effects. The longitudinal character of the model allows for tree height imputation for the current status of forests but also for future tree height prediction.

Decentralized adaptive neural network sliding mode position/force control of constrained reconfigurable manipulators

A decentralized adaptive neural network sliding mode position/force control scheme is proposed for constrained reconfigurable manipulators. Different from the decentralized control strategy in multi-manipulator cooperation, the proposed decentralized position/force control scheme can be applied to series constrained reconfigurable manipulators. By multiplying each row of Jacobian matrix in the dynamics by contact force vector, the converted joint torque is obtained. Furthermore, using desired information of other joints instead of their actual values, the dynamics can be represented as a set of interconnected subsystems by model decomposition technique. An adaptive neural network controller is introduced to approximate the unknown dynamics of subsystem. The interconnection and the whole error term are removed by employing an adaptive sliding mode term. And then, the Lyapunov stability theory guarantees the stability of the closed-loop system. Finally, two reconfigurable manipulators with different configurations are employed to show the effectiveness of the proposed decentralized position/force control scheme.

Longitudinal height-diameter curves for Norway spruce, Scots pine and silver birch in Norway based on shape constraint additive regression models

Background： Generalized height-diameter curves based on a re-parameterized version of the Korf function for Norway spruce （Piceo abies （L.） Karst.）, Scots pine （Pinus sylvestris L.） and silver birch （Betula pendula Roth） in Norwa are presented. The Norwegian National Forest Inventory （NFI） is used as data base for estimating the model parameters. The derived models are developed to enable spatially explicit and site sensitive tree height imputatio in forest inventories as well as future tree height predictions in growth and yield scenario simulations. Methods： Generalized additive mixed models （gamm） are employed to detect and quantify potentially non-linear effects of predictor variables. In doing so the quadratic mean diameter serves as longitudinal covariate since stand ag as measured in the NFI, shows only a weak correlation with a stands developmental status in Norwegian forests. Additionally the models can be locally calibrated by predicting random effects if measured height-diameter pairs are available. Based on the model selection of non-constraint models, shape constraint additive models （scare） were fit tc incorporate expert knowledge and intrinsic relationships by enforcing certain effect patterns like monotonicity. Results： Model comparisons demonstrate that the shape constraints lead to only marginal differences in statistical characteristics but ensure reasonable model predictions. Under constant constraints the developed models predict increasing tree heights with decreasing altitude, increasing soil depth and increasing competition pressure of a tree. / two-dimensional spatially structured effect of UTM-coordinates accounts for the potential effects of large scale spatial correlated covariates, which were not at our disposal. The main result of modelling the spatially structured effect is lower tree height prediction for coastal sites and with increasing latitude. The quadratic mean diameter affects both the level and the slope of the height-diameter curve and both effects are positive. Conclusions： In this investigation it is assumed that model effects in additive modelling of height-diameter curves which are unfeasible and too wiggly from an expert point of view are a result of quantitatively or qualitatively limited data bases. However, this problem can be regarded not to be specific to our investigation but more general since growth and yield data that are balanced over the whole data range with respect to all combinations of predictor variables are exceptional cases. Hence, scare may provide methodological improvements in several applications by combining the flexibility of additive models with expert knowledge.

A Dual-mode Nonlinear Model Predictive Control with the Enlarged Terminal Constraint Sets

Aiming at a class of nonlinear systems with multiple equilibrium points, we present a dual-mode model predictive control algorithm with extended terminal constraint set combined with control invariant set and gain schedule. Local LQR control laws and the corresponding maximum control invariant sets can be designed for finite equilibrium points. It is guaranteed that control invariant sets are overlapped each other. The union of the control invariant sets is treated as the terminal constraint set of predictive control. The feasibility and stability of the novel dual-mode model predictive control are investigated with both variable and fixed horizon. Because of the introduction of extended terminal constrained set, the feasibility of optimization can be guaranteed with short prediction horizon. In this way, the size of the optimization problem is reduced so it is computationally efficient. Finally, a simulation example illustrating the algorithm is presented.

Performance Analysis of GNSS/MIMU Tight Fusion Positioning Model with Complex Scene Feature Constraints

In order to meet the requirements of high-precision vehicle positioning in complex scenes,an observation noise adaptive robust GNSS/MIMU tight fusion model based on the gain matrix is proposed considering static zero speed,non-integrity,attitude,and odometer constraint models.In this model,the robust equivalent gain matrix is constructed by the IGG-Ⅲmethod to weaken the influence of gross error,and the on-line adaptive update of observation noise matrix is carried out according to the change of actual observation environment,so as to improve the solution performance of filtering system and realize high-precision position,attitude and velocity measurement when GNSS signal is unlocked.A real test on a road over 600 km demonstrates that,in about 100 km shaded environment,the fixed rate of GNSS ambiguity resolution in the shaded road is 10%higher than that of GNSS only ambiguity resolution.For all the test,the positioning accuracy can reach the centimeter level in an open environment,better than 0.6 m in the tree shaded environment,better than 1.5 m in the three-dimensional traffic environment,and can still maintain a positioning accuracy of 0.1 m within 10 s when the satellite is unlocked in the tunnel scene.The proposal and verification of the algorithm model show that low-cost MIMU equipment can still achieve high-precision positioning when there are scene feature constraints,which can meet the problem of high-precision vehicle navigation and location in the urban complex environment.

An Adaptive Nonhydrostatic Atmospheric Dynamical Core Using a Multi-Moment Constrained Finite Volume Method

An adaptive 2 D nonhydrostatic dynamical core is proposed by using the multi-moment constrained finite-volume(MCV) scheme and the Berger-Oliger adaptive mesh refinement(AMR) algorithm. The MCV scheme takes several pointwise values within each computational cell as the predicted variables to build high-order schemes based on single-cell reconstruction. Two types of moments, such as the volume-integrated average(VIA) and point value(PV), are defined as constraint conditions to derive the updating formulations of the unknowns, and the constraint condition on VIA guarantees the rigorous conservation of the proposed model. In this study, the MCV scheme is implemented on a height-based, terrainfollowing grid with variable resolution to solve the nonhydrostatic governing equations of atmospheric dynamics. The AMR grid of Berger-Oliger consists of several groups of blocks with different resolutions, where the MCV model developed on a fixed structured mesh can be used directly. Numerical formulations are designed to implement the coarsefine interpolation and the flux correction for properly exchanging the solution information among different blocks. Widely used benchmark tests are carried out to evaluate the proposed model. The numerical experiments on uniform and AMR grids indicate that the adaptive model has promising potential for improving computational efficiency without losing accuracy.

Constraining the Generalized and Superfluid Chaplygin Gas Models with the Sandage-Loeb Test

The Sandage Loeb （SL） test is a direct measurement of the cosmic expansion by probing the redshift drifts of quasi-stellar objects in the ＇redshift desert＇ of 2 〈 z 〈 5. In this work, we investigate its constraints on the unified dark energy and dark matter models including the generalized Chaplygin gas and the superfluid Chaplygin gas. In addition, type Ia supernovae （SNIa） data and the distance ratios derived from the cosmic microwave background radiation and baryon acoustic oscillation observations （CMB/BAO） are also used. We find that the mock SL data gives the tightest constraints on the model parameters and it can help to reduce the parameter regions allowed by the present SNIa＋CMB/BAO by about 75% when all datasets considered are combined. Thus the SL test is a worthy and long awaited measurement to probe effectively the cosmic expanding history and the properties of dark energy.

Disturbance rejection tube model predictive levitation control of maglev trains

Magnetic levitation control technology plays a significant role in maglev trains.Designing a controller for the levitation system is challenging due to the strong nonlinearity,open-loop instability,and the need for fast response and security.In this paper,we propose a Disturbance-Observe-based Tube Model Predictive Levitation Control(DO-TMPLC)scheme combined with a feedback linearization strategy for the levitation system.The proposed strategy incorporates state constraints and control input constraints,i.e.,the air gap,the vertical velocity,and the current applied to the coil.A feedback linearization strategy is used to cancel the nonlinearity of the tracking error system.Then,a disturbance observer is implemented to actively compensate for disturbances while a TMPLC controller is employed to alleviate the remaining disturbances.Furthermore,we analyze the recursive feasibility and input-to-state stability of the closed-loop system.The simulation results indicate the efficacy of the proposed control strategy.

地球物理技术向油藏开发延伸的现状与展望

地球物理技术近几十年来一直不断地向油气藏开发延伸,在面向目标成像、岩石物理、多约束反演、地震约束和驱动建模、油藏动态监测与表征以及多数据融合与智能化等多方面取得了长足的进步。地球物理技术向油藏开发延伸是地球物理技术本身发展的需要,也是更好开发油气藏的需要。与勘探地球物理相比,面向油藏开发领域的地球物理技术,由于数据条件、关键技术问题与目标的变化而形成了具有自身独特特征的技术体系。旨在简单总结和梳理该领域的发展和现状,目的是能和读者分享作者对未来技术发展的思考,希望对读者有所启发。

An Enhanced Conic Reformulation for Capacity-Constrained Assortment Optimization Under the Mixture of Multinomial Logit Model

In this work,we study the conic quadratic mixed-integer formulation for assortment optimization problem under the mixture of multinomial logit(MMNL)model.The MMNL model generalizes the widely studied multinomial logit choice model and can approximate any random utility model with an arbitrary additive error.An important operational decision problem in revenue management is assortment optimization problem,which aims to find a subset of products to make available to customers that maximizes the expected revenue of the retailer.It is known that assortment optimization problem under the MMNL model is NP-hard and inapproximable within any constant performance guarantee.Commonly used methods for solving such problem are heuristical approaches or customized combinatorial optimization approaches.In the meanwhile,studies related to global optimization approaches are relatively scarce.We propose an enhanced conic quadratic mixed-integer formulation for solving assortment optimization problem under the MMNL model with a higher computational efficiency.Furthermore,we conduct extensive numerical experiments to demonstrate that the proposed reformulation significantly outperforms the existing conic reformulations for assortment optimization under the MMNL model.

AN INEXACT PROXIMAL DC ALGORITHM FOR THE LARGE-SCALE CARDINALITY CONSTRAINED MEAN-VARIANCE MODEL IN SPARSE PORTFOLIO SELECTION

Optimization problem of cardinality constrained mean-variance(CCMV)model for sparse portfolio selection is considered.To overcome the difficulties caused by cardinality constraint,an exact penalty approach is employed,then CCMV problem is transferred into a difference-of-convex-functions(DC)problem.By exploiting the DC structure of the gained problem and the superlinear convergence of semismooth Newton(ssN)method,an inexact proximal DC algorithm with sieving strategy based on a majorized ssN method(siPDCA-mssN)is proposed.For solving the inner problems of siPDCA-mssN from dual,the second-order information is wisely incorporated and an efficient mssN method is employed.The global convergence of the sequence generated by siPDCA-mssN is proved.To solve large-scale CCMV problem,a decomposed siPDCA-mssN(DsiPDCA-mssN)is introduced.To demonstrate the efficiency of proposed algorithms,siPDCA-mssN and DsiPDCA-mssN are compared with the penalty proximal alternating linearized minimization method and the CPLEX(12.9)solver by performing numerical experiments on realword market data and large-scale simulated data.The numerical results demonstrate that siPDCA-mssN and DsiPDCA-mssN outperform the other methods from computation time and optimal value.The out-of-sample experiments results display that the solutions of CCMV model are better than those of other portfolio selection models in terms of Sharp ratio and sparsity.