Joint inversion of gravity and magnetic data for a two-layer model

Based on the synchronous joint gravity and magnetic inversion of single interface by Pilkington and the need of revealing Cenozoic and crystalline basement thickness in the new round of oil-gas exploration, we propose a joint gravity and magnetic inversion methodfor two-layer models by concentrating on the relationship between the change of thicknessI and position of the middle layer and anomaly and discuss the effects of the key parameters. Model tests and application to field data show the validity of this method.

Anisotropic time-dependent behaviors of shale under direct shearing and associated empirical creep models

Understanding the anisotropic creep behaviors of shale under direct shearing is a challenging issue.In this context,we conducted shear-creep and steady-creep tests on shale with five bedding orientations (i.e.0°,30°,45°,60°,and 90°),under multiple levels of direct shearing for the first time.The results show that the anisotropic creep of shale exhibits a significant stress-dependent behavior.Under a low shear stress,the creep compliance of shale increases linearly with the logarithm of time at all bedding orientations,and the increase depends on the bedding orientation and creep time.Under high shear stress conditions,the creep compliance of shale is minimal when the bedding orientation is 0°,and the steady-creep rate of shale increases significantly with increasing bedding orientations of 30°,45°,60°,and 90°.The stress-strain values corresponding to the inception of the accelerated creep stage show an increasing and then decreasing trend with the bedding orientation.A semilogarithmic model that could reflect the stress dependence of the steady-creep rate while considering the hardening and damage process is proposed.The model minimizes the deviation of the calculated steady-state creep rate from the observed value and reveals the behavior of the bedding orientation's influence on the steady-creep rate.The applicability of the five classical empirical creep models is quantitatively evaluated.It shows that the logarithmic model can well explain the experimental creep strain and creep rate,and it can accurately predict long-term shear creep deformation.Based on an improved logarithmic model,the variations in creep parameters with shear stress and bedding orientations are discussed.With abovementioned findings,a mathematical method for constructing an anisotropic shear creep model of shale is proposed,which can characterize the nonlinear dependence of the anisotropic shear creep behavior of shale on the bedding orientation.

Investigating the elliptic anisotropy of identified particles in p-Pb collisions with a multi-phase transport model

The elliptic azimuthal anisotropy coefficient(v_(2))of the identified particles at midrapidity(|η|<0.8)was investigated in p-Pb collisions at√s_(NN)=5.02 TeV using a multi-phase transport model(AMPT).The calculations of differential v_(2)based on the advanced flow extraction method of light flavor hadrons(pions,kaons,protons,andΛ)in small collision systems were extended to a wider transverse momentum(p_(T))range of up to 8 GeV/c for the first time.The string-melting version of the AMPT model provides a good description of the measured p_(T)-differential v_(2)of the mesons but exhibits a slight deviation from the baryon v_(2).In addition,we observed the features of mass ordering at low p_(T)and the approximate number-of-constituentquark(NCQ)scaling at intermediate p_(T).Moreover,we demonstrate that hadronic rescattering does not have a significant impact on v_(2)in p-Pb collisions for different centrality selections,whereas partonic scattering dominates in generating the elliptic anisotropy of the final particles.This study provides further insight into the origin of collective-like behavior in small collision systems and has referential value for future measurements of azimuthal anisotropy.

Asymptotic solution of a weak nonlinear model for the mid-latitude stationary wind field of a two-layer barotropic ocean

A weak nonlinear model of a two-layer barotropic ocean with Rayleigh dissipation is built.The analytic asymptotic solution is derived in the mid-latitude stationary wind field,and the physical meaning of the corresponding problem is discussed.

Effect of the Coefficient on the Performance of A Two-Layer Boussinesq- Type Model

The coefficients embodied in a Boussinesq-type model are very important since they are determined to optimize the linear and nonlinear properties.In most conventional Boussinesq-type models,these coefficients are assigned the specific values.As for the multi-layer Boussinesq-type models with the inclusion of the vertical velocity,however,the effect of the different values of these coefficients on linear and nonlinear performances has never been investigated yet.The present study focuses on a two-layer Boussinesq-type model with the highest spatial derivatives being 2 and theoretically and numerically examines the effect of the coefficient on model performance.Theoretical analysis show that different values for(0.13≤α≤0.25)do not have great effects on the high accuracy of the linear shoaling,linear phase celerity and even third-order nonlinearity for water depth range of 0<kh≤10(k is wave number and h is water depth).The corresponding errors using different values are restricted within 0.1%,0.1%and 1%for the linear shoaling amplitude,dispersion and nonlinear harmonics,respectively.Numerical tests including regular wave shoaling over mildly varying slope from deep to shallow water,regular wave propagation over submerged bar,bichromatic wave group and focusing wave propagation over deep water are conducted.The comparison between numerical results using different values of,experimental data and analytical solutions confirm the theoretical analysis.The flexibility and consistency of the two-layer Boussinesq-type model is therefore demonstrated theoretically and numerically.

Anisotropic dynamic permeability model for porous media

Based on the tortuous capillary network model,the relationship between anisotropic permeability and rock normal strain,namely the anisotropic dynamic permeability model(ADPM),was derived and established.The model was verified using pore-scale flow simulation.The uniaxial strain process was calculated and the main factors affecting permeability changes in different directions in the deformation process were analyzed.In the process of uniaxial strain during the exploitation of layered oil and gas reservoirs,the effect of effective surface porosity on the permeability in all directions is consistent.With the decrease of effective surface porosity,the sensitivity of permeability to strain increases.The sensitivity of the permeability perpendicular to the direction of compression to the strain decreases with the increase of the tortuosity,while the sensitivity of the permeability in the direction of compression to the strain increases with the increase of the tortuosity.For layered reservoirs with the same initial tortuosity in all directions,the tortuosity plays a decisive role in the relative relationship between the variations of permeability in all directions during pressure drop.When the tortuosity is less than 1.6,the decrease rate of horizontal permeability is higher than that of vertical permeability,while the opposite is true when the tortuosity is greater than 1.6.This phenomenon cannot be represented by traditional dynamic permeability model.After the verification by experimental data of pore-scale simulation,the new model has high fitting accuracy and can effectively characterize the effects of deformation in different directions on the permeability in all directions.

Two-Layer Non-Hydrostatic Model for Generation and Propagation of Interfacial Waves

When pycnocline thickness of ocean density is relatively small, density stratification can be well represented as a two-layer system. In this article, a depth integrated model of the two-layer fluid with constant density is considered,and a variant of the edge-based non-hydrostatic numerical scheme is formulated. The resulting scheme is very efficient since it resolves the vertical fluid depth only in two layers. Despite using just two layers, the numerical dispersion is shown to agree with the analytical dispersion curves over a wide range of kd, where k is the wave number and d the water depth. The scheme was tested by simulating an interfacial solitary wave propagating over a flat bottom, as well as over a bottom step. On a laboratory scale, the formation of an interfacial wave is simulated,which also shows the interaction of wave with a triangular bathymetry. Then, a case study using the Lombok Strait topography is discussed, and the results show the development of an interfacial wave due to a strong current passing through a sill.

Cooperative Caching Strategy Based on Two-Layer Caching Model for Remote Sensing Satellite Networks

In Information Centric Networking(ICN)where content is the object of exchange,in-network caching is a unique functional feature with the ability to handle data storage and distribution in remote sensing satellite networks.Setting up cache space at any node enables users to access data nearby,thus relieving the processing pressure on the servers.However,the existing caching strategies still suffer from the lack of global planning of cache contents and low utilization of cache resources due to the lack of fine-grained division of cache contents.To address the issues mentioned,a cooperative caching strategy(CSTL)for remote sensing satellite networks based on a two-layer caching model is proposed.The two-layer caching model is constructed by setting up separate cache spaces in the satellite network and the ground station.Probabilistic caching of popular contents in the region at the ground station to reduce the access delay of users.A content classification method based on hierarchical division is proposed in the satellite network,and differential probabilistic caching is employed for different levels of content.The cached content is also dynamically adjusted by analyzing the subsequent changes in the popularity of the cached content.In the two-layer caching model,ground stations and satellite networks collaboratively cache to achieve global planning of cache contents,rationalize the utilization of cache resources,and reduce the propagation delay of remote sensing data.Simulation results show that the CSTL strategy not only has a high cache hit ratio compared with other caching strategies but also effectively reduces user request delay and server load,which satisfies the timeliness requirement of remote sensing data transmission.

Forward modeling of marine DC resistivity method for a layered anisotropic earth

Since the ocean bottom is a sedimentary environment wherein stratification is well developed, the use of an anisotropic model is best for studying its geology. Beginning with Maxwell＇s equations for an anisotropic model, we introduce scalar potentials based on the divergence-free characteristic of the electric and magnetic （EM） fields. We then continue the EM fields down into the deep earth and upward into the seawater and couple them at the ocean bottom to the transmitting source. By studying both the DC apparent resistivity curves and their polar plots, we can resolve the anisotropy of the ocean bottom. Forward modeling of a high-resistivity thin layer in an anisotropic half-space demonstrates that the marine DC resistivity method in shallow water is very sensitive to the resistive reservoir but is not influenced by airwaves. As such, it is very suitable for oil and gas exploration in shallowwater areas but, to date, most modeling algorithms for studying marine DC resistivity are based on isotropic models. In this paper, we investigate one-dimensional anisotropic forward modeling for marine DC resistivity method, prove the algorithm to have high accuracy, and thus provide a theoretical basis for 2D and 3D forward modeling.

Finite-difference numerical modeling with even-order accuracy in two-phase anisotropic media

To improve the accuracy of the conventional finite-difference method, finitedifference numerical modeling methods of any even-order accuracy are recommended. We introduce any even-order accuracy difference schemes of any-order derivatives derived from Taylor series expansion. Then, a finite-difference numerical modeling method with any evenorder accuracy is utilized to simulate seismic wave propagation in two-phase anisotropic media. Results indicate that modeling accuracy improves with the increase of difference accuracy order number. It is essential to find the optimal order number, grid size, and time step to balance modeling precision and computational complexity. Four kinds of waves, static mode in the source point, SV wave cusps, reflection and transmission waves are observed in two-phase anisotropic media through modeling.

Numerical modeling of PP- and PS-wave azimuthal anisotropy in HTI media

This study focuses on the factors that may affect the feasibility of performing elliptical anisotropy analysis on azimuthal PP and PS-wave data in HTI media, with the aim of using the modeling results as guidance in real seismic data application. Our results reveal that there is an offset limitation for both PP- and PS-waves in elliptical anisotropy fitting, and that PS-waves show a wider applicable offset range and larger observable azimuthal anisotropy than PP-waves. The major axis of the elliptical fit to the amplitudes of the R-component is perpendicular to the fracture strike, which is opposite to that in PP-wave analysis. The azimuthal interval travel time of PS-waves shows a nearly elliptical distribution and the major axis of the fit ellipse is perpendicular to the fracture strike, which is same as that in PP-wave analysis. For data within the applicable offset range, the anisotropic magnitude obtained from amplitude and travel time attributes of PP- and PS-waves exhibits a dependence on fracture density, and the major to minor axis ratio of the fit ellipse may be used to infer the relative distribution of fracture densities.

Brittleness index and seismic rock physics model for anisotropic tight-oil sandstone reservoirs

Brittleness analysis becomes important when looking for sweet spots in tightoil sandstone reservoirs. Hence, appropriate indices are required as accurate brittleness evaluation criteria. We construct a seismic rock physics model for tight-oil sandstone reservoirs with vertical fractures. Because of the complexities in lithology and pore structure and the anisotropic characteristics of tight-oil sandstone reservoirs, the proposed model is based on the solid components, pore connectivity, pore type, and fractures to better describe the sandstone reservoir microstructure. Using the model, we analyze the brittleness sensitivity of the elastic parameters in an anisotropic medium and establish a new brittleness index. We show the applicability of the proposed brittleness index for tight-oil sandstone reservoirs by considering the brittleness sensitivity, the rock physics response characteristics, and cross-plots. Compared with conventional brittleness indexes, the new brittleness index has high brittleness sensitivity and it is the highest in oil-bearing brittle zones with relatively high porosity. The results also suggest that the new brittleness index is much more sensitive to elastic properties variations, and thus can presumably better predict the brittleness characteristics of sweet spots in tight-oil sandstone reservoirs.

Study on Strength and Life of Anisotropic Single Crystal Blade - Part I: Crystallographic Constitutive Models and Applications

The single crystal blade is one of the key technologies for improving the performance, durability and reliability of aero-engines and ground gas-turbine engines. However, the anisotropic mechanical properties of the single crystal material makes a great deal of difficulties on the development and the application of the single crystal blade, which is a challenge for the engineering application of the single crystal superalloy and the theoretic bases of the application. Some researches on the strength analysis and the life prediction of the anisotropic single crystal blade were carried out by the authors' research team. They are as follows. The crystallographic constitutive models for the plastic and the creep behaviors and the method of the rupture life prediction were established and verified. The tensile or the creep experiments for DD3 single crystal alloy with different orientations under different temperatures and different tensile rates or under different temperatures and different stress levels were carried out. The experimental data and the anisotropic properties at intermediate and high temperatures revealed by the experiments are significant for the application of the single crystal alloy. In addition, the experimental research for a kind of single crystal blade was also made. As the application of the researches the strength analysis and the life prediction were carried out for the single crystal blade of a certain aero-engine. In this part, the constitutive models and their applications are described, and the experimental research work will be described in part II.

Critical behaviours and magnetic properties of three-dimensional bond and anisotropy dilution Blume-Capel model in the presence of an applied field

This paper studies the critical behaviours and magnetic properties of three-dimensional bond and anisotropy dilution Blume-Capel model （BCM） in the presence of an applied field within the effective field theory. The trajectory of tricritical point, reentrant transitions and degenerate patterns of anisotropy are obtained both for the bond and the anisotropy dilutions. The global phase diagrams demonstrate unusually reentrant phenomena. The temperature dependences of magnetization curves undergo remarkable spin glass behaviour at low temperatures, and transform from ferromagnetism to paramagnetism at high temperature in applied fields. Temperature dependence of magnetic susceptibility curve is in qualitative agreement with experimental result.

A physical model study of shale seismic responses and anisotropic inversion

The seismic responses of the shale properties are critical for shale gas reservoir evaluation and production.It has been widely reported that the clay minerals have substantial influences on the seismic wave anisotropy and brittleness.Hence,knowing the seismic responses of the clay-rich shales and estimation of shale elastic properties are significant for the shale gas industry.A physical model containing two groups of shale blocks as the target formations is constructed in laboratory.The group S contains six shale blocks with different clay contents,and the group N contains six shale blocks with different porosity.The acquired 2D seismic data is used to analyze the seismic responses of two corresponding seismic lines.An anisotropic three-term inversion method is applied to one of the seismic inline to estimate the elastic properties the target shale blocks.The inversed attributes can be used to reveal the effects of shale clay contents.This study shows the substantial significance of using a physical model to observe the seismic responses of shale properties.The inversion results indicate that the anisotropic three-term inversion method could provide accurate results of elastic properties as well as the P-wave anisotropy parameter for shale formations.

Evaluating the mechanical properties of anisotropic shale containing bedding and natural fractures with discrete element modeling

Natural fracture data from one of the Carboniferous shale masses in the eastern Qaidam Basin were used to establish a stochastic model of a discrete fracture network and to perform discrete element simulation research on the size efect and mechanical parameters of shale.Analytical solutions of fctitious joints in transversely isotropic media were derived,which made it possible for the proposed numerical model to simulate the bedding and natural fractures in shale masses.The results indicate that there are two main factors infuencing the representative elementary volume(REV)size of a shale mass.The frst and most decisive factor is the presence of natural fractures in the block itself.The second is the anisotropy ratio:the greater the anisotropy is,the larger the REV.The bedding angle has little infuence on the REV size,whereas it has a certain infuence on the mechanical parameters of the rock mass.When the bedding angle approaches the average orientation of the natural fractures,the mechanical parameters of the shale blocks decrease greatly.The REV representing the mechanical properties of the Carboniferous shale masses in the eastern Qaidam Basin were comprehensively identifed by considering the infuence of bedding and natural fractures.When the numerical model size is larger than the REV,the fractured rock mass discontinuities can be transformed into equivalent continuities,which provides a method for simulating shale with natural fractures and bedding to analyze the stability of a borehole wall in shale.

Prediction of brittleness based on anisotropic rock physics model for kerogen-rich shale

The construction of a shale rock physics model and the selection of an appropriate brittleness index （B/） are two significant steps that can influence the accuracy of brittleness prediction. On one hand, the existing models of kerogen-rich shale are controversial, so a reasonable rock physics model needs to be built. On the other hand, several types of equations already exist for predicting the BI whose feasibility needs to be carefully considered. This study constructed a kerogen-rich rock physics model by performing the self- consistent approximation and the differential effective medium theory to model intercoupled clay and kerogen mixtures. The feasibility of our model was confirmed by comparison with classical models, showing better accuracy. Templates were constructed based on our model to link physical properties and the BL Different equations for the BI had different sensitivities, making them suitable for different types of formations. Equations based on Young＇s Modulus were sensitive to variations in lithology, while those using Lame＇s Coefficients were sensitive to porosity and pore fluids. Physical information must be considered to improve brittleness prediction.

Reducing the anisotropy of a Brazilian disc generated in a bonded-particle model

The Brazilian test is a widely used method for determining the tensile strength of rocks and for calibrating parameters in bondedparticle models（BPMs）. In previous studies, the Brazilian disc has typically been trimmed from a compacted rectangular specimen. The present study shows that different tensile strength values are obtained depending on the compressive loading direction. Several measures are proposed to reduce the anisotropy of the disc. The results reveal that the anisotropy of the disc is significantly influenced by the compactibility of the specimen from which it is trimmed. A new method is proposed in which the Brazilian disc is directly generated with a particle boundary, effectively reducing the anisotropy. The stiffness（particle and bond） and strength（bond） of the boundary are set at less than and greater than those of the disc assembly, respectively,which significantly decreases the stress concentration at the boundary contacts and prevents breakage of the boundary particle bonds. This leads to a significant reduction in the anisotropy of the disc and the discreteness of the tensile strength. This method is more suitable for carrying out a realistic Brazilian test for homogeneous rock-like material in the BPM.

An anisotropic constitutive model of geomaterials based on true triaxial testing and its application

Series of testing on coarse grained soils were carried out with a true triaxial testing apparatus. The loads were applied from the major principal and minor principal directions, respectively, to simulate the construction and water impounding process of a rock fill dam. The stress and strain relationships induced by the different loading methods were investigated. A remarkable stress-induced anisotropy under complex stress state was observed. Contrary to popular assumptions in traditional numerical analysis and constitutive models, it was found that different elastic modulus and Poisson ratio exist in different principal directions in rock fill dams. From the testing results, an anisotropic constitutive model based on Duncan-Chang nonlinear model is presented to overcome the limitations of axi-symmetric assumptions in conventional triaxial experiments and constitutive models. Both models were then applied in FEM analysis of an under-construction earth core high rock soil filled dam with the focus on hydraulic fracturing. The study reveals the major biases that exist when numerical analysis and constitutive models do not give serious consideration to the intermediate principal stress and anisotropy effects in soil rock built structures.

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.