Gas sand distribution prediction by prestack elastic inversion based on rock physics modeling and analysis

Seismic inversion is one of the most widely used technologies for reservoir prediction. Many good results have been obtained but sometimes it fails to differentiate the lithologies and identify the fluids. However, seismic prestack elastic inversion based on rock physics modeling and analysis introduced in this paper is a significant method that can help seismic inversion and interpretation reach a new quantitative （or semi-quantitative） level from traditional qualitative interpretation. By doing rock physics modeling and forward perturbation analysis, we can quantitatively analyze the essential relationships between rock properties and seismic responses and try to find the sensitive elastic properties to the lithology, porosity, fluid type, and reservoir saturation. Finally, standard rock physics templates （RPT） can be built for specific reservoirs to guide seismic inversion interpretation results for reservoir characterization and fluids identification purpose. The gas sand distribution results of the case study in this paper proves that this method has unparalleled advantages over traditional post-stack methods, by which we can perform reservoir characterization and seismic data interpretation more quantitatively and efficiently.

Perfectly matched layer-absorbing boundary condition for finite-element time-domain modeling of elastic wave equations

The perfectly matched layer （PML） is a highly efficient absorbing boundary condition used for the numerical modeling of seismic wave equation. The article focuses on the application of this technique to finite-element time-domain numerical modeling of elastic wave equation. However, the finite-element time-domain scheme is based on the second- order wave equation in displacement formulation. Thus, the first-order PML in velocity-stress formulation cannot be directly applied to this scheme. In this article, we derive the finite- element matrix equations of second-order PML in displacement formulation, and accomplish the implementation of PML in finite-element time-domain modeling of elastic wave equation. The PML has an approximate zero reflection coefficients for bulk and surface waves in the finite-element modeling of P-SV and SH wave propagation in the 2D homogeneous elastic media. The numerical experiments using a two-layer model with irregular topography validate the efficiency of PML in the modeling of seismic wave propagation in geological models with complex structures and heterogeneous media.

Nonlinear elastic deformation of magnesium and cobalt by Preisach-Mayergoyz model

A classic hysteretic model, Preisach-Mayergoyz model （P-M model）, was used to calculate the nonlinear elastic deformation of magnesium （Mg） and cobalt （Co）. Mg and Co samples in cylinder shape were compressively tested by uniaxial test machine to obtain their stress—strain curves with hysteretic loops. The hysteretic loops do have two properties of P-M hysteretic systems： wiping out and congruency. It is proved that P-M model is applicable for the analysis of these two metals’ hysteresis. This model was applied on Mg at room temperature and Co at 300 ℃. By the P-M model, Co and Mg nonlinear elastic deformation can be calculated based on the stress history. The simulated stress—strain curves agree well with the experimental results. Therefore, the mechanical hysteresis of these two metals can be easily predicted by the classic P-M hysteretic model.

Sensitivity analysis for stochastic user equilibrium with elastic demand assignment model

This paper puts forward a rigorous approach for a sensitivity analysis of stochastic user equilibrium with the elastic demand （SUEED） model. First, proof is given for the existence of derivatives of output variables with respect to the perturbation parameters for the SUEED model. Then by taking advantage of the gradient-based method for sensitivity analysis of a general nonlinear program, detailed formulae are developed for calculating the derivatives of designed variables with respect to perturbation parameters at the equilibrium state of the SUEED model. This method is not only applicable for a sensitivity analysis of the logit-type SUEED problem, but also for the probit-type SUEED problem. The application of the proposed method in a numerical example shows that the proposed method can be used to approximate the equilibrium link flow solutions for both logit-type SUEED and probit-type SUEED problems when small perturbations are introduced in the input parameters.

On dynamics of elastic rod based on exact Cosserat model

The analysis of kinematics and dynamics of an elastic rod with circular cross section is studied on the basis of exact Cosserat model under consideration of the tension and shear deformation of the rod. The dynamical equations of a rod with arbitrary initial shape are established in general form. The dynamics of a straight rod under axial tension and torsion is discussed as an example. In discussion of static stability in the space domain the Greenhill criteria of stability and the Euler load are corrected by the influence of tension and shear strain. In analysis of dynamical stability in the time domain it is shown that the Lyapunov and Euler stability conditions of the rod in space domain are the necessary conditions of Lyapunov＇s stability in the time domain. The longitudinal, torsional and lateral vibrations of a straight rod based on exact model are discussed, and an exact formula of free frequency of lateral vibration is obtained. The free frequency formulas of various simplified models, such as the Rayleigh beam, the Kirchhoff rod, and the Timoshenko beam, can be seen as special cases of the exact formula under different conditions of simplification.

Driving Force Model for Non-pneumatic Elastic Wheel

In order to obtain good driving performance,a driving force model is presented for non-pneumatic elastic wheel.Brush model of pneumatic tyres is introduced and the deformations of elastic supports and tread are also taken into account.The longitudinal slip rate is redefined.The grounding pressure distribution of elastic wheels is analyzed and corrected according to speed,temperature and stiffness.Then rolling resistance equation is developed.Finally,simulation is conducted by software CarSim,and the results show that the estimated values are consistent with simulation values,especially at low longitudinal slip rate.The research can help to optimize design of non-pneumatic elastic wheel.

The Numerical Modeling of 3-D Elastic Wave Equation Using a High-Order,Staggered-Grid, Finite Difference Scheme

This article provides the application of the high-order, staggered-grid, finite-difference scheme to model elastic wave propagation in 3-D isotropic media. Here, we use second-order, tempo- ral- and high-order spatial finite-difference formulations with a staggered grid for discretization of the 3-D elastic wave equations of motion. The set of absorbing boundary conditions based on paraxial approximations of 3-D elastic wave equations are applied to the numerical boundaries. The trial re- sults for the salt model show that the numerical dispersion is decreased to a minimum extent, the accuracy high and diffracted waves abundant. It also shows that this method can be used for modeling wave propagation in complex media with the lateral variation of velocity.

Dynamics Modeling and Numerical Analysis of Rotor with Elastic Support/Dry Friction Dampers

The elastic support/dry friction damper is a type of damper which is used for active vibration control in a rotor system.To establish the analytical model of this type of damper,a two-dimensional friction model-ball/plate model was proposed.By using this ball/plate model,a dynamics model of rotor with elastic support/dry friction dampers was established and experimentally verified.Moreover,the damping performance of the elastic support/dry friction damper was studied numerically with respect to some variable parameters.The numerical study shows that the damping performance of the elastic support/dry friction damper is closely related to the stiffness distribution of the rotor-support system,the damper location,the pressing force between the moving and stationary disk,the friction coefficient,the tangential contact stiffness of the contact interface,and the stiffness of the stationary disk.In general,the damper should be located on an elastic support which has a large vibration amplitude in order to achieve a better damping performance,and the more vibration energy in this elastic support concentrates,the better performance of the damper will be.The larger the tangential contact stiffness of the contact interface,and the stiffness of the stationary disk are,the better performance of the damper will be.There will be an optimal value of the friction force at which the damper performs best.

A method to Estimate Dynamic Responses of VLFS Based on Multi-Floating-Module Model Connected by Elastic Hinges

Based on the elastic foundation beam theory and the multi-floating-module hydrodynamic theory,a novel method is proposed to estimate the dynamic responses of VLFS(Very Large Floating Structure).In still water,a VLFS can be simplified as an elastic foundation beam model or a multi-floating-module model connected by elastic hinges.According to equivalent displacement of the two models in static analysis,the problem of rotation stiffness of elastic hinges can be solved.Then,based on the potential flow theory,the dynamic responding analysis of multi-floatingmodule model under wave loads can be computed in ANSYS-AQWA software.By assembling the time domain analysis results of each module,the dynamic responses of the VLFS can be obtained.Validation of the method is conducted through a series of comparison calculations,which mainly includes a continuous structure and a three-part structure connected by hinges in regular waves.The results of this paper method show a satisfactory agreement with the experiment and calculation data given in relative references.

Mathematical modeling of elastic inverted pendulum control system

Inverted pendulums are important objects of theoretical investigation and experiment in the area of control theory and engineering. The researches concentrate on the rigid finite dimensional models which are described by ordinary differential equations (ODEs) .Complete rigidity is the approximation of practical models ; Elasticity should be introduced into mathematical models in the analysis of system dynamics and integration of highly precise controller. A new kind of inverted pendulum, elastic inverted pendulum was proposed, and elasticity was considered. Mathematical model was derived from Hamiltonian principle and variational methods, which were formulated by the coupling of partial differential equations (PDE) and ODE. Because of infinite dimensional, system analysis and control of elastic inverted pendulum is more sophisticated than the rigid one.

On the Elastic Vibration Model for High Length-Diameter Ratio Rocket with Attitude Control System

An elastic vibration model for high length diameter ratio spinning rocket with attitude control system which can be used for trajectory simulation is established. The basic theory of elastic dynamics and vibration dynamics were both used to set up the elastic vibration model of rocket body. In order to study the problem more conveniently, the rocket's body was simplified to be an even beam with two free ends. The model was validated by simulation results and the test data.

The Analysis of Wall Shear Stress Modulated by Acute Exercise in the Human Common Carotid Artery with an Elastic Tube Model

Assessment of the magnitude and pattern of wall shear stress(WSS)in vivo is the prerequisite for studying the quantitative relationship between exercise-induced WSS and arterial endothelial function.In the previous studies,the calculation of the WSS modulated by exercise training was primarily based upon the rigid tube model,which did not take non-linear effects of vessel elastic deformation into consideration.In this study,with an elastic tube model,we estimated the effect of a bout of 30-minute acute cycling exercise on the WSS and the flow rate in the common carotid artery according to the measured inner diameter,center-line blood flow velocity,heart rates and the brachial blood pressures before and after exercise training.Furthermore,the roles of exerciseinduced arterial diameter and blood flow rate in the change of WSS were also determined.The numerical results demonstrate that acute exercise significantly increases the magnitudes of blood flow rate and WSS.Moreover,the vessel elastic deformation is a non-negligible factor in the calculation of the WSS induced by exercise,which generates greater effects on the minimum WSS than the maximum WSS.Additionally,the contributions of exercise-induced variations in blood flow rate and diameter are almost identical in the change of the mean WSS.

Hybrid absorbing boundary condition for threedimensional elastic wave modeling

Edge reflections are inevitable in numerical modeling of seismic wavefields, and they are usually attenuated by absorbing boundary conditions. However, the commonly used perfectly matched layer （PML） boundary condition requires special treatment for the absorbing zone, and in three-dimensional （3D） modeling, it has to split each variable into three corresponding variables, which increases the computing time and memory storage. In contrast, the hybrid absorbing boundary condition （HABC） has the advantages such as ease of implementation, less computation time, and near-perfect absorption; it is thus able to enhance the computational efficiency of 3D elastic wave modeling. In this study, a HABC is developed from two-dimensional （2D） modeling into 3D modeling based on the I st Higdon one way wave equations, and a HABC is proposed that is suitable for a 3D elastic wave numerical simulation. Numerical simulation results for a homogenous model and a complex model indicate that the proposed HABC method is more effective and has better absorption than the traditional PML method.

Analytical Model of Elastic Modulus and Coefficient of Thermal Expansion for 2.5D C/SiC Composite

To make better use of 2.5D C/SiC composites in industry, it is necessary to understand the mechanical properties. A finite element model＇of 2.5D composites is established, by considering the fiber undulation and the porosity in 2.5D C/SiC composites. The fiber direction of warp is defined by cosine function to simulate the undulation of warp, and based on uniform strain assumption, analytical model of the elastic modulus and coefficient of thermal expansion （CTE） for 2.5D C/SiC composites were established by using dual- scale model. The result is found to correlate reasonably well with the predicted results and experimental results. The parametric study also demonstrates the effects of the fiber volume fraction, distance of warp yarn, and porosity in micro-scale on the mechanical properties and the coefficients of thermal expansion.

Naturally fractured hydrocarbon reservoir simulation by elastic fracture modeling

Accurate fluid flow simulation in geologically complex reservoirs is of particular importance in construction of reservoir simulators.General approaches in naturally fractured reservoir simulation involve use of unstructured grids or a structured grid coupled with locally unstructured grids and discrete fracture models.These methods suffer from drawbacks such as lack of flexibility and of ease of updating.In this study,I combined fracture modeling by elastic gridding which improves flexibility,especially in complex reservoirs.The proposed model revises conventional modeling fractures by hard rigid planes that do not change through production.This is a dubious assumption,especially in reservoirs with a high production rate in the beginning.The proposed elastic fracture modeling considers changes in fracture properties,shape and aperture through the simulation.This strategy is only reliable for naturally fractured reservoirs with high fracture permeability and less permeable matrix and parallel fractures with less cross-connections.Comparison of elastic fracture modeling results with conventional modeling showed that these assumptions will cause production pressure to enlarge fracture apertures and change fracture shapes,which consequently results in lower production compared with what was previously assumed.It is concluded that an elastic gridded model could better simulate reservoir performance.

Study of elastic wave propagation in two phase anisotropic media by numerical modeling of pseudospectral method

When there exists anisotropy in underground media elastic parameters of the observed coordinate possibly do not coincide with that of the natural coordinate. According to the theory that the density of potential energy, dissipating energy is independent of the coordinate, the relationship of elastic parameters between two coordinates is derived for two-phase anisotropic media. Then, pseudospectral method to solve wave equations of two-phase anisotropic media is derived. At last, we use this method to simulate wave propagation in two-phase anisotropic media four types of waves are observed in the snapshots, i.e., fast P wave and slow P wave, fast S wave and slow S wave. Shear wave splitting, SV wave cusps and elastic wave reflection and transmission are also observed.

Dynamic Elastic Modulus of Cement Paste at Early Age based on Nondestructive Test and Multiscale Prediction Model

This paper introduced a nondestructive testing method to evaluate the dynamic elastic modulus of cement paste. Moreover, the effect of water-cement ratio and conventional admixtures on the dynamic elastic modulus of cement paste was investigated, in which three kinds of admixtures were taken into account including viscosity modifying admixture （VMA）, silica.fume （SF）, and shrinkage-reducing admixture （SRA）. The experimental results indicate that the dynamic elastic modulus of cement paste increases with decreasing water-cement ratio. The addition of SF increases the dynamic elastic modulus, however, the overdosage of VMA causes its reduction. SRA reduces the dynamic elastic modulus at early age without affecting it in later period. Finally, a multiscale micromechanics approach coupled with a hydration model CEMHYD3D and percolation theory is utilized to predict the elastic modulus of cement paste, and the predictive results by the model are in accordance with the experimental data.

High Energy ^-pp and pp Elastic Scatterings in QCD Inspired Model

We propose QCD inspired model to calculate ^-pp and pp elastic scatterings at high energies in this paper. A calculation for total cross section of ^-pp and pp is performed in which the contributions from gluon-gluon, quark-quark, and gluon-quark interactions are included. Our results show that the QCD inspired model gives a perfect fit to experimental data of total cross section both for ^-pp and pp elastic scatterings at the whole energy region where experimental data existed at FNAL and CERN.

A Modified Monte Carlo Model of Speckle Tracking of Shear Wave Induced by Acoustic Radiation Force for Acousto-Optic Elasticity Imaging

A modified Monte Carlo model of speckle tracking of shear wave propagation in scattering media is proposed. The established Monte Carlo model mainly concerns the variations of optical electric field and speckle. The two- dimensional intensity distribution and the time evolution of speckles in different probe locations are obtained. The fluctuation of speckle intensity tracks the acoustic-radiation-force shear wave propagation, and especially the reduction of speckle intensity implies attenuation of shear wave. Then, the shear wave velocity is estimated quantitatively on the basis of the time-to-peak algorithm and linear regression processing. The results reveal that a smaller sampling interval yields higher estimation precision and the shear wave velocity is estimated more efficiently by using speckle intensity difference than by using speckle contrast difference according to the estimation error. Hence, the shear wave velocity is estimated to be 2.25 m/s with relatively high accuracy for the estimation error reaches the minimum （0.071）.

Elastic Constants of Superconducting MgB2 from Molecular Dynamics Simulations with Shell Model

The elastic constants of superconducting MgB2 are calculated using a molecular dynamics method (MD)with shell model. The lattice parameters, five independent elastic constants, equations of state (EOS), Debye temperature, and bulk modulus of MgB2 are obtained. Meanwhile, the dependence of the bulk modulus B, the lattice parameters a and c, and the unit cell volume V on the applied pressure are presented. It is demonstrated that the method introduced here can well reproduce the experimental results with a reasonable accuracy.