Clogging caused by coupled grain migration and compaction effect during groundwater recharge for unconsolidated sandstone reservoir in groundwater-source heat pump

In unconsolidated sandstone reservoirs,presence of numerous movable grains and a complex grain size composition necessitates a clear understanding of the physical clogging process for effective groundwater recharge in groundwater-source heat pump systems.To investigate this,a series of seepage experiments was conducted under in situ stress conditions using unconsolidated sandstone samples with varying grain compositions.The clogging phenomenon arises from the combined effects of grain migration and compaction,wherein the migration of both original and secondary crushed fine-grain particles blocks the seepage channels.Notably,grain composition influences the migration and transport properties of the grains.For samples composed of smaller grains,the apparent permeability demonstrates a transition from stability to decrease.In contrast,samples with larger grains experience a skip at the stability stage and directly enter the decrease stage,with a minor exception of a slight increase observed.Furthermore,a unique failure mode characterized by diameter shrinkage in the upper part of the sample is observed due to the combined effects of grain migration and in situ stress-induced compaction.These testing results contribute to a better understanding of the clogging mechanism caused by the coupled effects of grain migration and compaction during groundwater recharge in unconsolidated sandstone reservoirs used in groundwater-source heat pump systems.

Experimental study on the moving characteristics of fine grains in wide grading unconsolidated soil under heavy rainfall

The initiation mechanism of debris flow is regarded as the key step in understanding the debrisflow processes of occurrence, development and damage. Moreover, migration, accumulation and blocking effects of fine particles in soil will lead to soil failure and then develop into debris flow. Based on this hypothesis and considering the three factors of slope gradient, rainfall duration and rainfall intensity, 16 flume experiments were designed using the method of orthogonal design and completed in a laboratory. Particle composition changes in slope toe, volumetric water content, fine particle movement characteristics and soil failure mechanism were analyzed and understood as follows: the soil has complex, random and unstable structures, which causes remarkable pore characteristics of poor connectivity, non-uniformity and easy variation. The major factors that influence fine particle migration are rainfall intensity and slope. Rainfall intensity dominates particle movement, whereby high intensity rainfall induces a large number of mass movement and sharp fluctuation, causing more fine particles to accumulate at the steep slope toe. The slope toe plays an important role in water collection and fine particleaccumulation. Both fine particle migration and coarse particle movement appears similar fluctuation. Fine particle migration is interrupted in unconnected pores, causing pore blockage and fine particle accumulation, which then leads to the formation of a weak layer and further soil failure or collapses. Fine particle movement also causes debris flow formation in two ways: movement on the soil surface and migration inside the soil. The results verify the hypothesis that the function of fine particle migration in soil failure process is conducive for further understanding the formation mechanism of soil failure and debris flow initiation.

Quantitative description of characteristics of high-capacity channels in unconsolidated sandstone reservoirs using in situ production data

Quantitative description of the high-capacity channels in unconsolidated sandstone reservoirs, into which water was injected to improve oil recovery, is a hot topic in the field of reservoir development. This paper presents a novel approach to describing quantitatively the characteristics of connected macropores in unconsolidated sandstone reservoirs using in situ production data. Based on physical simulation for formation mechanisms of high capacity channels and interwell tracer test data, a mathematical model was established to describe high-capacity channels by grey correlation theory, flow mechanism of fluid in porous media and reservoir engineering, and a program was developed to describe quantitatively the channel characteristics. The predicted results were consistent with field monitoring data （80%）, so this model could be economically and effectively used to identify high-capacity channels.

Numerical simulation of dewatering thick unconsolidated aquifers for safety of underground coal mining

With an increase of mining the upper limits under unconsolidated aquifers, dewatering of the bottom aquifer of the Quaternary system has become a major method to avoid water and sand inrushes.In the 8th District of the Taiping Coal Mine in south-western Shandong province, the bottom aquifer of the Quaternary system is moderate to excellent in water-yielding capacity.The base rock above the coal seam is very thin in the concealed coal field of the Carboniferous and Permian systems.Therefore, a comprehensive dewatering plan from both the ground surface and the panel was proposed to lower the groundwater level in order to ensure mining safety.According to the hydrogeologic conditions of the 8th District, we established a numerical model so that we could simulate the groundwater flow in the dewatering process.We obtained the simulation parameters from previous data using backward modeling, such as the average coefficient of permeability of 12 m/d and the elastic storage coefficient of 0.002.From the same model, we predicted the movement of groundwater and water level variables and obtained the visible effect of the dewatering project.Despite the overburden failure during mining, no water and/or sand inrush occurred because the groundwater level in the bottom aquifer was lowered to a safe water level.

Numerical simulation of surface movement laws under different unconsolidated layers thickness

Based on specific geology and mining conditions of certain coal working face in China, a series of numerical models under different unconsolidated layers thickness were respectively established by employing FLAC3D. The relationship between the unconsolidated layers thickness and surface movement laws was studied. Maximum surface subsidence, Maximum horizontal displacement and surface subsidence degree were obtained. Contours of surface subsidence/horizontal displacement and curves were drawn. Some laws of surface subsidence/horizontal displacement were analyzed. The role of the unconsolidated layers in surface subsidence was revealed. It is significant to predict surface subsidence of thick unconsolidated layers for coal mine and take effective measures to control surface subsidence.

Study on Active Faults in Quaternary Unconsolidated Sediments by Microstructural Method

Active faults in Quaternary unconsolidated sediments are important indicators to determine paleoearthquake events. They can be studied by macroscopic geological survey, but some problems have been encountered, such as invisible active faults. The authors found an approach to solve these problems by microstructural observation. Firstly, oriented original-state samples of Quaternary unconsolidated sediments in active fault zones are collected. Secondly, the samples are consolidated in laboratory. Thirdly, microstructural slides are ground in three-dimension coordinate system. Lastly, microstructures are observed in the microscopic field. By this method, the movement properties of active faults can be determined in lack of the support of macroscopic data. The invisible faults in Quaternary unconsolidated sediments can be found and described. The mechanisms of die-out faults can be also studied. Moreover, the boundaries between different unconsolidated sediments in engineering sites can be studied to judge whether they are active faults or not.

A Novel Simulation Framework for Predicting the Formation Parameters Variation in Unconsolidated Sandstone Reservoir

After long-term waterflooding in unconsolidated sandstone reservoir, the high-permeability channels are easy to evolve, which leads to a significant reduction in water flooding efficiency and a poor oilfield development effect. The current researches on the formation parameters variation are mainly based on the experiment analysis or field statistics, while lacking quantitative research of combining microcosmic and macroscopic mechanism. A network model was built after taking the detachment and entrapment mechanisms of particles in unconsolidated sandstone reservoir into consideration. Then a coupled mathematical model for the formation parameters variation was established based on the network modeling and the model of fluids flowing in porous media. The model was solved by a finite-difference method and the Gauss-Seidel iterative technique. A novel field-scale reservoir numerical simulator was written in Fortran 90 and it can be used to predict 1) the evolvement of high-permeability channels caused by particles release and migration in the long-term water flooding process, and 2) well production performances and remaining oil distribution. In addition, a series of oil field examples with inverted nine-spot pattern was made on the new numerical simulator. The results show that the high-permeability channels are more likely to develop along the main streamlines between the injection and production wells, and the formation parameters variation has an obvious influence on the remaining oil distribution.

Quantitative prediction of oil saturation of unconsolidated sandstone reservoir based on time-lapse seismic “relative difference method”: Taking Zeta oil field in West Africa as an example

In view of the disadvantage that the absolute difference of time-lapse seismic(the difference between monitoring data and base data) is not only related to the change of oil saturation, but also closely related to the thickness of reservoir, a time-lapse seismic "relative difference method"(the ratio of monitoring data to base data) not affected by the thickness of reservoir but only related to the change of fluid saturation, is proposed through seismic forward modeling after fluid displacement simulation. Given the same change of fluid saturation, the absolute difference of time-lapse seismic conforms to the law of "tuning effect" and seismic reflection of "thin bed", and the remaining oil prediction method based on absolute difference of time-lapse seismic is only applicable to the reservoirs with uniform thickness smaller than the tuning thickness or with thickness greater than the tuning thickness. The relative difference of time-lapse seismic is not affected by reservoir thickness, but only related to the change of fluid saturation. It is applicable to all the deep-sea unconsolidated sandstone reservoirs which can exclude the effect of pressure, temperature, pore type and porosity on seismic. Therefore, the relation between the relative difference of time-lapse seismic and the change of fluid saturation, which is obtained from seismic forward modeling after Gassmann fluid displacement simulation, can be used to quantitatively predict the change of reservoir water saturation and then the distribution of the remaining oil. The application of this method in deep sea Zeta oil field in west Africa shows that it is reasonable and effective.

Fracturing and Displacement Types of Active Faults in Unconsolidated Accumulative Materials and Their Implication for Timing Palaeoseismic Rupture Events

Seismic rupture produced by active faulting causes ground surface fracturing and seriously damages buildings.However,the fracture traces are usually indistinct or non-visible because of complicated deformations in Quaternary unconsolidated materials,so their upper parts can hardly be discovered on exploratory trenches.The aim of this paper is to study the types of fault rupture and displacement in different loose deposit beds from a joint view of mega-,meso-,and microscopic scales,and to find the mark and method for timing the palaeoseismic rupture on exploratory trenches.

Productivity Testing Design Method of Multi-Factor Control for Unconsolidated Sandstone Gas Reservoir

Reservoir safety, testing-string safety, and flow control are key factors that should be considered in deep-water unconsolidated sandstone gas well testing work system. Combined with the feature of testing reservoir, pipe string type and sea area, the required minimum testing flow rate during cleaning up process, as well as minimum test flow rate without hydrate generation, pipe string erosion critical production, the maximum testing flow rate without destroying sand formation and the minimum output of meeting the demand of development was analyzed;based on the above critical test flow rates, testing working system is designed. Field application showed that the designed work system effectively provided good guidance for field test operations;no sand production or hydrate generation happened during the test process;the test parameter evaluated the reservoir accurately;the safe and efficient test operation was achieved.

A new method to optimize the fracture geometry of a frac-packed well in unconsolidated sandstone heavy oil reservoirs

The worldwide proven recoverable reserves of conventional oil are less than the amount of the heavy oil.Owing to weakly consolidated formation,sand production is an important problem encountered during oil production in heavy oil reservoirs,for which frac-pack technique is one of the most common treatments.Hence,how to obtain the optimal fracture geometry is the key to increasing well production and preventing sand.Due to the faultiness that current optimization of the fracture geometry only depends on well productivity,fracture-flow fraction was used to describe the contribution of the fracture collecting and conducting fluids from the reservoir.The higher the fracture-flow fraction,the more likely bilinear flow pattern occurs,thus leading to smaller flow resistance and better results in oil productivity and sand prevention.A reservoir numerical simulation model was established to simulate the long-term production dynamic of a fractured well in rectangular drainage areas.In order to reach the aim of increasing productivity meanwhile preventing sand,a new method based on Unified Fracture Design was developed to optimize the fracture geometry.For a specific reservoir and a certain amount of proppant injected to the target layer,there exits an optimal dimensionless fracture conductivity which corresponds to the maximum fracture-flow fraction,accordingly we can get the optimal fracture geometry.The formulas of the optimal fracture geometry were presented on square drainage area conditions,which are very convenient to apply.Equivalent Proppant Number was used to eliminate the impact of aspect ratios of rectangular drainage area,then,the same method to optimize the fracture geometry as mentioned for square drainage areas could be adopted too.

Double-Frequency Microseisms on the Thick Unconsolidated Sediments in Eastern and Southeastern Coasts of United States: Sources and Applications on Seismic Site Effect Evaluation

This study presents a systematic analysis of double-frequency(DF) microseisms recorded on the unconsolidated sediments in the eastern and southeastern coasts of United States. For all recordings, the site effect parameters(predominant frequency(f_(0)), amplification factor and unconsolidated sediment thickness(UST)) are obtained by Nakamura method and the DF spectra are classified into five groups in terms of the DF peak patterns and the recording locations relative to the coastline. The frequencies and energy levels of the DF peaks in horizontal direction and the amplification factors are associated with the UST which is resulted from seismic site effect. By polarization analysis, the primary vibration directions of the DF peaks are identified and presented as great circles passing through the recording stations intersecting mainly along the continental slope. Correlation analyses of time histories of the DF energy and the ocean wave climate observed at buoys show that the low(<0.2 Hz) and high(>0.2 Hz) frequency DF microseisms are generated in the deep ocean and the continental shelf respectively. It is concluded that the continental slope plays a significant role in the generation of DF microseisms as it causes reflection of waves from the open ocean, initiating standing waves.

Fluid flow with compaction and sand production in unconsolidated sandstone reservoir

The fluid flow of unconsolidated sandstone reservoir can be affected by compaction and sand production which will damage the reservoir and affect oil well productivity.This study aims to measure how the two factors affect the fluid flow.Firstly,single-phase displacement test was applied to investigate how the permeability changed with compaction.Then two-phase displacement test assessed the influence of compaction on oil production.Finally,the characteristics of fluid flow with compaction and sand production were studied under different water content.The results demonstrate that the reduction of permeability with compaction is irreversible,which will result in lower productivity.In contrast,sand production can increase the permeability at mid and high water content,which slows down the decline of oil production.Generally,the oil well productivity is reduced because of compaction even with sand production,especially when the formation pressure drop varies from 2MPa to 4MPa.Consequently,advance water injection is necessary to keep the formation pressure and oil production during oilfield development of unconsolidated sandstone reservoir.Simultaneously,the study can provide theoretical basis and references for the similar reservoirs.

Roof pre-blasting to prevent support crushing and water inrush accidents

Support crushing and water inrush when mining under an unconsolidated confined aquifer in the Qidong Coal Mine was prevented by roof pre-blasting. The mechanism and applicable conditions for this method have been studied. The results show that when an overburden structure that may cause support crushing and a water inrush accident exists the weakening of the primary key stratum, which thereby reduces its weighting step, roof pre-blasting is both feasible and effective. If the position of the primary key stratum can be moved upward to exceed 10 times the mining height the possibility of support crushing and water inrush disaster caused by key stratum compound breakage will be lowered. The overburden structure of the number 7121 working face was considered during the design of a technical proposal involving roof pre-blasting. After comprehensively analyzing the applicability of roof pre-blasting the resulting design prevented support crushing and water inrush disasters from happening at the number 7121 working face and laid a solid foundation for mining safely.

Porothermoelastic response of an oil sand formation subjected to injection and micro-fracturing in horizontal wells

Stimulation of unconsolidated formations via horizontal wells has seen its vast implementation in the recent development of heavy oil reservoir to save the time and cost of preheating the reservoir before the steam-assisted gravity drainage(SAGD)process.A mathematical approach was proposed in this research that fully couples the hydraulic,mechanical and thermal responses of unconsolidated sandstone formations and also applies failure criteria for describing either shear dilation or tensile parting mechanism that generates microcracks.The approach was implemented to predict the porothermoelastic response of a pair of SAGD wells subject to injection and subsequent micro-fracturing using hot water.It was found that the predicted bottom hole pressures(BHPs)match closely with the field observed data.An elliptical dilation zone developed around the dual wells with relatively high pore pressure,porosity,permeability and temperature,implying good interwell hydraulic communication between both wells.The activation of microcracks dramatically accelerated the dissipation of pore pressure across the entire formation depth and also facilitated heat convection in between the dual wells,though to a lesser extent.In summary,the approach provides a convenient means to assist field engineers in the optimization of injection efficiency and evaluation of interference among multiple horizontal wells.

Factors influencing physical property evolution in sandstone mechanical compaction:the evidence from diagenetic simulation experiments

In order to analyze the factors influencing sandstone mechanical compaction and its physical property evolution during compaction processes, simulation exper- iments on sandstone mechanical compaction were carried out with a self-designed diagenetic simulation system. The experimental materials were modem sediments from dif- ferent sources, and the experiments were conducted under high temperature and high pressure. Results of the exper- iments show a binary function relation between primary porosity and mean size as well as sorting. With increasing overburden pressure during mechanical compaction, the evolution of porosity and permeability can be divided into rapid compaction at an early stage and slow compaction at a late stage, and the dividing pressure value of the two stages is about 12 MPa and the corresponding depth is about 600 m. In the slow compaction stage, there is a good exponential relationship between porosity and overburden pressure, while a good power function relationship exists between permeability and overburden pressure. There is also a good exponential relationship between porosity and permeability. The influence of particle size on sandstone mechanical compaction is mainly reflected in the slowcompaction stage, and the influence of sorting is mainly reflected in the rapid compaction stage. Abnormally high pressure effectively inhibits sandstone mechanical com- paction, and its control on sandstone mechanical com- paction is stronger than that of particle size and sorting. The influence of burial time on sandstone mechanical compaction is mainly in the slow compaction stage, and the porosity reduction caused by compaction is mainly con- trolled by average particle size.

Hydrochemistry of the Natural Low p H Groundwater in the Coastal Aquifers near Beihai, China

Natural weak acidic groundwater occurs in the unconfined and confined aquifers consisting of Quaternary and Neogene unconsolidated sediments near Beihai in southern Guangxi, China. Under natural conditions the groundwater has low TDS(less than 200 mg L-1) and low concentrations of trace elements(less than 100 μg L-1) with a deceasing tend in contents of the Lanthanides(rare earth elements, less than 1 μg L-1) towards higher atomic number. The groundwater ranges in p H from 3.33 to 7.0 with an average value of 5.12(even lower than that of local rainwater, 5.88). p H values in the groundwater are a bit higher in rainy seasons than those in dry seasons and do not show significant increasing or decreasing trend with time. The average p H value in groundwater in the confined aquifers is even a bit lower than that in the unconfined aquifer. Comprehensive analyses of the groundwater environment suggest that H+ in the groundwater may be derived from dissociation of H2CO3, release of the absorbed H3O+ in clay layers and the acidity of rainwater. The H2CO3 in the groundwater may be formed by dissolution of CO2(g). Minerals in the unconsolidated sediment are predominated by quartz with small amount of clay minerals. The sediments undergoing a long-term weathering contain low levels of soluble constitutes. Lack of alkaline substances in the groundwater system is also helpful in the accumulation of acidity of the groundwater.

Effects of Moisture Content and Plasticity Index on Duncan-Chang Model Parameters of Hydraulic Fill Soft Soil

In order to explore the effects of moisture content and plasticity index on Duncan-Chang model parameters?K,n,?C?and?Rf,?we selected 8 groups of soft soil with water content of 69.1%?-?94.3% and plasticity index of 32.2?-?54.1 for triaxial unconsolidated undrained shear test. The results show that?Cuu,?K?and?n?values all showed a downward trend, and?Rf?variation was not obvious with the increase of moisture content. The variation rule of each parameter is not obvious with the increase of plasticity index. When moisture content is constant,?Cuu?and?n?values do not change much,?K?increases with the increase of plasticity index within the range of 70%?-?80% moisture content, and does not change much with the increase of plasticity index when moisture content is greater than 80%,?Rf?has no obvious rule.?When the plasticity index is constant,?Cuu,?Kand?n?decrease with the increase of moisture content,?Rf?has no obvious rule. The maximum value of?Cuu?is 20.18?kPa, the minimum is 3.72?kPa, and the maximum to minimum ratio is 5.42. The maximum value of?K?is 0.517, the minimum is 0.022, and the maximum to minimum ratio is 23.5. The maximum value of?n?is 1.198, the minimum is 0.150, and the maximum to minimum ratio is 7.99. The maximum value of?Rf?is 0.872, the minimum is 0.679, and the maximum to minimum ratio is 1.28.

Electrical Resistivity Imaging of Suspected Seepage Channels in an Earthen Dam in Zaria, North-Western Nigeria

To determine and map the subsurface conditions of a dam, a 2D electrical resistivity tomography study was carried out within the two flanks of Zaria dam at Shika. This was done to ascertain if the variations in the volume of water content in the dam is due to an anomalous seepage beneath the subsurface or seasonal effects. On the basis of the interpretation of the acquired data, various zones of relatively uniform resistivity values were mapped and identified. The first zone is characterized by moderate resistivity values of 150 - 600 ohm-m. It represents unsaturated topsoil with thicknesses varying from 1 - 4.5 m. The second (intermediate depth) resistivity zone, with values ranging from 5 - 100 ohm-m and thickness varying from 3.5 - 10 m, represents a silt clay layer with high moisture content. The third resistivity zone represents fairly weathered granite and is characterized by relatively high resistivity values ranging from 700 - 6000 ohm-m. The available borehole log data correlated well with the pseudo-sections in relation to the obtained resistivity values and depth. Zones of relatively low resistivity within the bedrock are interpreted to represent potential seepage pathways. Hence, this geophysical method can be successfully used to delineate and map these seepage pathways within the subsurface of the earth dam.

A two-step multi-frequency receiver function inversion method for shallow crustal S-wave velocity structure and its application across the basin-mountain range belts in Northeast China

A shallow crustal velocity structure(above 10 km depth) is essential for understanding the crustal structures and deformation and assessing the exploration prospect of natural resources, and also provides priori information for imaging deeper crustal and mantle structure. Passive-source seismic methods are cost-effective and advantageous for regional-scale imaging of shallow crustal structures compared to active-source methods. Among these passive methods, techniques utilizing receiver function waveforms and/or body-wave amplitude ratios have recently gained prominence due to their relatively high spatial resolution. However, in basin regions, reverberations caused by near-surface unconsolidated sedimentary layers often introduce strong non-uniqueness and uncertainty, limiting the applicability of such methods. To address these challenges, we propose a two-step inversion method that uses multi-frequency P-RF waveforms and P-RF horizontal-to-vertical amplitude ratios. Synthetic tests indicate that our two-step inversion method can mitigate the non-uniqueness of the inversion and enhance the stability of the results. Applying this method to teleseismic data from a linear seismic array across the sedimentary basins in Northeast China, we obtain a high-resolution image of the shallow crustal S-wave velocity structure along the array. Our results reveal significant differences between the basins and mountains. The identification of low-velocity anomalies(<2.8 km s^(-1)) at depths less than 1.0 km beneath the Erlian Basin and less than 2.5 km beneath the Songliao Basin suggests the existence of sedimentary layers. Moreover, the high-velocity anomalies(~3.4–3.8 km s^(-1)) occurring at depths greater than 7 km in the Songliao Basin may reflect mafic intrusions emplaced during the Early Cretaceous. Velocity anomaly distribution in our imaging result is consistent with the location of the major faults, uplifts, and sedimentary depressions, as well as active-source seismic results. This application further validates the effectiveness of our method in constraining the depth-dependent characteristics of the S-wave velocity in basins with unconsolidated sedimentary cover.