Dynamic Field Division of Hydrocarbon Migration,Accumulation and Hydrocarbon Enrichment Rules in Sedimentary Basins

Hydrocarbon distribution rules in the deep and shallow parts of sedimentary basins are considerably different, particularly in the following four aspects. First, the critical porosity for hydrocarbon migration is much lower in the deep parts of basins： at a depth of 7000 m, hydrocarbons can accumulate only in rocks with porosity less than 5%. However, in the shallow parts of basins （i.e., depths of around 1000 m）, hydrocarbon can accumulate in rocks only when porosity is over 20%. Second, hydrocarbon reservoirs tend to exhibit negative pressures after hydrocarbon accumulation at depth, with a pressure coefficient less than 0.7. However, hydrocarbon reservoirs at shallow depths tend to exhibit high pressure after hydrocarbon accumulation. Third, deep reservoirs tend to exhibit characteristics of oil （-gas）-water inversion, indicating that the oil （gas） accumulated under the water. However, the oil （gas） tends to accumulate over water in shallow reservoirs. Fourth, continuous unconventional tight hydrocarbon reservoirs are distributed widely in deep reservoirs, where the buoyancy force is not the primary dynamic force and the caprock is not involved during the process of hydrocarbon accumulation. Conversely, the majority of hydrocarbons in shallow regions accumulate in traps with complex structures. The results of this study indicate that two dynamic boundary conditions are primarily responsible for the above phenomena： a lower limit to the buoyancy force and the lower limit of hydrocarbon accumulation overall, corresponding to about 10%-12% porosity and irreducible water saturation of 100%, respectively. These two dynamic boundary conditions were used to divide sedimentary basins into three different dynamic fields of hydrocarbon accumulation： the free fluid dynamic field, limit fluid dynamic field, and restrain fluid dynamic field. The free fluid dynamic field is located between the surface and the lower limit of the buoyancy force, such that hydrocarbons in this field migrate and accumulate under the influence of, for example, the buoyancy force, pressure, hydrodynamic force, and capillary force. The hydrocarbon reservoirs formed are characterized as ＂four high,＂ indicating that they accumulate in high structures, are sealed in high locations, migrate into areas of high porosity, and are stored in reservoirs at high pressure. The basic features of distribution and accumulation in this case include hydrocarbon migration as a result of the buoyancy force and formation of a reservoir by a caprock. The limit fluid dynamic field is located between the lower limit of the buoyancy force and the lower limit of hydrocarbon accumulation overall; the hydrocarbon migrates and accumulates as a result of, for example, the molecular expansion force and the capillary force. The hydrocarbon reservoirs formed are characterized as ＂four low,＂ indicating that hydrocarbons accumulate in low structures, migrate into areas of low porosity, and accumulate in reservoirs with low pressure, and that oil（-gas）-water inversion occurs at low locations. Continuous hydrocarbon accumulation over a large area is a basic feature of this field. The restrain fluid dynamic field is located under the bottom of hydrocarbon accumulation, such that the entire pore space is filled with water. Hydrocarbons migrate as a result of the molecular diffusion force only. This field lacks many of the basic conditions required for formation of hydrocarbon reservoirs： there is no effective porosity, movable fluid, or hydrocarbon accumulation, and potential for hydrocarbon exploration is low. Many conventional hydrocarbon resources have been discovered and exploited in the free fluid dynamic field of shallow reservoirs, where exploration potential was previously considered to be low. Continuous unconventional tight hydrocarbon resources have been discovered in the limit fluid dynamic field of deep reservoirs; the exploration potential of this setting is thought to be tremendous, indicating that future exploration should be focused primarily in this direction.

Fluid Dynamic Field in BozhongDepression, Bohai Bay Basin

The data from regional geology, boreholes, geophysics and tests are integrated to analyze the fluid dynamic field in the Bozhong depression, Bohai Bay basin. The current geothermal gradient is determined to be about 2.95 /100 m by integrating 266 drill-stem test (DST) measurements and comparing with the global average value. The paleogeothermal gradients are calculated from the homogenization temperatures of saline inclusions, which vary both laterally and vertically. The data from sonic logs, well tests and seismic velocities are used to investigate the pressure variations in the study area. The mudstone compaction is classified as three major types: normal compaction and normal pressure, under-compaction and overpressure, and past-compaction and under-overpressure. The current pressure profile is characterized by normal pressure, sight pressure and intense overpressure from top to bottom The faults, unconformity surfaces and interconnecting pores constitute a complex network of vertical and horizontal fluid flows within the depression. The fluid potential energy profiles present a 'double-deck' structure. The depocenters are the area of fluids supply, whereas the slopes and uplifts are the main areas of fluids charge.

Fast impurity solver for dynamical mean field theory based on second order perturbation around the atomic limit

This paper proposes an impurity solver for the dynamical mean field theory （DMFT） study of the Mott insulators, which is based on the second order perturbation of the hybridization function. After careful benchmarking with quantum Monte Carlo results on the anti-ferromagnetic phase of the Hubbard model, it concludes that this impurity solver can capture the main physical features in the strong coupling regime and can be a very useful tool for the LDA （local density approximation） ＋ DMFT studies of the Mort insulators with long range order.

MEASUREMENT OF WELDING DYNAMIC DISPLACEMENT FIELDS BY ELECTRONIC SPECKLE PATTERN INTERFERENCE

In order to effectively control the stress and distortion which produced in welding process, the dynamic change laws of displacement field is the most important factor. The characteristics of the welding dynamic displacement field is high temperature, high strain velocity, thus ordinary methods such as resistance strain gauge or Moiré method can not be used for the measurement of the zone of high temperature. Speckle interference method has the merits of non-contact, resistance to the disturbance of impure lights, high accuracy of measurement (half of wavelength).The paper represents the measurement of dynamic displacement field of argon-arcspot welding, by which it shows that the method of speckle interference is feasible for the measurement of welding dynamic displacement.

Flow Field and Temperature Field in GaN-MOCVD Reactor Based on Computational Fluid Dynamics Modeling

Metal organic chenlical vapor deposition （AIOCVD） growth systems arc one of the. main types of equipment used for growing single crystal materials, such as GaN. To obtain fihn epitaxial materials with uniform performanee, the flow field and ternperature field in a GaN-MOCVD reactor are investigated by modeling and simulating. To make the simulation results more consistent with the actual situation, the gases in the reactor are considered to be compressible, making it possible to investigate the distributions of gas density and pressure in the reactor. The computational fluid dynamics method is used to stud,v the effects of inlet gas flow velocity, pressure in the reactor, rotational speed of graphite susceptor, and gases used in the growth, which has great guiding~ significance for the growth of GaN fihn materials.

Double Ionization Dynamics of Molecular Hydrogen in Ultrashort Intense Laser Fields

We experimentally investigate the double ionization pulses. The total kinetic energy release of the two of molecular hydrogen subjected to ultrashort intense laser coincident H＋ ions, which provides a diagnosis of different processes to double ionization of H2, is measured for two different pulse durations, i.e., 25 and 5 fs, and various laser intensities. It is found that, for the long pulse duration （i.e., 25 fs）, the double ionization occurs mainly via two processes, i.e., the charge resonance enhanced ionization and recollision-induced double ionization. Moreover, the contributions from these two processes can be significantly modulated by changing the laser intensity. In contrast, for a few-cycle pulse of 5 fs, only the recollsion-induced double ionization survives, and in particular, this process could be solely induced by the first-return reeollision at appropriate laser intensities, providing an efficient way to probe the sub-laser-cycle molecular dynamics.

Saturated sodium chloride solution under an external static electric field: A molecular dynamics study

The behavior of saturated aqueous Na Cl solutions under a constant external electric field（E） was studied by molecular dynamics（MD） simulation. Our dynamic MD simulations indicated that the irreversible nucleation process towards crystallization is accelerated by a moderate E but retarded or even prohibited under a stronger E, which can be understood by the competition between self-diffusion and drift motion. The former increases with E, thereby accelerating the nucleation process, whereas the latter pulls oppositely charged ions apart under a stronger E, thereby decelerating nucleation.Additionally, our steady-state MD simulations indicated that a first-order phase transition occurs in saturated solutions at a certain threshold Ec. The magnitude of Ec increases with concentration because larger clusters form more easily when the solution is more concentrated and require a stronger E to dissociate.

Field structure at mode Ⅲ dynamically propagating crack tip in elastic-viscoplastic materials

An elastic-viscoplastic mechanics model is used to investigate asymptotically the mode Ⅲ dynamically propagating crack tip field in elastic-viscoplastic materials. The stress and strain fields at the crack tip possess the same power-law singularity under a linear-hardening condition. The singularity exponent is uniquely determined by the viscosity coefficient of the material. Numerical results indicate that the motion parameter of the crack propagating speed has little effect on the zone structure at the crack tip. The hardening coefficient dominates the structure of the crack-tip field. However, the secondary plastic zone has little influence on the field. The viscosity of the material dominates the strength of stress and strain fields at the crack tip while it does have certain influence on the crack-tip field structure. The dynamic crack-tip field degenerates into the relevant quasi-static solution when the crack moving speed is zero. The corresponding perfectly-plastic solution is recovered from the linear-hardening solution when the hardening coefficient becomes zero.

Flow Field Characteristics for Parachute-projectile System

In the current research for parachute flow field nowadays,the size of parachutes in previous research are so large compared with their carriers that the effects of the carriers wake flow to parachute are always neglected.Different from such large parachutes,the parachute size in this paper is on the same magnitude with the carrier,thus,the carrier can obviously affect the parachute flow field.In this paper,flow field characteristics of small parachute for projectile decelerating are researched through two approaches,namely,computational fluid dynamics（CFD） simulation and wind tunnel tests.Three parachutes with various sizes are chosen for study.Firstly,the CFD simulation of flow field around these parachutes is carried out,and then the CFD simulation of parachute-projectile systems is executed.According to the simulation results,the phenomenon is observed that in the simulations of parachutes there are two vortex-rings at the wind shadow of parachutes,however,in the second simulations of parachute-projectile systems,two additional vortex-rings emerge inside the parachutes.Due to these two inner vortex-rings,the pressure inside parachutes decreases.As a result,the drag of parachute in simulation of parachute-projectile systems is about 20% smaller compared with the prior one.In order to verify the numerical results of CFD simulations,wind tunnel tests are employed.In terms of the data of the wind tunnel tests,the CFD simulation for flow field characteristics is reasonable and feasible.The results of both CFD simulation and wind tunnel tests demonstrated the influence of projectile wake flow to parachute drag can not be neglected if the parachute size is on the same magnitude with projectile.The influence to parachute drag from the ratio of projectile diameter to parachute diameter is also analyzed both in CFD simulations and wind tunnel tests.The approach combined CFD simulation and wind tunnel tests proposed can be used to guide the design of such parachute whose size is on the same magnitude with carrier.

ReaxFF molecular dynamics study on oxidation behavior of 3C-SiC:Polar face effects

The oxidation of nanoscale 3C-SiC involving four polar faces（C（100）, Si（100）, C（111）, and Si（111）） is studied by means of a reactive force field molecular dynamics（Reax FF MD） simulation. It is shown that the consistency of 3C-SiC structure is broken over 2000 K and the low-density carbon chains are formed within SiC slab. By analyzing the oxygen concentration and fitting to rate theory, activation barriers for C（100）, Si（100）, C（111）, and Si（111） are found to be 30.1,35.6, 29.9, and 33.4 k J·mol^-1. These results reflect lower oxidative stability of C-terminated face, especially along [111] direction. Compared with hexagonal polytypes of SiC, cubic phase may be more energy-favorable to be oxidized under high temperature, indicating polytype effect on SiC oxidation behavior.

Evolution of windblown sand flux and dune field——Trans-scale modeling and simulation

Windblown sand flux and dune field evolving toward the oasis have been a common ecological and environmental threat confronted by many countries.Meanwhile,it is also a kind of complex dynamical process involving multiple temporal and spatial scales which is still out of accurate description through current field observations.Available models and reliable quantitative simulations are of significant value to predict the spreading rate of desertification and provide an optimal design for sand prevention.This paper presents a 'triple-jump' method to realize quantitative simulations to the formation and evolution of an aeolian dune field from an arbitrary initial configuration.Simulated results achieve a satisfactory agreement with observations qualitatively and quantitatively,which also reveal the characteristics and dynamical behaviors of dunes and dune field.Such a paradigm is of a good level of generality,which provides an exploratory probe into the subject of multi-scale physics.

Application of FLUENT on fine-scale simulation of wind field over complex terrain

The state-of-art Computational Fluid Dynamics （CFD） codes FLUENT is applied in a fine-scale simulation of the wind field over a complex terrain. Several numerical tests are performed to validate the capability of FLUENT on describing the wind field details over a complex terrain. The results of the numerical tests show that FLUENT can simulate the wind field over extremely complex terrain, which cannot be simulated by mesoscale models. The reason why FLUENT can cope with extremely complex terrain, which can not be coped with by mesoscale models, relies on some particular techniques adopted by FLUENT, such as computer-aided design （CAD） technique, unstructured grid technique and finite volume method. Compared with mesoscale models, FLUENT can describe terrain in much more accurate details and can provide wind simulation results with higher resolution and more accuracy.

An experimental study of the dynamic features of shadow areas of caustics in response to loading/unloading fracture

Using a plexiglass sample and by means of real-time holographic interferometry and shadow optical method of caustics, the different features of dynamic variation in stress (strain) field, plastic area and nucleation zone (shadow area) when the sample fractures during loading (loading-fracture) and unloading (unloading-fracture) are studied visually. The results show that the strain nuclei (zones with dense fringes) appear first at the tips of prefabricated cracks at low stress, and then the shadow areas of caustics form with the increase of load. These nuclei and shadow areas can become larger, or smaller, when the process of loading, or unloading, goes on. When the stress is kept within a certain range, the shadow areas of caustics can become larger and smaller alternatively with repeated loading and unloading (cyclic loading). However, when loading and unloading at high stress, in particular when the macrofracture is about to appear, the variations of the shadow areas of caustics are irreversible and quite different. The shadow areas of caustics expand rapidly at an increasing speed when loading-fracture appears. In contrast, the shadow areas of caustics expand at a lower speed when unloading-fracture appears; besides, there is a circular shadow in front of the sharp-angle shaped area.

Earth gravity field solution with combining CHAMP and GRACE data

Satellite gravity data fusion with multi-type and huge-amount is one of the hot topics in physical geodesy. After a brief review of dynamic approach, the CHAMP-only and GRACE-only gravity fields by using HL-SST and LL-SST data from 2003 to 2009 are recovered respectively. An combination strategy of CHAMP and GRACE data by using Helmert variance component estimation （VCE） is proposed based on normal equation level fusion. Three gravity field models with 150° and order by CHAMP-only data, GRACE-only data and combining CHAMP and GRACE data from 2003 to 2009 are recovered. The comparisons between our recovered models and those latest released models were performed. The external accuracy validations using marine gravity anomalies from DTU13 products and height anomalies from GPS/leveling data are also conducted in this paper. The results show that long-term CHAMP data do contribute to the accuracy improvement of gravity field solution. The accuracy of the combined model using CHAMP and GRACE data is better than those of the individuals and comparative to the models published by international groups.

Genetic algorithm-finite element method inversion of the factors determining the recent tectonic stress field of part of East Asia area

Genetic algorithm finite element method (GA FEM) is applied to the study of tectonic stress field of part of East Asia area. From the observed stress distribution, 2 D elastic plane stress inversion is made to deduce the boundary forces and investigate controlling factors. It is suggested that the continent continent collision is the dominant factor controlling the Chinese tectonic stress field. The ocean continent convergence along the subduction zone is an important factor. There exists tensile boundary force along the marginal sea.

Molecular modeling of oscillating GHz electric field influence on the kinesin affinity to microtubule

Kinesin is a microtubule-associated motor protein which can respond to the external electric field due to its polarity.Using a molecular dynamics simulation method, the effect of such a field on the affinity of kinesin to the α β-tubulin is investigated in this study. To consider kinesin affinity, the system is exposed to an electric field of 0.03 V/nm with frequency values of 1, 2,..., 9, and 10 GHz. It is found that the applied electric field can change kinesin affinity to the microtubule.These changes could perturb the normal operation of kinesin, such as the processive motility of kinesin on the microtubule.

NUMERICAL SIMULATION BY COMPUTATIONAL FLUID DYNAMICS AND EXPERIMENTAL STUDY ON STIRRED BIOREACTOR WITH PUNCHED IMPELLER

Instantaneous flow field and temperature field of the two-phase fluid are measured by particle image velocimetry （PIV） and steady state method during the state of onflow. A turbulent two-phase fluid model of stirred bioreactor with punched impeller is established by the computational fluid dynamics （CFD）, using a rotating coordinate system and sliding mesh to describe the relative motion between impeller and baffles. The simulation and experiment results of flow and temperature field prove their warps are less than 10% and the mathematic model can well simulate the fields, which will also provide the study on optimized-design and scale-up of bioreactors with reference value.

As Regards the Speed in a Medium of the Electromagnetic Radiation Field

The velocity of the electromagnetic radiation in a perfect dielectric, containing no charges and no conduction currents, is explored and determined on making use of the Lorentz transformations. Besides the idealised blackbody radiation, whose vacuum propagation velocity is the universal constant c, being this value independent of the observer, there is another behaviour of electromagnetic radiation, we call inertial radiation, which is characterized by an electromagnetic inertial density , and therefore, it happens to be described by a time-like Poynting four-vector field which propagates with velocity . is found to be a relativistic invariant expressible in terms of the relativistic invariants of the electromagnetic field. It is shown that there is a rest frame, where the Poynting vector is equal to zero. Both phase and group velocities of the electromagnetic radiation are evaluated. The wave and eikonal equations for the dynamics of the radiation field are formulated.

Bioactive superparamagnetic iron oxide-gold nanoparticles regulated by a dynamic magnetic field induce neuronal Ca^(2+)influx and differentiation

Treating neurodegenerative diseases,e.g.,Alzheimer’s Disease,remains a significant challenge due to the limited neuroregeneration rate in the brain.The objective of this study is to evaluate the hypothesis that external magnetic field(MF)stimulation of nerve growth factor functionalized superparamagnetic iron oxide-gold(NGF-SPIO-Au)nanoparticles(NPs)can induce Ca^(2+)influx,membrane depolarization,and enhance neuron differentiation with dynamic MF(DMF)outperforming static MF(SMF)regulation.We showed the that total intracellular Ca^(2+)influx of PC-12 cells was improved by 300%and 535%by the stimulation of DMF(1 Hz,0.5 T,30min)with NGF-SPIO-Au NPs compared to DMF alone and SMF with NGF-SPIO-Au NPs,respectively,which was attributed to successive membrane depolarization.Cellular uptake performed with the application of sodium azide proved that DMF enhanced cellular uptake of NGF-SPIO-Au NPs via endocytosis.In addition,DMF upregulated both the neural differentiation marker(β3-tubulin)and the cell adhesive molecule(integrin-β1)with the existence of NGF-SPIO-Au NPs,while SMF did not show these effects.The results imply that noninvasive DMF-stimulated NPs can regulate intracellular Ca^(2+)influx and enhance neuron differentiation and neuroregeneration rate.

Microscopic mechanism study and process optimization of dimethyl carbonate production coupled biomass chemical looping gasification system

Biomass chemical looping gasification technology is one of the essential ways to utilize abundant biomass resources.At the same time,dimethyl carbonate can replace phosgene as an environmentfriendly organic material for the synthesis of polycarbonate.In this paper,a novel system coupling biomass chemical looping gasification with dimethyl carbonate synthesis with methanol as an intermediate is designed through microscopic mechanism analysis and process optimization.Firstly,reactive force field molecular dynamics simulation is performed to explore the reaction mechanism of biomass chemical looping gasification to determine the optimal gasification temperature range.Secondly,steady-state simulations of the process based on molecular dynamics simulation results are carried out to investigate the effects of temperature,steam to biomass ratio,and oxygen carrier to biomass ratio on the syngas yield and compositions.In addition,the main energy indicators of biomass chemical looping gasification process including lower heating value and cold gas efficiency are analyzed based on the above optimum parameters.Then,two synthesis stages are simulated and optimized with the following results obtained:the optimal temperature and pressure of methanol synthesis stage are 150℃ and 4 MPa;the optimal temperature and pressure of dimethyl carbonate synthesis stage are 140℃ and 0.3 MPa.Finally,the pre-separation-extraction-decantation process separates the mixture of dimethyl carbonate and methanol generated in the synthesis stage with 99.11%purity of dimethyl carbonate.Above results verify the feasibility of producing dimethyl carbonate from the perspective of multi-scale simulation and realize the multi-level utilization of biomass resources.