Online First

Articles in press have been peer-reviewed and accepted, which are not yet assigned to volumes/issues, but are citable by Digital Object Identifier (DOI).
Display Method:
Experimental study on static and dynamic water dissolution characteristics of red layer gypsum interlayer
ZHONG Zhibin, FENG Jie, LV Lei, ZHOU Qijian, LI Sijia, XUE Changrui
, Available online  , doi: 10.11932/karst2024y046
Abstract:
Gypsum interlayers with uneven thickness are widely distributed in the red layer sand and mudstone in Sichuan Basin. In the process of groundwater infiltration or erosion, these gypsum interlayers will inevitably produce dissolution reaction and form a large number of holes. Coupled with the special physical and mechanical properties of red layer, it is easy to reduce the bearing capacity of foundation, uneven settlement deformation and even damage of engineering structures. Therefore, it is necessary to carry out more in-depth research on the dissolution characteristics of gypsum interlayer in red layer.Based on an avionics hub project in Chongqing, this paper takes the gypsum sandwich rock at dam foundation as the research object, and uses a self-designed dynamic water dissolution device to carry out the dissolution test. The dissolution process of gypsum in the test can be characterized by regularly monitoring the concentration changes of $ {\mathrm{C}\mathrm{a}}^{2+}\mathrm{、}{{\mathrm{S}\mathrm{O}}_{4}^{2-}} $ ions in the solution. The dissolution characteristics of gypsum interlayer in red layer at different flow rates were investigated.At the end of the test, different degree of dissolution grooves appeared in the gypsum interlayers under different flow conditions. Among them, the grooves on the surface of the sample with static water corrosion are the smallest, only 2~3mm, while the grooves produced by sample corrosion in the moving water environment reach the deepest 12mm. It can be found that the dissolution rate of gypsum in the moving water environment has been greatly improved. The concentrations of $ {\mathrm{C}\mathrm{a}}^{2+}\mathrm{、}{{\mathrm{S}\mathrm{O}}_{4}^{2-}} $ in the solution were measured after the test, and it was found that the contents of both ions showed an increasing trend with time. Comparing the results of hydrostatic and dynamic water dissolution tests, it can be seen that the slope of $ {\mathrm{C}\mathrm{a}}^{2+} $ content change curve of experimental groups increases gradually with the increase of flow velocity, indicating that higher flow velocity has a certain promotion effect on the chemical dissolution of gypsum. Through the processing and conversion of the test data, it can be found that the average dissolution rate under different flow velocity environment shows a rapid decline in a short time, and then gradually tends to be stable. At the initial time of the corrosion test, the average chemical corrosion rate of the samples under different flow rates showed that the faster the flow rate, the greater the corrosion rate. With the experiment, the dissolution rate decreases rapidly and the decreasing trend in different flow velocity environments is different to some extent. The greater the flow velocity is, the greater the dissolution rate is. Compared with the dynamic water corrosion test, the average chemical corrosion rate of the samples in the static water environment is smaller in both value and range. Combined with the change curve of the average chemical dissolution rate of each sample, it can be found that the average dissolution rate of gypsum increases with the acceleration of the flow rate, which proves that the larger flow rate plays a role in promoting the dissolution of gypsum interlayer. The trend of instantaneous chemical dissolution rate curve of each dynamic water dissolution test group was basically consistent with the trend of average chemical dissolution rate curve, showing a law of decreasing with time, and all curves showed a sudden drop on the third day of the test, and then gradually decreased and tended to be stable. According to the curve of the dissolution amount of each sample with the flow rate, it can be found that the chemical dissolution and mechanical submergence of gypsum interlayer will be strengthened when the flow rate increases.
Geological significance, prediction, and application of structural fractures in tight limestone reservoirs of the Taiyuan Formation in Ordos Basin
FU Xunxun, ZHAO Weibo, PEI Wenchao, ZHAO Xiaohui, SHI Lei, DONG Guodong, WANG Huaichang, LIU Xiaopeng, LU Zixin, YIN Liangliang, GAO Xing, CHEN Juanping
, Available online  , doi: 10.11932/karst2024y040
Abstract:
In Ordos Basin,the Taiyuan Formation mainly consists of a marine-terrestrial transitional sedimentary system.The northern part featured shallow water delta deposits dominated by terrigenous clastic rocks, while the southern part consists of coastal marine limestone deposits. The Shenmu Gas Field has been discovered in the northern terrigenous clastic rock strata. The southern limestone strata has a wide distribution, covering approximately 14 square kilometers, and exhibits considerable thickness, ranging from 5 to 30 meters.However, the limestone reservoir is dense with an average porosity of 2.1%, an average permeability of 0.21mD, and an average throat radius of 0.12μm. For a long time, exploration had not focused on it due to its dense reservoir lithology.Recent years, further research has indicated that the high-yield gas wells are often related to the development degree of structural fractures.However, there is a lack of systematic analysis of whether the quality of the limestone reservoir and the natural gas enrichment relate to these structural development.To solve the above difficulties, the following steps have been studied and adopted in this study.Considering the tectonic setting,the structural fracture characteristics the of the Taiyuan limestone were firstly analyzed by field profiles, core observations, and well log data.Secondly the impact of structural fractures on the porosity of dense limestone reservoirs was observed by uniaxial loading rock core nuclear magnetic resonance (NMR) technology and CT scanning for the first time. Furthermore, the fluid inclusion was investigated to analyze the conducing role of the fracture networks during gas accumulation process.Lastly, based on the work mentioned above,to predict the development of structural fractures and favorable blocks exploration, an improved Gaussian curvature method (IGC) was proposed,according to the specific geological conditions, and it was applied in well deployment in the Hengshan area.The results show that: 1) In the study area, regional tectonic stress field is mainly influenced by the Yanshanian Movement,and the main tectonic framework was mainly established by it with NW-SE compressive stress.Compared with low-angle fractures and oblique fractures, high-angle tension fractures have low filling levels and good effectiveness.it often penetrates, cuts, and dislocates other types of fracture improving the limestone reservoir properties. In addition, the formed fracture network system may facilitate the efficient migration of hydrocarbons; 2) Rock physics experiments had manifested that structural fractures greatly contribute to the physical properties of the tight limestone reservoirs.With the development of fractures, the average porosity doubled from 2.1% to 4.2%. The contribution rate of fractures to total porosity ranged from 14.3%~72.7%, with an average of 43.4%. Among the three reservoir rock types, the porosity of powder crystal limestone increased the most, followed by algae bound limestone, and mud crystal limestone had the smallest increase; Combined with typical well reservoir-forming temperature, burial history, and thermal evolution history in the study area,inclusion analysis made clear that the widely developed fracture networks formed an effective transport system, and hydrocarbons effectively migrated by it during the main reservoir-forming period from the Late Jurassic to the Early Cretaceous.That is to say, generated from the main coal-bearing source rocks 8# coal and 5# coal located at the top and bottom of the limestone reservoir,the natural gas entered the limestone reservoir by the effective transport system.And supplemented by a cap rock of the marine mudstone on the top of the Taiyuan Formation, a superior "sandwich" type reservoir formed;4) The improved Gaussian curvature method (IGC) can predict the development degree of favorable structural fractures. The higher the IGC value, the more developed the high-angle tension fractures, and the more favorable the conditions for natural gas enrichment.This method was a debut in the Hongshan area to support well deployment. YT1H*,the first risk exploration well, achieved a gas production of 540,000 cubic meters per day, and marked a breakthrough in the exploration of the Taiyuan limestone in the Ordos Basin. In the past 2 years, 9 industrial gas wells have been discovered in this area, 5 of which have output above 100,000 cubic meters per day. Among them, YT4H* achieved an output of 1,088,000 cubic meters per day.The above exploration achievement demonstrated the effectiveness of the research and provided a reference for re reacquainting tight carbonate reservoirs. It is expected that this study can further guide the large-scale exploration and development of the Taiyuan Formation limestone in the Ordos Basin.
Study on ecological risk threshold of farmland soil in typical karst region of Guangxi: A case study of Cd
Huang Chunyang, Huang Juan, Zhong Xiaoyu, Li Jie, Li Mingyu
, Available online  , doi: 10.11932/karst2024y032
Abstract:
Cadmium (Cd) is a non-essential element that is harmful to human health, and it is a primary pollutant in agricultural soils. It is important to assess an ecological risk threshold of Cd in soil-crop systems. The soils have previously been shown to be anomalously enriched in Cd through geogenic processes in the karst region,Southwestern China. Therefore, a valuable method that could accurately evaluate ecological risk threshold of Cd in agricultural soils in karst areas needs to be proposed immediately. The rice is the most widely cultivated cereal crop that is responsible for the largest human dietary exposure to Cd in Daxin county.In this study, we systematically collected 102 groups of crops of rice and its root soil from a typical karst area in Guangxi Province. The results showed that 90.2% of Cd in soil samples exceeded the screening value in the National Environmental Quality Standard for Soil (GB15618-2018), the Cd content in rice seeds planted in those fields, however, were quite low, and only 8.82% samples exceeded the limit for Contaminants in Food (GB 2762-2022). The results revealed a spatial mismatch in Cd pollution between the soil and rice grains in the study area, and the main results are as follows: (i) slight soil pollution and safe rice; (ii) slight soil pollution and serious rice pollution and (iii) safe soil and serious rice pollution. In the study area, the higher the pH value of rice seed growth soil, the weaker the migration ability of Cd, and the less Cd absorbed or transferred by rice seed. Although the soil Cd form is mainly dominated by the content of bioavailable state, the soil pH is mainly weakly alkaline, which inhibits the absorption of soil Cd by plants, and will not affect food security if the soil pH does not change.It is one of the key reasons for the low Cd concentration in rice with high Cd soil geological background. The current soil environmental quality standard is difficult to accurately evaluate the level of soil Cd pollution in the study area, so the soil environmental quality benchmark should be proposed according to the actual situation.Therefore, based on the species sensitivity distribution method, this study calculated the safety limit for protecting 95% of rice seed soil Cd, and evaluated the actual Cd pollution level in the study area.Based on SSD method, we deduced the pollution risk screening values as 0.4 mg·kg-1, 1.08 mg·kg-1 and 6.4 mg·kg-1 corresponding to soil Cd that protected 95% rice from exceeding the standard under pH 5.5~6.5, 6.5~7.5 and 7.5~8.5, respectively.Compared with the soil environmental quality limits stipulated by GB15618-2018, it can be found that the national standard is not rigorous enough and too broad in the environment of acidic soil (5.57.5) soil environment.With the increase of pH value, the correct rate of derived value increased from 23.5% to 92.1%, while the correct rate of screening value given by national standard decreased from 47.1% to 6.3%. The correct rate of derived value according to the evaluation results increased significantly with the increase of pH, from 23.5% to 92.1%. Therefore, this result could be considered a more accurate assessment threshold, and provide a scientific reference for the safe production of Cd-contaminated rice paddies in typical karst areas of Guangxi.
Discrete fracture network seepage analysis based on MODFLOW
WANG Jinli, CHEN Xi, ZHANG Zhicai, KANG Jianrong, HU Jinshan
, Available online  , doi: 10.11932/karst2024y043
Abstract:
Karst bedrock fractures are difficult to model in hydrological research due to their heterogeneous and anisotropic. Groundwater can mainly flow in the fractures and conduits of the karst aquifer system, and the bedrock voids and micro-fractures mainly play the role of water storage, so the ideal hydrological model of the karst area is the continuum-fracture (conduit) model. As the main path of groundwater flow, the connectivity of fractures has a great influence on the seepage of fracture media. The fracture-continuum model (FC) generated by combining the merits of the continuum model and the discrete fracture network model (DFN) can realistically describe the binary structure of karst aquifer system in principle, and the diffusion of matrix as well as the water exchange between matrix and fracture can also be considered, but there is a problem that the water exchange between matrix and fracture medium is difficult to determine.This study is based on the statistical distribution of fracture geometry and hydraulic parameters, the two-dimensional fracture network original maps were generated using Monte Carlo random simulation technology and MATLAB program with the same distribution as actual rock fractures. The adjacency matrix of undirected graph of graph theory is used to characterize the intersecting relationship between fractures in the fracture network. Using percolation theory to remove fracture segment which cannot conduct fluid, such as the 1st connected maps were generated by removing the isolated fractures (fracture clusters), and the 2nd connected maps generated by removing fractures with single nodes, as well as the 3rd connected maps generated by means of cutting fracture dead ends. The DFN model in a single-phase, saturated, steady-state flow system was established based on the connected graph of fracture network, applying the cubic law and the continuity equation of seepage. Mapping of single fractures of the connected graph onto finite difference grids, FC seepage model was built by combining DFN approach and MODFLOW software. When deleting the fracture dead ends, some nodes may be lost due to computer trade-off, as a result, the 2nd connected graph and the 3rd connected graph are used for DFN and FC approaches, respectively. The corresponding simulation programs were developed, with which reproducing cross fracture test conditions, and constructing two multi-fractures cases to analyze the impact of finite difference grid resolutions and fracture dead end on stable seepage simulation accuracy. The DFN flow solution is used as a metric to evaluate the correctness and effectiveness of FC model code writing.Research results show that the simulated total flow and head of the fine grid are in good agreement with the theoretical values of DFN under the two test conditions. The node head and total grid flow simulated by the FC model of two cases are overestimated or underestimated in the coarse grid. Although a longer path of grid flow Lg from cell to cell is corrected by increasing cell hydraulic conductivity, the coarse grid based on MODFLOW cannot accurately analyze the inclined fracture length Lf on the grid. If the selected coarse grid size is greater than the distance between two adjacent fractures, a part of the two adjacent fractures will merge in the grid, thereby shortening the actual flow path. Secondly, the overestimation or underestimation of the inclined fracture length Lf.. Inaccurate analysis of the flow path will lead to an increase in the estimation error of flow rate and head, especially in the densely fractured areas near the center of the model domain. Refining the grid can greatly eliminate this error. The fracture dead-end has non-negligible effect on MODFLOW grid flow. In this study, the MODFLOW-based FC model can be used to solve for pressure distributions both within the interconnected network of fractures and the rock matrix, and could also be a solution for the water exchange between the rock matrix and the fractures is difficult to determine. This study can achieve the expression of heterogeneity and anisotropy of fracture hydraulic conductivities in numerical calculation units and will promote the understanding of MODFLOW simulation of fracture flow based on connectivity.Furthermore, The size of the MODFLOW grid should be determined in combination with the size of the study area, which should be small enough to capture the detailed features in the flow, and large enough to reduce the computational time and cost. Generally, the grid size is determined by a trial algorithm. When FC approach is applied to the actual watershed, the spatial variability of topography and landform factors should also be considered, such as coupling surface elevation DEM data.
, Available online  
Abstract: