裂隙/孔隙型弹性介质中地震波场模拟与叠前逆时偏移

With the development of oil/gas seismic exploration, seismic survey for fracture/porosity type reservoir is becoming more and more important. As for China, since it has over 60% store of low porosity and low permeability oil/gas reservoir, it’s more urgent to validly describe fracture/porosity type oil/gas trap and proposing the related, developed seismic technique. To achieve mapping fracture/porosity region and its development status, it demands profound understanding of seismic wave propagation discipline in complex fractured/pored media. Meanwhile, it has profound scientific significance and applied worth to study forward modeling of fracture/porosity type media and pre-stacked reverse time migration. Especially, pre-stacked reverse-time migration is the lead edge technique in the field of seismology and seismic exploration. In this paper, the author has summarized the meaning, history and the present state of numerical simulation of seismic propagation in fractured/pored media and seismic exploration of fractured/pored reservoirs. Extensive Dilatancy Anisotropy (EDA) model is selected as media object in this work. As to forward modeling, due to local limitation of solving spatial partial derivative when using finite-difference and finite-element method, the author turns to pseudo-spectral method (PSM), which is based on the global characteristic of Fourier transform to simulate three-component elastic wave-field. Artifact boundary effect reduction and simulation algorithm stability are also discussed in the work. The author has completed successfully forward modeling coding of elastic wave-field and numerical simulation of two-dimensional and three-dimensional EDA models with different symmetric axis. Seismic dynamic and kinematical properties of EDA media are analyzed from time slices and seismic records of wave propagation. As to pre-stacked reverse-time migration for elastic wave-field in fractured/pored media, based on the successful experience in forward modeling results with PSM, the author has studied pre-stacked reverse-time depth-domain migration technique using PSM of elastic wave-field in two dimensional EDA media induced by preferred fracture/pore distribution. At the same time, different image conditions will bring up what kind of migration result is detailed in this paper. The author has worded out software for pre-stacked reverse-time depth-domain migration of elastic wave-field in EDA media. After migration processing of a series of seismic shot gathers, influences to migration from different isotropic and anisotropy models are described in the paper. In summary, following creative research achievements are obtained:  Realizing two-dimensional and three-dimensional elastic wave-field modeling for fractured/pored media and related software has been completed.  Proposed pre-stacked reverse-time depth-domain migration technique using PSM of elastic wave-field.  Through analysis of the seismic dynamic and kinematical properties of EDA media, the author made a conclusion that collection of multi-component seismic data can provide important data basis for locating and describing the fracture/pore regions and their magnitudes and the preferred directions.  Pre-stacked reverse-time depth-domain migration technique has the ability to reconstruct complex geological object with steep formations and tilt fracture distribution. Neglecting seismic anisotropy induced by the preferred fracture/pore distribution, will lead to the disastrous imaging results.

地震波模拟：非规则网格离散的数值解法、吸收边界及其实际应用

This dissertation presents a series of irregular-grid based numerical technique for modeling seismic wave propagation in heterogeneous media. The study involves the generation of the irregular numerical mesh corresponding to the irregular grid scheme, the discretized version of motion equations under the unstructured mesh, and irregular-grid absorbing boundary conditions. The resulting numerical technique has been used in generating the synthetic data sets on the realistic complex geologic models that can examine the migration schemes. The motion equation discretization and modeling are based on Grid Method. The key idea is to use the integral equilibrium principle to replace the operator at each grid in Finite Difference scheme and variational formulation in Finite Element Method. The irregular grids of complex geologic model is generated by the Paving Method, which allow varying grid spacing according to meshing constraints. The grids have great quality at domain boundaries and contain equal quantities of nodes at interfaces, which avoids the interpolation of parameters and variables. The irregular grid absorbing boundary conditions is developed by extending the Perfectly Matched Layer method to the rotated local coordinates. The splitted PML equations of the first-order system is derived by using integral equilibrium principle. The proposed scheme can build PML boundary of arbitrary geometry in the computational domain, avoiding the special treatment at corners in a standard PML method and saving considerable memory and computation cost. The numerical implementation demonstrates the desired qualities of irregular grid based modeling technique. In particular, (1) smaller memory requirements and computational time are needed by changing the grid spacing according to local velocity; (2) Arbitrary surfaces and interface topographies are described accurately, thus removing the artificial reflection resulting from the stair approximation of the curved or dipping interfaces; (3) computational domain is significantly reduced by flexibly building the curved artificial boundaries using the irregular-grid absorbing boundary conditions. The proposed irregular grid approach is apply to reverse time migration as the extrapolation algorithm. It can discretize the smoothed velocity model by irregular grid of variable scale, which contributes to reduce the computation cost. The topography. It can also handle data set of arbitrary topography and no field correction is needed.

二氧化碳在孔隙咸水层中多相流动数值模拟研究

There has been a growing concern about the use of fossil fuels and its adverse effects on the atmospheric greenhouse and ecological environment. A reduction in the release rate of CO2 into the atmosphere poses a major challenge to the land ecology of China. The most promising way of achieving CO2 reduction is to dispose of CO2 in deep saline aquifers. Deep aquifers have a large potential for CO2 sequestration in geological medium in terms of volume and duration. Through the numerical simulation of multiphase flow in a porous media, the transformation and motion of CO2 in saline aquifers has been implemented under various temperature and hydrostatic pressure conditions, which plays an important role to the assessment of the reliability and safety of CO2 geological storage. As expected, the calculated results can provide meaningful and scientific information for management purposes. The key problem to the numerical simulation of multiphase flow in a porous media is to accurately capture the mass interface and to deal with the geological heterogeneity. In this study, the updated CE/SE (Space and time conservation element and solution element) method has been proposed, and the Hybrid Particle Level Set method (HPLS) has extended for multiphase flows in porous medium, which can accurately trace the transformation of the mass interface. The benchmark problems have been applied to evaluate and validate the proposed method. In this study, the reliability of CO2 storage in saline aquifers in Daqingzi oil field in Sunlong basin has been discussed. The simulation code developed in this study takes into account the state for CO2 covering the triple point temperature and pressure to the supercritical region. The geological heterogeneity has been implemented, using the well known geostatistical model (GSLIB) on the base of the hard data. The 2D and 3D model have been set up to simulate the CO2 multiphase flow in the porous saline aquifer, applying the CE/SE method and the HPLS method .The main contents and results are summarized as followings. (1) The 2D CE/SE method with first and second –order accuracy has been extended to simulate the multiphase flow in porous medium, which takes into account the contribution of source and sink in the momentum equation. The 3D CE/SE method with the first accuracy has been deduced. The accuracy and efficiency of the proposed CE/SE method have been investigated, using the benchmark problems. (2) The hybrid particle level set method has been made appropriate and extended for capturing the mass interface of multiphase flows in porous media, and the numerical method for level set function calculated has been formulated. (3) The closed equations for multiphase flow in porous medium has been developed, adept to both the Darcy flow and non-Darcy flow, getting over the limitation of Reynolds number to the calculation. It is found that Darcy number has a decisive influence on pressure as well as velocity given the Darcy number. (4) The new Euler scheme for numerical simulations of multiphase flows in porous medium has been proposed, which is efficient and can accurately capture the mass interface. The artificial compressibility method has been used to couple the velocities and pressure. It is found that the Darcy number has determinant effects on the numerical convergence and stability. In terms of the different Darcy numbers, the coefficient of artificial compressibility and the time step have been obtained. (5) The time scale of the critical instability for critical CO2 in the saline aquifer has been found, which is comparable with that of completely CO2 dissolved saline aquifer. (6) The concept model for CO2 multiphase flows in the saline aquifer has been configured, based on the temperature, pressure, porosity as well as permeability of the field site .Numerical simulation of CO2 hydrodynamic trapping in saline aquifers has been performed, applying the proposed CE/SE method. The state for CO2 has been employed to take into account realistic reservoir conditions for CO2 geological sequestration. The geological heterogeneity has been sufficiently treated , using the geostatistical model. (7) It is found that the Rayleigh-Taylor instability phenomenon, which is associated with the penetration of saline fluid into CO2 fluid in the direction of gravity, has been observed in CO2 multiphase flows in the saline aquifer. Development of a mushroom-type spike is a strong indication of the formation of Kelvin-Helmholtz instability due to the developed short wavelength perturbations present along the interface and parallel to the bulk flow. Additional key findings: the geological heterogeneity can distort the flow convection. The ascending of CO2 can induce the persistent flow cycling effects. The results show that boundary conditions of the field site have determinant effects on the transformation and motion of CO2 in saline aquifers. It is confirmed that the proposed method and numerical model has the reliability to simulate the process of the hydrodynamic trapping, which is the controlling mechanism for the initial period of CO2 storage at time scale of 100 years.

地幔对流理论及其在解释地球物理观测的应用

With the improving of mantle convection theory, the developing of computing method and increasing of the measurement data, we can numerically simulate more clearly about the effects on some geophysical observed phenomenons such as the global heat flow and global lithospheric stress field in the Earth's surface caused by mantle convection, which is the primary mechanism for the transport of heat from the Earth's deep interior to its surface and the underlying force mechanism of dynamics in the Earth.Chapter 1 reviews the historical background and present research state of mantle convection theory.In Chapter 2, the basic conception of thermal convection and the basic theory about mantle flow.The effects on generation and distribution of global lithospheric stres s field induced by mantle flow are the subject of Chapter 3. Mantle convection causes normal stress and tangential stresses at the bottom of the lithosphere, and then the sublithospheric stress field induces the lithospheric deformation as sixrface force and results in the stress field within the lithosphere. The simulation shows that the agreement between predictions and observations is good in most regions. Most of subduction zones and continental collisions are under compressive. While ocean ridges, such as the east Pacific ridge, the Atlantic ridge and the east African rift valley, are under tensile. And most of the hotspots preferentially occur in regions where calculated stress is tensile. The calculated directions of the most compressive principal horizontal stress are largely in accord with that of the observation except for some regions such as the NW-Pacifie subduction zone and Qinghai-Tibet Plateau, in which the directions of the most compressive principal horizontal stress are different. It shows that the mantel flow plays an important role in causing or affecting the large-scale stress field within the lithosphere.The global heat flow simulation based on a kinematic model of mantle convection is given in Chapter 4. Mantle convection velocities are calculated based on the internal loading theory at first, the velocity field is used as the input to solve the thermal problem. Results show that calculated depth derivatives of the near surface temperature are closely correlated to the observed surface heat flow pattern. Higher heat flow values around midocean ridge systems can be reproduced very well. The predicted average temperature as a function of function of depth reveals that there are two thermal boundary layers, one is close to the surface and another is close to the core-mantle boundary, the rest of the mantle is nearly isothermal. Although, in most of the mantle, advection dominates the heat transfer, the conductive heat transfer is still locally important in the boundary layers and plays an important role for the surface heat flow pattern. The existence of surface plates is responsible for the long wavelength surface heat flow pattern.In Chapter 5, the effects on present-day crustal movement in the China Mainland resulted from the mantle convection are introduced. Using a dynamic method, we present a quantitative model for the present-day crustal movement in China. We consider not only the effect of the India-Eurasia collision, the gravitational potential energy difference of the Tibet Plateau, but also the contribution of the shear traction on the bottom of the lithosphere induced by the global mantle convection. The comparison between our results and the velocity field obtained from the GPS observation shows that our model satisfactorily reproduces the general picture of crustal deformation in China. Numerical modeling results reveal that the stress field on the base of the lithosphere induced by the mantle flow is probably a considerable factor that causes the movement and deformation of the lithosphere in continental China with its eflfcet focuing on the Eastern China A numerical research on the small-scale convection with variable viscosity in the upper mantle is introduced in Chapter 6. Based on a two-dimensional model, small-scale convection in the mantle-lithosphere system with variable viscosity is researched by using of finite element method. Variation of viscosity in exponential form with temperature is considered in this paper The results show that if viscosity is strongly temperature-dependent, the upper part of the system does not take a share in the convection and a stagnant lid, which is identified as lithosphere, is formed on the top of system because of low temperature and high viscosity. The calculated surface heat flow, topography and gravity anomaly are associated well with the convection pattern, namely, the regions with high heat flow and uplift correspond to the upwelling flow, and vice versa.In Chapter 7, we give a brief of future research subject: The inversion of lateral density heterogeneity in the mantle by minimizing the viscous dissipation.

A multi-scale atomistic-continuum modelling of crack propagation in a two-dimensional macroscopic plate

A novel multi-scale seamless model of brittle-crack propagation is proposed and applied to the simulation of fracture growth in a two-dimensional Ag plate with macroscopic dimensions. The model represents the crack propagation at the macroscopic scale as the drift-diffusion motion of the crack tip alone. The diffusive motion is associated with the crack-tip coordinates in the position space, and reflects the oscillations observed in the crack velocity following its critical value. The model couples the crack dynamics at the macroscales and nanoscales via an intermediate mesoscale continuum. The finite-element method is employed to make the transition from the macroscale to the nanoscale by computing the continuum-based displacements of the atoms at the boundary of an atomic lattice embedded within the plate and surrounding the tip. Molecular dynamics (MD) simulation then drives the crack tip forward, producing the tip critical velocity and its diffusion constant. These are then used in the Ito stochastic calculus to make the reverse transition from the nanoscale back to the macroscale. The MD-level modelling is based on the use of a many-body potential. The model successfully reproduces the crack-velocity oscillations, roughening transitions of the crack surfaces, as well as the macroscopic crack trajectory. The implications for a 3-D modelling are discussed.

Finite-Element Time-Step Coupled Generator, Load, AVR and Brushless Exciter Modelling

This paper presented results from a details and comprehensive simulation using finite element method of the practical operation of an electrical machine. The results it displayed have been used in practice to design more efficient equipment.

Quality factor assessment of subwavelength resonant cavities at FIR frequencies

Subwavelength resonators at FIR are presented and studied. The structures consist of 1D cavities formed between a metallized (silver) surface and a metamaterial surface comprising a periodic array of silver patches on a silver-backed silicon substrate. The concept derives from recent discoveries of artificial magnetic conductors (AMC). By studying the currents excited on the metamaterial surface by a normally incident plane wave, the nature of the emerging resonant phenomena and the physical mechanism underlying the AMC operation are investigated. Full wave simulations, based on finite element method and time-domain transmission line modelling technique, have been carried out to demonstrate the effective AMC boundary condition and prove the possibilities for subwavelength cavities. The quality factor of the resonant cavities is assessed as a function of the cavity profile. It is demonstrated that the quality factor drops to about 1/8 of the half-wavelength value for lambda/8 resonant cavity.

Analysis of compressive membrane action in concrete slabs

This paper summarises the results obtained from non-linear finite-element analysis (NLFEA) of a series of reinforced-concrete one-way slabs with various boundary conditions representative of a bridge deck slab strip in which compressive membrane action governs the structural behaviour. The application of NLFEA for the optimum analysis and design of in-plane restrained concrete slabs is explored. An accurate material model and various equation solution methods were assessed to find a suitable finite-element method for the analysis of concrete slabs in which arching action occurs. Finally, the results from the NLFEA are compared and validated with those from various experimental test data. Significantly, the numerical analysis was able to model the arching action that occurred as a result of external in-plane restraint at the supports and which enhanced the ultimate strength of the slab. The NLFEA gave excellent predictions for the ultimate load-carrying capacity and far more accurate predictions than those obtained using standard flexural or elastic theory.

Finite Element Application to Analysis, Testing, and Manufacturing of Satellite Structures

Designing satellite structures poses an ongoing challenge as the interaction between analysis, experimental testing, and manufacturing phases is underdeveloped. Finite Element Analysis for Satellite Structures: Applications to Their Design, Manufacture and Testing explains the theoretical and practical knowledge needed to perform design of satellite structures. By layering detailed practical discussions with fully developed examples, Finite Element Analysis for Satellite Structures: Applications to Their Design, Manufacture and Testing provides the missing link between theory and implementation.
Computational examples cover all the major aspects of advanced analysis; including modal analysis, harmonic analysis, mechanical and thermal fatigue analysis using finite element method. Test cases are included to support explanations an a range of different manufacturing simulation techniques are described from riveting to shot peening to material cutting. Mechanical design of a satellites structures are covered in three steps: analysis step under design loads, experimental testing to verify design, and manufacturing.
Stress engineers, lecturers, researchers and students will find Finite Element Analysis for Satellite Structures: Applications to Their Design, Manufacture and Testing a key guide on with practical instruction on applying manufacturing simulations to improve their design and reduce project cost, how to prepare static and dynamic test specifications, and how to use finite element method to investigate in more details any component that may fail during testing.

Reliability Assessment of 3-D Space Frame Structures Applying Stochastic Finite Element Analysis

Throughout design development of satellite structure, stress engineer is usually challenged with randomness in applied loads and material properties. To overcome such problem, a risk-based design is applied which estimates satellite structure probability of failure under static and thermal loads. Determining probability of failure can help to update initially applied factors of safety that were used during structure preliminary design phase. These factors of safety are related to the satellite mission objective. Sensitivity-based analysis is to be implemented in the context of finite element analysis (probabilistic finite element method or stochastic finite element method (SFEM)) to determine the probability of failure for satellite structure or one of its components.

Transient Heat Conduction of Variable Stiffness Composite Laminate

Recent research on Variable Stiffness (VS) laminates, which are constructed by steering the fiber orientation as a spatial function of location, have shown to improve laminate performance under mechanical loads. Two distinct cases of stiffness variation can be achieved either by variation of the fiber orientation in the direction of the global x-axis, or perpendicular to it. In the present paper, thermal analysis of a VS laminate is performed to study the effect of steering fibers on transient heat conduction under uniform heat flux using finite element method. The goal of the present paper is a parametric study of the effect of variable stiffness properties on transient response including time to reach steady state and temperature profile. Also, stress resultants and maximum stress location are investigated under different boundary conditions. A FEM algorithm is applied to exactly incorporate the boundary conditions for stress resultant analysis.

Transient Heat Conduction of Variable Stiffness Composite Laminate

Recent research on Variable Stiffness (VS) laminates, which are constructed by steering the fiber orientation as a spatial function of location, have shown to improve laminate performance under mechanical loads. Two distinct cases of stiffness variation can be achieved either by variation of the fiber orientation in the direction of the global x-axis, or perpendicular to it. In the present paper, thermal analysis of VS laminate is performed to study the effect of steering fibers on transient heat conduction under uniform heat flux using finite element method. The goal of the present paper is a parametric study of the
effect of variable stiffness properties on transient response including time to reach steady state and temperature profile. Also, stress resultants and maximum stress location are investigated under different boundary conditions. A FEM algorithm is applied to exactly incorporate the boundary conditions.

The development of a surface defect machining method for hard turning processes

Hard turning (HT) is a material removal process employing a combination of a single point cutting tool and high speeds to machine hard ferrous alloys which exhibit hardness values over 45 HRC. In this paper, a surface defect machining (SDM) method for HT is proposed which harnesses the combined advantages of porosity machining and pulsed laser pre-treatment processing. From previous experimental work, this was shown to provide better controllability of the process and improved quality of the machined surface. While the experiments showed promising results, a comprehensive understanding of this new technique could only be achieved through a rigorous, in depth theoretical analysis. Therefore, an assessment of the SDM technique was carried out using both finite element method (FEM) and molecular dynamics (MD) simulations.
FEM modelling was used to compare the conventional HT of AISI 4340 steel (52 HRC) using an Al2O3 insert with the proposed SDM method. The simulations showed very good agreement with the previously published experimental results. Compared to conventional HT, SDM provided favourable machining outcomes, such as reduced shear plane angle, reduced average cutting forces, improved surface roughness, lower residual stresses on the machined surface, reduced tool–chip interface contact length and increased chip flow velocity. Furthermore, a scientific explanation of the improved surface finish was revealed using a state-of-the-art MD simulation model which suggested that during SDM, a combination of both the cutting action and rough polishing action help improve the machined surface finish.

Finite element modeling of debonding failures in FRP-strengthened RC beams: A dynamic approach

This paper is concerned with the finite element simulation of debonding failures in FRP-strengthened concrete beams. A key challenge for such simulations is that common solution techniques such as the Newton-Raphson method and the arc-length method often fail to converge. This paper examines the effectiveness of using a dynamic analysis approach in such FE simulations, in which debonding failure is treated as a dynamic problem and solved using an appropriate time integration method. Numerical results are presented to show that an appropriate dynamic approach effectively overcomes the convergence problem and provides accurate predictions of test results.

Experimental and discrete element modelling of cohesive iron ore fines

In this study, the behaviour of iron ore fines with varying levels of adhesion was investigated using a confined compression test and a uniaxial test. The uniaxial test was conducted using the semi-automated uniaxial EPT tester in which the cohesive strength of a bulk solid is evaluated from an unconfined compression test following a period of consolidation to a pre-defined vertical stress. The iron ore fines were also tested by measuring both the vertical and circumferential strains on the cylindrical container walls under vertical loading in a separate confined compression tester - the K0 tester, to determine the lateral pressure ratio. Discrete Element Method simulations of both experiments were carried out and the predictions were compared with the experimental observations. A recently developed DEM contact model for cohesive solids, an Elasto-Plastic Adhesive model, was used. This particle contact model uses hysteretic non-linear loading and unloading paths and an adhesion parameter which is a function of the maximum contact overlap. The model parameters for the simulations are phenomenologically based to reproduce the key bulk characteristics exhibited by the solid. The simulation results show a good agreement in capturing the stress history dependent behaviour depicted by the flow function of the cohesive iron ore fines while also providing a reasonably good match for the lateral pressure ratio observed during the confined compression K0 tests. This demonstrates the potential for the DEM model to be used in the simulation of bulk handling applications.