Mechanisms of salt rejection at the ice-liquid interface during the freezing of pore fluids in the seasonal frozen soil area

Seasonal frozen soil accounts for about 53.50%of the land area in China.Frozen soil is a complex multiphase system where ice,water,soil,and air coexist.The distribution and migration of salts in frozen soil during soil freezing are notably different from those in unfrozen soil areas.However,little knowledge is available about the process and mechanisms of salt migration in frozen soil.This study explores the mechanisms of salt migration at the ice-liquid interface during the freezing of pore fluids through batch experiments.The results are as follows.The solute concentrations of liquid and solid phases at the ice-liquid interface(C*_(L),C*_(S))gradually increased at the initial stage of freezing and remained approximately constant at the middle stage.As the ice-liquid interface advanced toward the system boundary,the diffusion of the liquid phase was blocked but the ice phase continued rejecting salts.As a result,C*_(L)and C*_(S)rapidly increased at the final stage of freezing.The distribution characteristics of solutes in ice and the liquid phases before C*_(L)and C*_(S)became steady were mainly affected by the freezing temperature,initial concentrations,and particle-size distribution of media(quartz sand and kaolin).In detail,the lower the freezing temperature and the better the particle-size distribution of media,the higher the solute proportion in the ice phase at the initial stage of freezing.Meanwhile,the increase in concentration first promoted but then inhibited the increase of solutes in the ice phase.These results have insights and scientific significance for the tackling of climate change,the environmental protection of groundwater and soil,and infrastructure protection such as roads,among other things.

Linear numerical calculation method for obtaining critical point,pore fluid,and framework parameters of gas-bearing media

Up to now, the primary method for studying critical porosity and porous media are experimental measurements and data analysis. There are few references on how to numerically calculate porosity at the critical point, pore fluid-related parameters, or framework-related parameters. So in this article, we provide a method for calculating these elastic parameters and use this method to analyze gas-bearing samples. We first derive three linear equations for numerical calculations. They are the equation of density p versus porosity Ф, density times the square of compressional wave velocity p Vp^2 versus porosity, and density times the square of shear wave velocity pVs^2 versus porosity. Here porosity is viewed as an independent variable and the other parameters are dependent variables. We elaborate on the calculation steps and provide some notes. Then we use our method to analyze gas-bearing sandstone samples. In the calculations, density and P- and S-velocities are input data and we calculate eleven relative parameters for porous fluid, framework, and critical point. In the end, by comparing our results with experiment measurements, we prove the viability of the method.

The effect of pore fluid on seismicity: a computer model

The influence of fluid on seismicity of a computerized system is analyzed in this paper. The diffusion equation of fluid in a crustal fault area is developed and used in the calculation of a spring slide damper model. With mirror imagin boundary condition and three initial conditions, the equation is solved for a dynamic model that consists of six seismic belts and eight seismogenous sources in each belt with both explicit algorithm and implicit algorithm. The analysis of the model with water sources shows that the implicit algorithm is better to be used to calculate the model. Taking a constant proportion of the pore pressure of a broken element to that of its neighboring elements, the seismicity of the model is calculated with mirror boundary condition and no water source initial condition. The results shows that the frequency and magnitude of shocks are both higher than those in the model with no water pore pressure, which provides more complexity to earthquake prediction.

Evolution of pore fluid of low porosity and permeability reservoir in Qinnan Sag:a case study of 29-2 structurein Qinhuangdao area

The southeast depression of Qinnan Sag is a potential oil-gas exploration region in the Bohai Sea area. With the analysis of large quantity of rock thin sections,scanning electron microscope and the physical property data of reservoir,the authors studied the petrological characteristics and the evolution of pore and fluid of sandstone in the deeper strata in 29-2 structure in Qinhuangdao area. The results show that the evolutionary tendency of Paleogene sandstone reservoir porous fluid in research area is changed from alkaline porous fluid to acidic porous fluid,and back to alkaline porous fluid. There are three stages of reservoir porous evolution in Qinhuangdao area,namely sharp decrease in porosity due to mechanical compaction,increase in porosity because of corrosion and dissolution,and remarkable reduction owing to carbonate cementation.

Conductivity characteristics of landslide considering porosity,saturation,temperature and ion concentration

High-density resistivity imaging method is widely used in landslide monitoring.The resistivity of rock and soil is closely related to factors,such as porosity,moisture content,saturation and temperature.In this study,the resistivity test was designed to investigate the influence of physical factors and pore solution components on the resistivity of landslide soil.Experimental and analytical results find that both moisture content and volumetric water content varies greatly under the same resistivity.At different temperatures,soil resistivity exhibits great changes.Under the same temperature,the ion concentration and species in pore solutions have great influence on soil resistivity.Based on the test results and grey correlation analysis,this study established a resistivity model by considering porosity,saturation,temperature and ion concentration.The study lays a foundation for the high-density resistivity method to measure the moisture content of landslides.

Fluid sensitivity study of elastic parameters in lowmedium porosity and permeability reservoir rocks

In this article,based on the acoustic measurements of core samples obtained from the low to medium porosity and permeability reservoirs in the WXS Depression,the densities and P and S wave velocities of these core samples were obtained.Then based on these data,a series of elastic parameters were computed.From the basic theory and previous pore fluid research results,we derived a new fluid identification factor（F）.Using the relative variations,Ag/w and Ao/w,of the elastic parameters between gas and water saturated samples and between oil and water saturated samples,λρ,σHSFIF,Kρ,λρ-2μρ,and F as quantitative indicators,we evaluate the sensitivity of the different fluid identification factors to identify reservoir fluids and validate the effects by crossplots.These confirm that the new fluid identification factor（F） is more sensitive for distinguishing oil and water than the traditional method and is more favorable for fliud identification in low to medium porosity and permeability reservoirs.

Seismic reflection and transmission coefficients of a single layer sandwiched between two dissimilar poroelastic solids

The seismic reflection and transmission characteristics of a single layer sandwiched between two dissimilar poroelastic solids saturated with two immiscible viscous fluids are investigated. The sandwiched layer is modeled as a porous solid with finite thickness. The propagation of waves is represented with potential functions. The displacements of particles in different phases of the aggregate are defined in terms of these potential functions. Due to the presence of viscosity in pore fluids, the reflected and transmitted waves are inhomogeneous in nature, i.e., with different directions of propagation and attenuation. The closed-form analytical expressions for reflection and transmission coefficients are derived theoretically for appropriate boundary conditions. These expressions are calculated as a non-singular system of linear algebraic equations and depend on the various parameters involved in this non-singular system. Hence,numerical examples are studied to determine the effects of various properties of the sandwich layer on reflection and transmission coefficients. The essential features of layer thickness, incident direction, wave frequency, liquidsaturation and capillary pressure of the porous layer on reflection and transmission coefficients are depicted graphically and discussed. The analysis shows that reflection and transmission coefficients are strongly associated with incident direction and various properties of the porous layer.

Brittleness characteristics of tight oil siltstones

Rock brittleness directly affects reservoir fracturing and its evaluation is essential for establishing fracturing conditions prior to reservoir reforming. Dynamic and static brittleness data were collected from siltstones of the Qingshankou Formation in Songliao Basin. The brittle-plastic transition was investigated based on the stress-strain relation. The results suggest that the brittleness indices calculated by static elastic parameters are negatively correlated with the stress drop coefficient and the brittleness index B2, defined as the average of the normalized Young＇s modulus and Poisson＇s ratio, is strongly correlated with the stress drop. The brittleness index B2, Young＇s modulus, and Poisson＇s ratio correlate with the brittle minerals content; that is, quartz, carbonates, and pyrite. We also investigated the correlation between pore fluid and porosity and dynamic brittle characteristic based on index B2. Pore fluid increases the plasticity of rock and reduces brittleness; moreover, with increasing porosity, rock brittleness decreases. The gas-saturated siltstone brittleness index is higher than that in oil- or water-saturated siltstone; the difference in the brittleness indices of oil- and water-saturated siltstone is very small. By comparing the rock mechanics and ultrasonic experiments, we find that the brittleness index obtained from the rock mechanics experiments is smaller than that obtained from the ultrasonic experiments; nevertheless, both decrease with increasing porosity as well as their differences. Ultrasonic waves propagate through the rock specimens without affecting them, whereas rock mechanics experiments are destructive and induce microcracking and porosity increases; consequently, the brittleness of low-porosity rocks is affected by the formation of internal microcrack systems.

Quantitative evaluation of geological fluid evolution and accumulated mechanism:in case of tight sandstone gas field in central Sichuan Basin

Tight gas exploration plays an important part in China’s unconventional energy strategy.The tight gas reservoirs in the Jurassic Shaximiao Formation in the Qiulin and Jinhua Gas Fields of central Sichuan Basin are characterized by shallow burial depths and large reserves.The evolution of the fluid phases is a key element in understanding the accumulation of hydrocarbons in tight gas reservoirs.This study investigates the fluid accumulation mechanisms and the indicators of reservoir properties preservation and degradation in a tight gas reservoir.Based on petrographic observations and micro-Raman spectroscopy,pure CH4 inclusions,pure CO2 inclusions,hybrid CH4–CO2 gas inclusions,and N2-rich gas inclusions were studied in quartz grains.The pressure–volume–temperature–composition properties(PVT-x)of the CH4 and CO2 bearing inclusions were determined using quantitative Raman analysis and thermodynamic models,while the density of pure CO2 inclusions was calculated based on the separation of Fermi diad.Two stages of CO2 fluid accumulation were observed:primary CO2 inclusions,characterized by higher densities(0.874–1.020 g/cm3)and higher homogenization temperatures(>210°C)and secondary CO2 inclusions,characterized by lower densities(0.514–0.715 g/cm3)and lower homogenization temperatures:~180–200°C).CO2 inclusions with abnormally high homogenization temperatures are thought to be the result of deep hydrothermal fluid activity.The pore fluid pressure(44.0–58.5 MPa)calculated from the Raman shift of C–H symmetric stretching(v1)band of methane inclusions is key to understanding the development of overpressure.PT entrapment conditions and simulation of burial history can be used to constrain the timing of paleo-fluid emplacement.Methane accumulated in the late Cretaceous(~75–65 Ma),close to the maximum burial depth during the early stages of the Himalayan tectonic event while maximum overpressure occurred at~70 Ma,just before uplift.Later,hydrocarbon gas migrated through the faults and gradually displaced the early emplaced CO2 in the reservoirs accompanied by a continuous decrease in overpressure during and after the Himalayan event,which has led to a decrease in the reservoir sealing capabilities.The continuous release of overpressure to present-day conditions indicates that the tectonic movement after the Himalayan period has led to a decline in reservoir conditions and sealing properties.

Study on the Triggering Effect of the October, 2005 Pakistan M_W7.6 Earthquake on Its Aftershocks

The influences upon aftershocks of Coulomb failure stress change (CFSC) generated by the main-shock of the October 8, 2005, Pakistan earthquake are calculated and analyzed. The following factors are included in the calculation: (1) the difference between the pore fluid pressure and the medium elastic constant in the fault plane area and those of its surrounding medium; (2) the tectonic stress direction of the seismic source area; (3) the aftershock failure mechanism of aftershocks is calculated by stacking the tectonic stress with the stress change generated by the main-shock. Our study, which includes many factors, fits fairly well with the aftershock distribution. It indicates that most of the aftershocks were triggered by the Pakistan main-shock that occurred on October 8, 2005.

Lattice Boltzmann simulation of fluid flow through coal reservoir's fractal pore structure

The influences of fractal pore structure in coal reservoir on coalbed methane(CBM) migration were analyzed in detail by coupling theoretical models and numerical methods.Different types of fractals were generated based on the construction thought of the standard Menger Sponge to model the 3D nonlinear coal pore structures.Then a correlation model between the permeability of fractal porous medium and its pore-size-distribution characteristics was derived using the parallel and serial modes and verified by Lattice Boltzmann Method(LBM).Based on the coupled method,porosity(ф),fractal dimension of pore structure(Db),pore size range(rmin,rmax) and other parameters were systematically analyzed for their influences on the permeability(ф) of fractal porous medium.The results indicate that:① the channels connected by pores with the maximum size(rmax) dominate the permeability,approximating in the quadratic law;② the greater the ratio of r max and r min is,the higher is;③ the relationship between D b and follows a negative power law model,and breaks into two segments at the position where Db ≌2.5.Based on the results above,a predicting model of fractal porous medium permeability was proposed,formulated as k=cfrnmax,where C and n(approximately equal to 2) are constants and f is an expression only containing parameters of fractal pore structure.In addition,the equivalence of the new proposed model for porous medium and the Kozeny-Carman model k=Crn was verified at Db =2.0.

Uncoupled state space solution to layered poroelastic medium with anisotropic permeability and compressible pore fluid

This paper presents an uncoupled state space solution to three-dimensional consolidation of layered poroelastic medium with anisotropic permeability and compressible pore fluid.Starting from the basic equations of poroelastic medium,and introducing intermediate variables,the state space equation usually comprising eight coupled state vectors is uncoupled into two sets of equations of six and two state vectors in the Laplace-Fourier transform domain.Combined with the continuity conditions between adjacent layers and boundary conditions,the uncoupled state space solution of a layered poroelastic medium is obtained by using the transfer matrix method.Numerical results show that the anisotropy of permeability and the compressibility of pore fluid have remarkable influence on the consolidation behavior of poroelastic medium.

ACOUSTOELASTIC THEORY FOR FLUID-SATURATED POROUS MEDIA

Based on the finite deformation theory of the continuum and poroelastic theory, the aeoustoelastic theory for fluid-saturated porous media （FSPM） in natural and initial coordi- nates is developed to investigate the influence of effective stresses and fluid pore pressure on wave velocities. Firstly, the assumption of a small dynamic motion superimposed on a largely static pre- deformation of the FSPM yields natural, initial, and final configurations, whose displacements, strains, and stresses of the solid-skeleton and the fluid in an FSPM particle could be described in natural and initial coordinates, respectively. Secondly, the subtraction of initial-state equations of equilibrium from the final-state equations of motion and the introduction of non-linear constitu- rive relations of the FSPM lead to equations of motion for the small dynamic motion. Thirdly, the consideration of homogeneous pre-deformation and the plane harmonic form of the small dynamic motion gives an acoustoelastic equation, which provides analytical formulations for the relation of the fast longitudinal wave, the fast shear wave, the slow shear wave, and the slow longitudinal wave with solid-skeleton stresses and fluid pore-pressure. Lastly, an isotropic FSPM under the close-pore jacketed condition, open-pore jacketed condition, traditional unjacketed condition, and triaxial condition is taken as an example to discuss the velocities of the fast and slow shear waves propagating along the direction of one of the initial principal solid-skeleton strains. The detailed discussion shows that the wave velocities of the FSPM are usually influenced by the effective stresses and the fluid pore pressure. The fluid pore-pressure has little effect on the wave velocities of the FSPM only when the components of the applied initial principal solid-skeleton stresses or strains are equal, which is consistent with the previous experimental results.