• Included in CSCD
  • Chinese Core Journals
  • Included in WJCI Report
  • Included in Scopus, CA, DOAJ, EBSCO, JST
  • The Key Magazine of China Technology

2022 Vol. 41, No. 5

Display Method:
Exploration of superficial soil-rock structure for karst area based on frequency domain electromagnetic method
CHENG Ping, CHENG Qinbo, CHEN Xi, LIU Jintao, ZHANG Zhicai, GAO Man
2022, 41(5): 675-683. doi: 10.11932/karst20220501
Abstract:
Soil thickness as a key hydrological and ecological factor, is distributed extremely uneven in karst areas. The explorations of the soil thickness and soil-rock interfaces are still challenge. In this study, a 1-D electrical conductivity (EC) inversion model was developed for the frequency domain electromagnetic method (EMI) based on Maxwell's equation system. The visualization of EC distribution for soil profile in karst area could be realized by the model. The inversion model was validated according to the given ideal stratigraphic EC data and applied to two detection lines and the three exposed profiles at a karst depression. The detected apparent electrical conductivities and soil-rock interfaces from EMI are further compared with those from high-density electrical method. The results showed that the inversion from EMI can capture the variation of ECat ideal strata and interfaces of soil layer and underlying limestone and dolomite bedrock at field. The distribution of soil thickness could be estimated with the detected soil-rock interface. However, there is relatively low accuracy when the model is applied to detect trench (trough) at the small scale or the soil-rock interfaces with mudstone.
Application of self-potential and high-density resistivity method to the water exploration in karst terrain of middle-lower reaches of Xijiang River
CHEN Yixiang, HUANG Qibo, QIN Xiaoqun, HAN Kai, XIAO Qiong, MIAO Ying, DU Chengliang, HE Dehuang
2022, 41(5): 684-697. doi: 10.11932/karst20220502
Abstract:
Based on hydrodynamic conditions of the shallow karst groundwater seepage, we mainly study the natural polarization phenomenon of geological structures in water-enriched karst induced by the action of groundwater seepage and distribution characteristics of the natural electric field therefrom. By comprehensive prospecting with the combination of self-potential method and high-density electronic resistivity method, we analyse and evaluate the application effect of the combination on the exploration of karst groundwater in different hydrogeological environments of middle and lower reaches of Xijiang River. In this study, we also integrate geophysical prospecting results with the data derived from exploration & production wells in the survey project of hydrogeological environment initiated and supported by China Geological Survey. Study results show that induced electrical methods, such as transient electromagnetic method, audio magnetotelluric sounding method, ground nuclear magnetic resonance method, induced polarization method, sonic frequency geoelectric field method, VLF electromagnetic method, etc., are seriously affected by electromagnetic fields like power transmission and communication network. The high-density resistivity method shows a strong anti-interference ability and high signal-to-noise ratio. But its exploration depth is relatively shallow, and it is easily disturbed by ore bodies with high conductivity such as charry limestones or karst minerals. The self-potential method is used in the situation of polarization of geological structures in water-enriched karst caused by groundwater migration. Anomalies of the natural electric field are closely related to the groundwater recharge and migration. The complementation of self-potential method and high-density electronic resistivity method in groundwater exploration can contribute to the mutual corroboration of prospecting results, hence improving prospecting resolution and precision.
Study on the microstructure and variation law of carbonate rock dissolution based on scanning electron microscopy and CT imaging technology
DENG Yulin, GUO Xulei, LUO Mingming, CHEN Xiangyong, KUANG Ye, ZHOU Hong
2022, 41(5): 698-707. doi: 10.11932/karst20220503
Abstract:
The dissolution of carbonate rocks is the most important process of karst development, and it is also the basis for the processes affecting karst geomorphological evolution, hydrochemistry of karst water, karst water system development, karst carbon sinks and global carbon cycle. Scholars have explored the dissolution mechanism of carbonate rocks and the factors that affect the development of carbonate karst through indoor simulated dissolution and in-situ experiments, but they mostly focus on the effect of chemical components of rocks on karst processes. In field surveys, it was found that not only limestone formations, but also some dolomite formations have developed karst. It is not comprehensive enough to explain the solubility of carbonate rocks only by chemical composition. The pore structure of the rock also plays an important role,but the influence mechanism of rock pore structure on karst development is still unclear. In order to study the variation law and control effect of carbonate rock pores under dissolution experiment,in this paper, we selects four types of carbonate rocks (Sinian Dengying Formation-medium fine-grained dolomite, Cambrian Tianheban Formation-oolitic limestone, fine-grained limestone, and Qinjiamiao Formation-dolomitic breccia)in the Three Gorges area to carry out dissolution experiments. At the same time, scanning electron microscopy and CT imaging methods were used to test and analyze the dissolution characteristics and pore structure of rocks before and after dissolution reaction. The results show that, (1)The dissolution of carbonate rocks always occurs in soluble minerals with low lattice energy and develops along the rhombohedral cleavage planes and weak parts of mineral crystals,showing the selective dissolution of soluble minerals and pore structures. Dissolution preferentially dissolves soluble minerals, and insoluble components either fall off or protrude from the rock surface.Due to the control of rhombic cleavage on the surface of dolomite crystals, the dissolution will develop into dissolution fractures along the cleavage, and at the same time, dissolution pores appear on the surface of dolomite crystals, showing a "honeycomb" dissolution characteristic; Lattice defects, cleavage, and the weaker plasmonic and molecular contact forces at the edge of the crystal are the places where the dissolution preferentially occurs, and the dissolution strength at the rock pores is slightly higher than that at other parts. (2) The pore size of carbonate rocks controls the development of dissolved pores and the evolution of dissolved pores during rock dissolution.For oolitic limestone and fine-grained limestone samples, because their pore diameters are mainly concentrated in 18-80 μm and the pore distribution is uniform,after the dissolution reaction, the small pores further increased, only a small amount of pores were enlarged by dissolution, and the area with increased pore size was mainly concentrated in 40-50 μm,and the pores were mostly developed on the surface of the rock. The pores on the surface of the medium-fine-grained dolomite are relatively developed and the pore size is large, about 20% of which are larger than 100 μm, after the dissolution transformation, the proportion of pores with a pore size of more than 50 μm in the rock increased to 78%. It is shown that the dissolution of rocks with small pore size (40-60 μm) mainly occurs on the surface of the sample, and small dissolved pores are often developed. The dissolution of large pore size (>80 μm) rocks mainly occurs in the pore walls and there are signs of dissolution into the rock. After dissolution transformation, the pore throat radius and connectivity of medium-fine-grained dolomite show a certain increase. (3) The dissolution process of carbonate rocks is jointly affected by chemical composition and pore structure, and is less affected by rock porosity. The pore size distribution is the key factor affecting the dissolution rate of carbonate rocks. The porosity of the samples before dissolution is between 1.97 and 10%. Under the experimental conditions (20 ℃, 1 atm, 0.1 mol·L−1 HCl solution), the dissolution rate of carbonate rocks is between 0.497 g·d−1 and 0.598 g·d−1. The mass loss of dolomitic breccia through dissolution is the largest due to the influence of cement shedding; the medium-fine crystal dolomite with more developed surface pores and larger pore size is second;fine-grained limestone with dense grains and small pores has the slowest dissolution rate.Samples with smaller pore diameters are greatly affected by the stagnant boundary layer, and the dissolution medium is likely to form a local saturation state inside the pores; the samples with larger pore diameters can expand into the rock during the dissolution process of the "water-rock" reaction, which is affected by the stagnant boundary, but the effect of the layer is less and the water flow conditions are better than the small pore size samples, so the dissolution rate is higher. This study has certain guiding significance for the understanding of the differential karstification mechanism and karst development law of carbonate rocks.
Response characteristics of forward modeling of 3D high-density resistivitymethod on different devices in the fault-water-filled cave
Yang Yanfang, Ju Hejian, Gan Fuping, Cheng Yang, Wang Yong
2022, 41(5): 708-717. doi: 10.11932/karst20220505
Abstract:
Low blocking layers and water-filled caves in high resistivity carbonate rocks are the main target bodies for groundwater exploration in karst areas. Three-dimensional imaging by high-density resistivity can obtain rich and high-precision geological anomaly information in any direction in the whole underground space of the survey area through multi-line and multi-direction observation, which is one of the most commonly used methods for water exploration of carbonate fissure caves. According to the principle of finding low resistance in high resistance, the water occurrence model of carbonate fissure karst cave in overlying karst areas is simplified into a combined geoelectric model of fault-karst cave, and the physical property parameters of various bodies are assigned according to the resistivity characteristics of common media in the karst area of south China. The specific parameters of the model are as follows: the surface overlying layer is cultivated soil with a thickness of 3 m and a resistivity of 100 Ω·m. In the fractured limestone with an underlying resistivity of 4,000 Ω·m, a tensile water-rich fault, with a dip angle of 78° and a width of 10m, is developed with a resistivity of 800 Ω·m. Water-filled cave 1, with a diameter of 15 m and a central burial depth of 17.5 m, is developed in the hanging wall of the fault, and Water-filled cave 2, with a diameter of 10m and a central burial depth of 15 m, is developed in the footwall. The resistivity of both caves is 25 Ω·m. Given the burial depth, resolution, working efficiency and other factors of the target body, the mesh of the 3D model is divided into x×y=50×30 electrodes. The number of layers is 12; the electrode distance is 5m; the measurement range of x direction is 0-245 m; the measurement range of y direction is 0-145 m. RES3D software has been used to perform forward and inverse simulation calculation of Winner, Schlumberger and Dipole-dipole acquisition devices. Inversion results within the range of detection depth are displayed in a three-dimensional manner, and XY horizontal slices are displayed on the model interface (Z=0 m, 3 m, 10 m, 25 m, and the maximum detection depth). To further explore the response characteristics of the model on different devices, XZ profile① across the karst cave and fault area (y=67.5 m), XZ profile② far away from the karst cave (y=7.5 m), and the sounding curves of different points (Ground sounding point I of Karst cave 1, Ground sounding point II of Karst cave 2, Ground sounding point III of Karst cave, and Ground sounding point IV far from the karst cave area) are extracted. The electrical characteristics of the three devices at and away from karst caves, the development depth and boundary of Karst caves 1 and 2, the inclination of faults, the burial depth at the top, the horizontal width, the inclination of the dip and the boundary of the hanging wall are analyzed. Through the three-dimensional forward and inverse calculation of the fault-water-filled karst cave, combining with the typical profile and electrical sounding curve, we analyze the response characteristics and laws of the target body in different devices and draw the following conclusions, (1) Within the detection depth, Winner, Schlumberger, Dipole-dipole devices can effectively identify the fault and the cave of the upper wall three times the size of electrode distance of Cave 1. The target body is a concentric circle of low resistance trap anomaly with green ribbon and color gradient. These three devices cannot distinguish the fault footwall of two times the size of the electrode distance of Cave 2. (2) Under the same model and observation conditions, the Dipole-dipole device has the strongest recognition ability for the target body. Three-dimensional inversion results can identify the lower boundary of the cave, and the left convex low resistance trap anomaly is formed in the profile. The curve types, inflection points and extreme points of different sounding points are closest to the model design. Inflection points correspond to the lithology of interface. And extreme points are located at the central depth of the geological body. Therefore, when using the three-dimensional high-density resistivity method for water exploration in the karst area of south China, we can prioritize Dipole-dipole observation method with the careful consideration of selecting the electrode distance. This study is of guiding significance for the selection of field observation mode and geological interpretation.
Karst hydrogeological problems and countermeasures of a proposed railway in plateau slope area of southeast Yunnan
ZHANG Hua, KANG Xiaobo, WANG Bo, CHAI Jinlong, ZHOU Cuiqiong, CAI Shuangle, HOU Xutao, HUANG Chenhui, PAN Xiaodong
2022, 41(5): 718-727. doi: 10.11932/karst20220506
Abstract:
The study area is located in the plateau slope area of southeast Yunnan at the southern margin of Yunnan-Guizhou Plateau, which is the surface watershed zone between Nanpanjiang River Basin and Honghe River Basin at the source of the Pearl River. There exist intense erosion and cutting, complicated topographic and geomorphic conditions and geological structure. The study area belongs to Gejiu sub-region in the stratigraphic region of South China. The widely exposed middle Triassic Gejiu formation and other aquifer groups are rich in water. Due to the complicated hydrogeological conditions, the study area is divided into 2 first-level water systems, 7 second-level water systems, 26 third-level water systems, and 12 fourth-level water systems. Karst is strongly developed with abundant precipitation and strong hydrodynamic effect. Underground river pipelines, depressions and sinkholes are also widely distributed. There are 9 large springs and 20 underground rivers (underflows). Through field investigation, observation, tracer tests and experiments, the influence of karst hydrogeological conditions and problems on the route selection of a proposed railway in the study area are analyzed and assessed under the guidance of geological and karstological theories. The study shows that 80% of the aquifers are the carbonate aquifer group with uneven water abundance from the intermediate to extremely strong degrees. Groundwater mainly exists in karst pipeline. Karst is intensively developed in the study area where the tertiary dissolution process creates the vertical distribution of multi-layer karst caves. At the depth of 100 m, strong karst caves are developed with at least 1 layer or at most 3-4 layers. Large springs, underground rivers and pipelines are also widely distributed. The length of karst pipelines ranges from 6 km to37.75 km. The maximum 6 m3s−1. tracer test results show that there exist underground dissolved pools, or lakes in the main pipeline of underground rivers. It is predicted that the middle line of the tunnel will suffer the largest water inflow, followed in turns by the south line, the north line and the through-line. The through-line is expected to experience relatively small risk of water inflow in tunnel excavation, and hence it is recommended in route selection. Geothermal resources are distributed near the contact zone between the Bozhushan granite body and the Cambrian strata in the west of Wenshan. The temperature of hot springs generally ranges from 45° to 50°. In the belt-shaped geothermal reservoir, temperatures are 70-110 ℃ by analogy of those in the neighboring area. The gray and thin medium dolomite intercalated with calcareous dolomite in the middle Cambrian Longha formation (Є2l) shows serious desertification. The surface sanding dolomite is broken with poor stability, and can be broken by hand. Consequently, accidents such as surrounding rock collapse and roof falling are quite possibly to occur in tunnel excavation. In this study, suggestions to solve karst geological problems in the study area are put forward. Several areas should be avoided in tunnel excavation, such as the areas with rich groundwater, especially those with strong water abundance, the areas with intense Karst development, especially the vertical development areas of underground and surface karst, the runoff area of the main pipeline of Daquan underground river where the route should be planned as far as possible above the recharge area or water level of the main pipeline, the geothermal anomaly area exposed by granites, and the areas distributed with the karst sandstone stratum. These suggestions are expected to reduce the karst hydrogeological hazards and the damage to groundwater environment that the proposed railway construction may face, and to provide reference for the engineering construction of karst areas with intense development.
Study on karst leakage caused by the combination of fault structure and interlayer corrosion zone of a hydropower station in southwest China
FENG Zhigang, WEI Guojian, ZHANG Hanmeng, XU Guo
2022, 41(5): 728-735. doi: 10.11932/karst2021y33
Abstract:
A hydropower station in southwest China is located in the middle part of the fold belt between the eastern Yunnan platform and the western Guizhou platform on the Yangzi platform. The largest fault in the near field area is F38, which is 61 km long and generally protrudes to the northwest in an arc shape. The dam site is located in the east-west structural belt on the south side of F38 fault, which is the core area of the northeast end of Dashuitang oblique anticline with an axial direction of NNE. At the dam site, the river flow direction is 25°-30°, and the valley is basically in symmetrical V shape. The dam area is dominated by carbonate rocks, and faults and karst are relatively developed. There are karst funnels, depressions, sinkholes and other surface karst forms in the upper level of planation surface on both sides of the dam site. Early investigation results of adits, drilling, foundation pit excavation and geophysical prospecting show that there are some karst caves, trough-like caves and interlayer corrosion in the riverbed and both sides, which are interconnected according to their extension trend. Since the reservoir starts to store water, continuous leakage has occurred in the 1,315 m drainage tunnel of the resistance body on the right bank of the dam site. With a rapid rise of the reservoir water level, a large number of muddy fillings gush out from each water inflow point. Then the water gradually becomes clear, the leakage gradually increases to about 1.9 m3·s−1, and the reservoir cannot store water normally. In order to find out the leakage path of reservoir water and take measures to reduce leakage, karst leakage research is carried out. 15 boreholes are arranged on the right bank of the dam site, with the borehole bottom elevation from 1,133.3 m to 1,249.9 m, and the borehole water pressure tests of 121 section are carried out below the bottom elevation of the dam foundation curtain. Results show that the permeability area of Lugeon value q≥3 Lu below the curtain bottom elevation is distributed in a band, mainly between 1,200 m and 1,260 m, where Lugeon value gradually decreases from q=3-75 Lu to 2-8.82 Lu toward the right bank. Inspection holes are arranged along the grouting curtain line of the grouting drainage tunnel at the elevations of 1,293 m and 1,340 m on the right bank, and the bottom elevation of the holes is about 1,200 m. The water pressure test, borehole television, the data of simple water pressure test and abnormal records during the grouting process are analyzed comprehensively. 8 inspection holes of the 1,293 m grouted drainage tunnel can be grouted normally, and no hydraulic connection with the water gushing channel of the 1,315 m downstream drainage tunnel is found. 6 inspection holes are arranged in the 1,340 m grouting drainage tunnel, and hundreds of grouting holes are arranged at intervals of 10 m to the left and 2 m to the right of borehole inspection-4. There are 2-3 rows of grouting holes in local sections, and the abnormal tunnel sections are densely grouted. The exploration results of 1,340 m grouting drainage tunnel show, (1) There is no large structural fracture zone developed in the left of inspection-4. Grouting is successfully completed in each hole, and no channel connected with water inrush at 1,315 m is found. (2) For the multiple boreholes between inspection-4 and inspection-3, porous water loss occurs when the drilling is at the elevation of 1,245-1,270 m, and the water pressure cannot be added to the predetermined pressure in the test. After the valves are closed at each water outlet of 1,315 m drainage tunnel, the water gushing phenomenon appears in many holes, especially holes near the right side of the inspection-4, which indicates that there is a leakage channel for reservoir water to bypass the curtain line at the inspection-4 position, and there is a strong hydraulic connection with the downstream 1,315 m drainage tunnel. (3) Most boreholes between inspection-3 and inspection-1 are drilled after the valves of water inlets of the 1,315 m drainage tunnel are closed. In the tunnel section of about 20 m from the left of fault F12 to inspection-2, the probability of borehole water inflow is relatively high when the drilling reaches 1,220-1,260 m, where there is a corrosion zone. Especially in the area 8 m from the left of fault F12, the water inflow phenomenon occurs in almost every borehole between the elevation about 1,240-1,260 m, which shows that the section from the left of fault F12 to inspection-2 is the main channel for the reservoir water to leak around the curtain to the 1,315 m drainage tunnel downstream. According to the water inflow conditions of each inspection hole and grouting hole of the 1,340 m grouting drainage tunnel, combined with the geological characteristics of the right bank resistance body, exploration results and rock permeability test results, the karst leakage path of 1,315 m water inflow is judged as follows: reservoir water infiltrates through the fault fracture zone F12, reaches the dissolution zone of the deep rock bedding plane, and then gradually lifts along the dissolution zone at a dip angle of 30° downstream, and bypasses the curtain line at the elevation of 1,220-1,260 m near inspection-4 and from the fault F12 to inspection-2. After encountering the geological weak zones such as vertical solution grooves, solution ditches and solution caves, reservoir water of the local area seeps upward and finally gushes out from the geological weak point of 1,315 m drainage tunnel. The treatment of karst leakage in this project is a dynamic interactive process of investigation and grouting. With the measures of grouting curtain reinforcement, the main leakage channel is blocked, and the reservoir water leakage is effectively controlled. The leakage of the right bank 1,315 m corridor gradually reduces to about 38 L·s−1. There is basically no water dripping on the tunnel wall, and the treatment meets the design requirements.the treatment meets the design requirements.e design requirements.
Feasibility study on dam and reservoir construction in the catchment area of complex karst water system: Taking Pucha Reservoir of Beipan River as an example
HU Daru, ZHENG Kexun, ZHAO Daiyao, CHEN Zhanheng
2022, 41(5): 736-745. doi: 10.11932/karst20220507
Abstract:
The topographic and geological conditions in the basin of Beipan River are complex with high mountains and deep valleys on both banks, strong karst development and deeply buried groundwater.Consequently, water resources are in great shortage in this area. The construction of water conservancy projects can effectively solve the problem of water shortage. However, karst leakage has become a difficult problem restricting the construction. By comprehensively using the methods of geological survey and mapping, drilling and geophysical exploration, hydrogeological test, karst water system analysis and groundwater balance analysis, this paper demonstrates that a leakage to the adjacent valley or along the karst pipeline deep under the dam foundation will not occur in the Pucha reservoir, but the possibility of solution crack leakage is great. The finite element method is used to simulate the solution crack leakage and analyze the engineering impact. Results show that with the increase of the dissolution rate of T1yn1-1 limestone, the anti-sliding stability coefficient of the dam foundation decreases slightly, and there is an inverse correlation between them. The regression equation is y1=−0.081x +2.678. The potential instability mode of the dam foundation is that the T1yn1-2 marl rock mass is sheared at the upstream, and the bottom is sheared along the gently inclined upstream with weak structural plane of rock debris mixed with mud in T1yn1-2 layer. With the increase of dissolution rate of T1yn1-1 limestone, the leakage of dam foundation increases significantly, and there is a positive correlation between them. The regression equation is y2=120.3x+224.8. The concentrated belt of solution crack is the main leakage area of dam foundation. When the concentrated belt develops horizontally and distributes intensively along the contact belt between T1yn1-1 limestone and T1yn1-2 marl, the anti-sliding stability coefficient of dam foundation will significantly reduce and the leakage of dam foundation will significantly increase.Therefore, the concentrated belt should be treated as a key area. When the concentrated belt is vertically developed, dispersed or mainly distributed in the area behind the dam, it has little impact on the anti-sliding stability and leakage of the dam foundation, and can be used as a secondary treatment area. Karst hydrogeological analysis and numerical simulation show that in complex karst areas, after groundwater is discharged from the surface in the potential catchment area, it mainly influxes in the form of runoff to the downstream river channel, and the vertical infiltration of water flow is relatively weak, so it is difficult to form karst pipelines bypassing the anti-seepage curtain and connecting the upstream and downstream in the deep part of the riverbed. The leakage form of dam foundation is mainly solution crack leakage, and its engineering impact is limited. Therefore, the catchment area of karst groundwater is suitable for dam and reservoir construction. In addition, according to the spatial distribution characteristics of the concentrated belt of solution crack, the targeted treatment of zoning grouting can improve the treatment efficiency and save investment.
Study on "three-section and four-layer" reinforcement technology of tunnel vault subsidence in the karst basement
QIN Jianwen, LI Haibo, PAN Guangming, TIAN Yueming
2022, 41(5): 746-758. doi: 10.11932/karst20220508
Abstract:
In order to effectively control the geological disaster of tunnel vault subsidence in the karst basement, the "three-section and four-layer" control technology is proposed, based on the full analysis of engineering geology and hydrogeological characteristics of the subsidence stratum. In this technology, the treatment area is divided into the key reinforcement section, the sub-key reinforcement section of the affected area and the pre-excavated reinforcement section. According to the treatment depth and order, the key reinforcement section of the subsidence area is divided into the top slurry barrier layer, the vault cushion layer of slurry-stopping, the vault reinforcement layer and the intermediate reinforcement layer. In terms of grouting mechanism, filling grouting and splitting-compacting grouting are adopted in both the top slurry barrier layer and the intermediate reinforcement layer. The alternate mode with filling-splitting-compacting grouting is implemented in the intermediate reinforcement layer, and the filling grouting mode is adopted in the vault cushion layer of slurry-stopping. Splitting-compacting grouting is conducted in both the pre-excavated reinforcement section and the sub-key reinforcement section in the affected area. In order to reduce the leakage of grout in the top slurry barrier layer and vault cushion layer of slurry-stopping, C-GT quick setting mixed grout is selected, and the initial setting time of the grout is controlled at 30-50 seconds, while the final grouting pressure is no more than 0.5 MPa and 0.8 MPa. In order to realize the effective diffusion of slurry in the vault reinforcement layer and intermediate reinforcement layer, the initial setting time of C-GT mixed slurry is extended to 50-90 seconds, even up to 150 seconds, and the final grouting pressure is up to 0.8-1.2 MPa. To improve the slurry diffusion and overall reinforcement strength of the pre-excavated reinforcement layer and the sub-key reinforcement layer in the affected area, the alternate and continuous grouting mode of single and double cement slurry is implemented based on the high consolidation strength of single cement slurry and the quick setting characteristics of C-GT mixed slurry. The above-mentioned technology is the key to the success of grouting reinforcement. Besides, the combined grout-stopping method is used to effectively deal with the segmented grout-stopping at the place where the hole wall collapses and shrinks, and the diameter of borehole is greatly oversized. Hence, the isolated grouting at different depths of those four layers is realized. The directional inclined hole is used to effectively avoid the hidden underground pipeline, and the grouting pipe is used as the pipe roof to improve the shear resistance of stratum. Combined with grouting, this technology can effectively improve the overall strength and stability of the reinforcement stratum by tunnel vault grouting. Research results indicate that "the three-section and four-layer" treatment scheme is the basis to effectively control the subsidence of soil overlying the tunnel. The combined grout-stopping technology in the hole is the guarantee to meet the requirement of segmented grouting at different depths. The combination of vertical hole and directional inclined hole and of filling grouting and splitting-compacting grouting, taking quick setting slurry as the primary part supplemented with single cement slurry, and the strict control of safe grouting parameters are the technical keys of safe and effective grouting. The application of the above technology has achieved good grouting reinforcement effect and hence should be popularized.
Karst development characteristics and distribution law in Huaihua basin, Hunan Province
WU Yuanbin, LIU Zhikui, YIN Renchao, YANG Jianxing, LUO Weiquan, LEI Mingtang, DAI Jianling, PAN Zongyuan
2022, 41(5): 759-772, 807. doi: 10.11932/karst2022y22
Abstract:
Huaihua basin is located at the center of the marginal areas of Hunan, Hubei, Chongqing, Guizhou and Guangxi. Shanghai-Kunming high-speed railway, Zhangjiajie-Jishou-Huaihua high-speed railway, Huaihua-Guilin high-speed railway (under planning), Chongqing-Huaihua high-speed railway and Huaihua-Shaoyang-Hengyang high-speed railway converge in Huaihua basin in the shape of Chinese character "mi", which makes Huaihua City a new national comprehensive transportation hub. Located in the zone of temperate subtropical monsoon climate, the study area is humid with abundant rainfall. The average rainfall over the year is 1,326.3 mm, with the maximum of 1,663 mm. Rainfall is mainly concentrated in April to August, accounting for more than 70% of the average annual precipitation. The landforms of the study area are mainly middle and low mountains and karst valleys. In the west is Zhongpo mountain-Shiziyan mountain range, and in the east is Liangshan mountain range. The overall topography is high in the northwest and southeast, and low in the middle, showing a basin form spreading in the NNE to NNE. The pure carbonate strata in the study area from old to new sequence are middle Dapu formation of Carboniferous (C2d), upper Maping formation (CPm), middle Qixia formation of Permian (P2q), Maokou formation (P2m), upper Wujiaping formation (P3w), second section of lower Daye formation of Triassic (T1d2), and lower Shimen formation of Cretaceous (K1s). The regional structure in the area belongs to the Neocathaysian tectonic system, which is located in the contact zone of thrust-fold belt between the Yuanma basin and the Xuefengshan. The main structure in the area is the north (north) east compressive structure, and the others are associated structures. The geological structure of Huaihua basin is very complex and karst is strongly developed. However, the previous studies mainly focus on the basic introduction of karst development characteristics in a small area. In addition, the 1∶50,000 regional geological-and-hydrogeological survey has not yet covered the study area. Therefore, the geological research in this area is generally not profound, and the existing geological data is difficult to meet the needs of urban planning and construction. Taking carbonate rocks in the area as the research object, a large number of previous geological data has been collected and analyzed. New geological investigation and drilling have been used to explore the space characteristics of strata and lithology combination. The characteristics of karst development in Huaihua basin have been systematically summarized from six aspects, namely, karst morphology, strata and lithology, tectonic conditions, base level of erosion, landform and hydrodynamic conditions, the contact zones of soluble and non-soluble rocks. Furthermore, Huaihua basin has been zoned according to the development degree of karst. The research findings show as follows, (1) The spatial distribution of soluble rocks in the study area is not uniform, mainly in covering karst areas and partially in buried-karst areas. Solution grooves, stone teeth and sinkholes are developed in the whole area, while the underground river is mainly developed in the tectonically dissolved peak ridge depression in the Baisha-Caimenxi area in the eastern part of the study area. (2) Karst development in the area is mainly controlled by the stratigraphic group, and the next controlling factor is lithology. The karst development intensity of each layer in the valley is ranked as thick-layer limestone (CPm, P2m, P2q, P3w)>dolomite (C2d)>red-bed gray conglomerate (K1s )>thin-layer limestone (T1d 2). However, influenced by geomorphology and hydrodynamic conditions, the karst development of dolomite (C2d ) in tectonically dissolved peak ridge depression is more intense than that in valley. (3) From high to low, there are three most developed segments of dissolution in the study area, with elevations of 500-600 m, 210-270 m and 190-210 m respectively. These three segments are respectively controlled by base level of historical erosion, base level of partial erosion in Luyang-Huaihua-Zhongfang valley and base level of partial erosion of Wushui river. Among them, the main karst forms in the segments at the depth of 210-270 m and of 190-210 m are karst springs, underground rivers, dissolved pores, solution cracks and caves under the overburden of valley. (4) The karst development is more intense in areas with strong alternating groundwater circulation, such as the banks of Wushui river, the fault zone, the contact zone between soluble and non-soluble rocks, and the piedmont alluvial area. The karst area is mainly distributed in Luyang-Chengdong new area and Baisha-Caimenxi area in the southeast of Luyang valley. If engineering activities such as groundwater exploitation, drainage of foundation pit, pile foundation construction, tunnel construction, etc. are planned to carry out in these areas, special risk assessment of geological disasters caused by karst collapse in construction site is needed.
Impact of controlling karst rocky desertification on soil particulate organic carbon and aggregate-associated organic carbon
LAN Jiacheng, WANG Junxian, WANG Shasha, QI Xue, LONG Qixia
2022, 41(5): 773-783. doi: 10.11932/karst20220509
Abstract:
The control measures of karst rocky desertification exert important influence on soil organic carbon (SOC) composition, and then affect the accumulation and stability of organic carbon. However, the effects of controlling karst rocky desertification on soil particulate organic carbon (POC), mineral-associated organic carbon (MOC), and their relationship between SOC and aggregate-associated organic carbon are still unclear. Huajiang karst gorge area is one of the most typic demonstration areas of controlling karst rocky desertification in Guizhou Province, Southwest China. Before the 1990s, this area underwent extensive land degradation, which led to the acceleration of SOC emissions. The local people have developed several well-known control measures of rocky desertification, among which we selected two-conversion of cropland to secondary forest and to Zanthoxylum plantation-as study objects. Given cropland as reference, soil was collected in the layers at the depth of 0-20 cm to analyze the impact of the two selected measures on SOC, POC, MOC and aggregate-associated organic carbon as well as their relationship. The results show that compared with cropland, the concentrations of SOC, POC, MOC and aggregate-associated organic carbon at the depth of 0-20 cm increase both in Zanthoxylum plantation and secondary forest. The concentrations of SOC and aggregate-associated organic carbon in secondary forest are significantly higher than those in Zanthoxylumplantation and cropland in the layers at the depth of 0-10 cm (P<0.05), but no significant difference is shown in the layers at the depth of 10-20 cm (P>0.05). The POC concentrations in both segments (0-10 cm and 10-20 cm) significantly increase inZanthoxylum plantation and secondary forest, but the MOC concentrations show no significant changes. POC/SOC ranging from 20.38% to 45.27% is significantly higher in Zanthoxylum plantation and secondary forest than that in cropland ( P<0.05). On the contrary, MOC/SOC in cropland is significantly higher than that inZanthoxylumplantation and secondary forest ( P<0.05). After the conversion of cropland toZanthoxylum plantation and secondary forest, the SOC concentrations have increased mainly due to the increase of POC concentrations. The contribution rate of the particle size bigger than 2 mm to SOC in Zanthoxylum plantation and secondary forest is significantly higher than that in cropland. However, the contribution rate of the particle size between 0.25-2 mm, between 0.053-0.25 mm and smaller than 0.053 mm respectively to SOC is significantly lower than that in cropland. The correlation analysis shows that POC and MOC are positively correlated with SOC and aggregate-associated organic carbon. Their correlations are listed as follows,secondary forest>Zanthoxylum plantation>cropland. After the conversion, the increase of SOC, POC and MOC is significantly correlated with the increase of aggregate-associated organic carbon (P<0.05), with higher correlation in the secondary forest. The control measures of rocky desertification have changed SOC and its physical composition and their relationship, thus promoting the accumulation of organic carbon.
Spatial variation of soil acidification in the karst area of northern Guangdong: A case in peak cluster depression and karst trough valley landforms
WU Yuheng, ZENG Faming, WEI Xinghu, LIU Shujuan, HUANG Jinguo
2022, 41(5): 784-795. doi: 10.11932/karst20220510
Abstract:
The process of soil acidification is very slow under the natural conditions, which generally takes decades or even hundreds of years. However, since the Industrial Revolution, human activities have greatly accelerated the process of soil acidification. The effect of acid rain is one of the main reasons for the process. The problems caused by soil acidification, such as barren soil, soil quality decline, enhanced heavy metal activities and others will affect soil microbial activities and the growth of plants, and then affect the quality of ecological environment in the whole region, hence threating the food safety and human health. Soil acidification has become an important ecological and environmental problem that has to be faced in the agriculture development of south China. Therefore, distribution of soil pH is the key to control measures for soil acidification in karst areas. In order to study the spatial variation of soil acidification in karst areas, an investigation was conducted to assess the status of soil pH in the karst trough landform and cluster-peak depression landform in northern Guangdong located in the acid rain area of south China. Methods such as normalized vegetation index (NDVI), normalized moisture index (NDWI), normalized building index (NDBI), slope aspect and gradient were used to analyze the driving mechanism of soil pH heterogeneity in cluster-peak depression and karst trough valley. The research questions are listed as follows, (1) What are the characteristics of spatial variation of soil pH in different karst landforms and what are the main influencing factors of the variation? (2) What are the driving factors and interactive effects of soil acidification in karst areas? The results show that pH values of the total soil samples range from 4.62 to 8.19 with the mean of 6.29±0.92 and the variance of 0.85; the soil pH values in karst trough valley range from 4.62 to 8.19 with the mean of 6.38±0.95 and the variance of 0.91; the pH values in cluster-peak depression range from 5.29 to 8.03 with the mean of 6.31±0.90 and the variance of 0.81; the pH values in non-karst area range from 4.71 to 7.55 with the mean of 5.99±0.92 and the variance of 0.73. Generally speaking, the average soil pH of total samples shows weak acidity, and no soil samples with strong acidity (pH<4.5) are detected, but alkaline soil samples (pH≥7.5) in some landforms are detected. The results indicate a substantial impact of acid rain on the soil pH of karst areas and the risk of soil acidification. According to the distribution frequency of soil pH values, the soil pH in the whole study area is mainly weakly acid and acid. The pH values (<6.5) of soil samples in different landforms are ranked as follows, karst trough valley (61.0%)<cluster-peak depression (69.2%)<non-karst area (75.0%), indicating that soil acidification in the study area has been widespread and its control is urgent. Through the analysis of the spatial characteristics of soil pH, the variation of the mean value of soil pH presents the following rules, (1) Distribution of the soil pH shows that soil parent materials directly affect the process of soil acidification. Under the conditions of different soil parent materials, the mean values of soil pH are ranked as, non-carbonate<dolomitic limestone<dolomite<limestone, indicating that the buffer capacity of carbonate rock is stronger than that of non-carbonate rock. (2) Soil pH values in the study area increase with the rise of gradient, but the factor of slope aspect does not show an obvious influence on distribution of soil pH. These results indicate that the influence of the gradient factor in the study area is stronger than that of the slope aspect. (3) In the karst trough valley, the mean values of soil pH in different land use types are ranked as, woodland>grassland>cultivated land, indicating that the land use type, in some degree, is influential in the spatial variation of soil pH. The karst environment is very complex. With the effect of acid rain, the soil parent materias, which may be the key factor to control the soil pH in the study area, directly affect the process of soil acidification. But distribution of soil pH spatial variation may be caused by multiple factors. In general, soil pH values in the study area show an increase with the rise of gradient. The influence of gradient on soil pH spatial variation is followed by that of soil parent materials. The interactive effect of gradient and soil parent materials shows great explanatory power for the spatial variation of soil pH in the study area. Therefore, soil acidification control is advisably to be carried out in northern Guangdong, and the gradient factor should be taken as one of the indicators. Besides the influence of acid rain, other exogenous factors such as nitrogen fertilizer, and industrial and mining activities, should also be paid attention to. With the integration of geostatistics, the geodetector is used to further analyze the driving factors of soil pH spatial variation and their interactive effects in the karst area in north Guangdong, based on the forecast of soil pH kriging spatial interpolation. However, the effect of rocky desertification and the scaling effects of different sizes of spatial statistical units on the geodetector model have not been fully explored. These limitations will be overcome in the future research.
Simulation of buffering process and carbon sink effect of lime soil on sulfuric acid rain
ZHAO Guangshuai, HUANG Qibo, ZHU Yinian, LI Tengfang, PU Zhenggong
2022, 41(5): 796-807. doi: 10.11932/karst20220504
Abstract:
After sulfuric acid rain settles to the surface and is buffered by lime soil, its participation in carbonate rock erosion and its impact on karst carbon sink are still unclear, which seriously restricts the accurate assessment of karst carbon sink effect in China. In order to clarify the buffering process and key control factors of lime soil to sulfuric acid rain, we conducted leaching experiments of sulfuric acid rain with pH=4.5 under different soil thicknesses. Results show that contents of Ca2+, Mg2+ and ${\rm{HCO}}_3^{-}$ in leaching solution decrease rapidly at the initial stage of leaching. When the leaching amount (after soil water reaches saturation) reaches 1,020 mL, the leaching loss of Ca2+, Mg2+ and ${\rm{HCO}}_3^{-}$ tends to be stable, and stable leaching amounts are 20 mg·L−1, 6 mg·L−1 and 40 mg·L−1 respectively. The same ion in leaching solution of different thicknesses of soil columns shows the same trend, which indicates that the lime soil buffering of acid rain may mainly occur in the 10 cm (surface) soil layer. Ca2+ and Mg2+ exchanged by acid rain H+ in leaching solution only accounts for a small part, and the soil water-soluble Ca2+ and Mg2+ is the main part of Ca2+ and Mg2+ in leaching solution. In an open system, atmospheric and soil CO2 dissolves in rainfall to form H2CO3, which not only increases carbon sink, but also exchanges Ca2+ and Mg2+ in soil with H+ generated by the dissociation of H2CO3. The leaching amount of Ca2+ and Mg2+as a result cannot be ignored. The special physical structure of topsoil may produce preferential flow, which will significantly reduce the exchangeable point of H+ in leachate, causing dissolution of carbonate minerals to assist the buffering of H+ in leachate. The karst area is the main place for carbon loss in terrestrial ecosystems, and the most intense carbon leaching occurs in surface soil. Exchangeable calcium, exchangeable magnesium and carbonate minerals in lime soil are principal reactants for buffering acid rain. The buffering capacity of exchangeable calcium and magnesium in lime soil with different thicknesses is greater than that of soil calcium carbonate minerals, and the former is 1.17-1.59 times as large as that of the latter. Under the same acidity and the same rainfall (after the soil water reaches saturation), the carbon sink generated by soil base ions participating in acid rain buffering is about 2.1 times as large as that generated by weathering buffering of carbonate minerals. Under the same rainfall conditions, the existence of lime soil can significantly increase carbon sink, and the carbon sink generated by lime soil with different thicknesses (≥10 cm) is roughly equal. According to this experiment and the rainfall data of Guilin, in the weathered residual soil area of pure limestone in Guilin (soil thickness greater than 10 cm), the participation of soil base ions in acid rain buffering can produce 0.59-0.93 mol·m−2 carbon sink flux per year.
Research progress of paleokarst oil and gas reservoirs
ZHANG Yu, ZHAO Lun, LI Changhai, ZHANG Xiangzhong
2022, 41(5): 808-824. doi: 10.11932/karst202204y10
Abstract:
The transformation of soluble rocks by paleokarst can form dissolution pores and other reservoir spaces. Especially for carbonate rocks, the development of paleokarst can greatly improve the reservoir performance of the formation and form karst reservoirs with good porosity and permeability. Paleokarst related oil and gas reservoirs are developed in many large and super-large carbonate oil and gas fields at home and abroad, such as Yurubcheno Oilfield in Russia, Yates Oilfield in the United States, Halfaya Oilfield in Iraq, Renqiu Oilfield, Tahe Oilfield, Maokou formation and Dengying formation in Sichuan Basin and Ordovician formation in Ordos Basin in China. As one of the important reservoir types in carbonate formation, the development and distribution of karst reservoirs are of great significance in the study on carbonate reservoirs. This paper systematically summarizes the previous studies on different types of karst reservoirs and the latest research progress of paleokarst oil and gas reservoirs, which are of important reference value for understanding the development and distribution of karst reservoirs in oil-gas basins in China. The following problems and results can be concluded in the study of paleokarst oil and gas reservoirs. There are few studies on early diagenetic karst reservoirs, and most of them focus on the formation mechanism and development characteristics of eogenetic karst reservoirs. In addition, there is no corresponding practice on how to use logging and seismic data to predict the distribution of this type of reservoir. There are many studies on supergene karst reservoirs, mainly focusing on paleogeomorphology restoration and vertical zonation. Meanwhile, many palaeogeomorphology restoration methods have been proposed, such as impression method, residual thickness method, stripping and filling method, sedimentological analysis method, sequence stratigraphy method, layer leveling method, double-interface method, well-seismic combination method and the new seismic geomorphology method. However, all of the present restoration methods of paleokarst geomorphology have their advantages and disadvantages in certain applicable scope, and can only be used to restore paleokarst geomorphology qualitatively and semi-quantitatively. How to precisely quantify the restoration of paleokarst geomorphology remains further studied. In terms of vertical zonation, many schemes have been proposed from the perspectives of runoff strength, karst space form, karst development intensity, and the characteristics of hydrodynamic action. In particular, a great deal of research and practice has been carried out on the characteristics of hydrodynamic action, and a variety of vertical zonation schemes have been put forward to divide karst vertically into 2-6 zones. But in general, although there are many schemes to divide the vertical zonation of paleokarst, a unified classification has not yet been formed. Besides, the boundary between the vertical zonation and the zonation is not clear, so it is unlikely to distinguish them quantitatively and accurately. From the case study of Sichuan Basin and Tarim Basin in China, buried karst reservoirs have significant identification marks in mineralogy, petrology and geochemistry. In the prediction of buried karst reservoir, the current research is mainly focused on migration channel and fluid distribution. From the perspective of reservoir formation, the research on buried karst reservoirs at present is mainly focused on sulfate thermal reduction and hydrothermal dissolution, but there is a lack of relevant research on buried karst reservoirs affected by organic acids. As a hotspot of karst reservoir research in recent years, the studies on fault-karst reservoirs are mainly focused on the development mechanism and characterization of fault solutions, etc. The recognition, characterization and prediction of fault-karst reservoirs mainly rely on seismic data, and are significantly affected by the resolution of seismic data. There are many internal and external factors affecting the development of paleokarst oil and gas reservoirs. The internal factors include lithology, development degree of primary pores, rock structure, mud and acid insoluble content, fractures, etc. The external factors mainly include sedimentary facies, rise and fall of sea levels, palaeogeomorphology, tectonic movement and paleoclimate. In terms of the prediction of paleokarst oil and gas reservoirs, the accurate prediction of karst reservoirs cannot be achieved with single geological or seismic data, because there are many factors affecting the development of karst reservoirs and many differences in the formation mechanism and development characteristics of different types of karst reservoirs. The comprehensive prediction method based on geology-log-seismic can be adopted as the mainstream of karst reservoir prediction in the future. With the continuous development of geophysics and artificial intelligence of oil and gas, the increasing new methods and technologies will be applied to karst reservoir prediction.
Study on the delimitation of affected zone of geological environment for karst underground engineering:taking Longgang district, Shenzhen City as an example
WU Shengtang, JIANG Xiaozhen, MA Xiao, TANG Zhen
2022, 41(5): 825-837. doi: 10.11932/karst20220511
Abstract:
The construction of a large number of urban underground rail transit and urban renewal projects has intensified the urban karst geological disasters in our country. As one of the first demonstration areas of national urban construction, Shenzhen has developed rapidly in terms of underground rail transit and urban renewal projects over the years, hence leading to frequent karst collapse disasters due to its location in the karst area. Therefore, the summary of the experience and lessons from the development of karst underground space is of great significance for the engineering construction in karst areas. In this thesis, a preliminary study on the delimitation of the affected zone of geological environment for karst underground engineering is conducted based on three events of karst ground collapse in Shenzhen. The affected zone of geological environment for karst underground engineering refers to the area where disasters are likely to happen due to the disturbance and damage of the rock and soil around the construction site during the construction process. The delimitation of the affected zone is not only conducive to the safe and smooth engineering construction, but also to the clear division of responsibility. The monitoring range of underground engineering construction in karst area usually reaches only tens of meters at the current stage. But when a disaster happens, the actual range influenced by engineering will often exceed hundreds of meters. Therefore, the current construction specification about the affected zone of geological environment for karst underground engineering is unreasonable and uncertain in some degree, and the relevant provisions are greatly challenged. For instance, it is stipulated that the engineering monitoring should be conducted within the plane range that is only 3 times as deep as the foundation pit during construction , and the description of the expansion of the monitoring range in the karst development area is not detailed. Therefore, the further research on the affected zone of geological environment for karst underground engineering is very necessary, so the actual affected range of the project can be effectively judged, and then the corresponding prevention and control measures can be taken. On the basis of fully mastering the regional geological background and the geological conditions of site engineering, a preliminary study is conducted in this thesis. Firstly, the actual affected range of the project is monitored and analyzed by the high-frequency monitoring technology of karst groundwater pressure. The monitoring scheme should be formulated according to local conditions. The monitoring frequency should capture the disturbance changes of regional karst groundwater with more than 3-month monitoring cycle. The results of monitoring and data measurement of water levels indicate that the affected range of the project is closely connected with the karst groundwater drawdown funnel. The obvious anisotropy of karst aquifer medium at each site is indicated in the groundwater flow field, which is mainly controlled by karst development and structure. The maximum affected ranges of the three projects are 560 m, 820 m and 850 m respectively. Then, the analysis on the formation mechanism of karst collapse is conducted. Results indicate the collapse mechanism. The excavation and precipitation of foundation pit leads to the change of groundwater hydrodynamic internally caused of strong karst development, the change of hydrodynamic conditions disturbs karst groundwater or gas, and the force generated by karst water disturbance acts on the overburden floor through karst pipeline or crack. As a result, the overburden soil mass collapses and loses gradually, until the roof becomes unstable and damaged for the insufficient collapse resistance. Finally, the Gehart’s empirical formula of influence radius of confined water pumping hole is used for reference, combined with engineering construction parameters, main factors of karst collapse and hydrogeological test parameters. The mixed research method of qualitative analysis and quantitative calculation is adopted to theoretically deduce the empirical calculation formula of the affected zone of geological environment for karst underground engineering. Research results indicate that the affected range of karst underground engineering is mainly in direct proportion to the permeability coefficient and engineering depth. However, it is inversely proportional to the thickness of soil layer. The deduced semi-quantitative theoretical formula is suitable for the calculation of a relatively thick aquifer between the Quaternary and the karst aquifer. If the depth of foundation pit is greater than that of rock surface and of confined karst water in underground engineering, this fomular can be used to quickly provide the basis for safety construction of underground engineering as well as urban disaster prevention and reduction in the karst area lacking groundwater monitoring data.
Carbon and oxygen isotope characteristics and paleoenvironmental significance of deep karst fracture-cave fillings in Huanjiang sag, Guangxi
DONG Hongqi, ZHANG Qingyu, LIANG Jiapeng, DAN Yong, LIANG Bin, NIE Guoquan, JI Shaocong
2022, 41(5): 838-846, 837. doi: 10.11932/karst20220512
Abstract:
In recent years, deep karst caves are commonly met in deep resource exploration and engineering construction. However, the deep-buried caves lead to the insufficient research on the genesis and development mechanism of deep karst caves. As good research samples, a large amount of complete drilling cores in deep karst caves are founded in Huanjiang sag in Guangxi. The karst morphology analysis of the drilling core in Well HD1-4 reveals that the deep karst in Huanjiang area is mainly composed of net cracks and holes expanding along cracks and dolomite honeycomb pores, and large karst caves are also developed, with the maximum height of 20 m. The types of fillings in deep karst caves are diverse. Through the observation of the whole well core, it is found that the deep karst cave fillings are of four characteristics. (1) Calcification growth,the mixed growth pattern of chemical and argillaceous substances reflects periodic growth in the cave, which may be associated with surface hydrological systems. (2) Flower-like growth pattern,white calcite with multi-stage growth can be identified in pores, but the growth direction and growth environment of each stage are different, hence forming flower-like growth pattern. (3) Primary chemical fillings of the hole,the hole is filled with calcite or dolomite, and the crystal form of calcite in part of the hole is good. (4) Argillaceous fillings in the cave, argillaceous fillings can be seen at 230 m, 432 m, 880 m and even 1,132 m in Well Huandi 1-4, which are far lower than the local discharge datum or the sea level. These argillaceous fillings may come from the surface mud seepage along the fracture or the mud beneath the ancient exposed surface. In this paper, the carbon and oxygen isotope analyses of deep karst cave fillings and surface karst cave fillings in Well HD1-4 drilling show the wide distribution of carbon and oxygen isotope of deep karst fracture. δ13C values are between −5.2‰ and −2‰ with the average value of −0.33‰. δ18O values are between −16.78‰ and −5.3‰ wtih the average value of −11.45‰. The values show the general negative skewness. The negative skewness of carbon and oxygen isotope of calcite fillings in the pores is the largest, and that of dolomite is the smallest. The carbon and oxygen isotope values of calcareous mudstone are the closest to those of modern atmospheric freshwater. Based on the analysis of geological conditions, four large-scale paleokarst processes and filling periods are founded in the formation and filling stages of deep karst in the Huanjiang area. (1) In the karst environment in contemporaneous period and penecontemporaneous peiord, the distribution range of δ18O is the same as the background value of the bedrock of carbonate rock, but δ13C changes greatly. The results show that the karst environment is similar to the sedimentary environment of carbonate rocks, which reflects the short-term exposure of karst after the deposition of carbonate rocks in the contemporaneous period, and the filling shows the characteristics of early precipitation. Karst space is mainly characterized by dissolution pores, which provides the basis for further karst development. (2) The karst environment of atmospheric fresh water in the hypergene period is affected by atmospheric fresh water δ13C, and the δ18O value shows a significant negative skewness. The fact that δ13C is less than −2.5‰ and δ18O is between −13 and −8.2‰ indicates the karst environment is an open system, and a large number of argillaceous substances infiltrate with surface water and fill karst caves, making a certain impact on the preservation of deep resources. (3) In the shallow-buried karst environment, the δ18O value of the filling is negative, which is more negatively skewed than that of Type I bedrock. The δ13C value is basically consistent with the bedrock value, and there is no obvious negative skewness. The δ13C value is between −2.0‰ and 2.5‰, and the δ18O is between −13.0‰ and −9.0‰. Due to the increase of temperature and pressure in the closed system during the burial peirod, there gradually precipitates and forms dolomite. (4) In the deep-buried or hydrothermal karst environment, the δ18O value of fillings is obviously negative, and high-temperature fluid flows up along the fault to form the karst dissolution space which becomes the reservoir place of various deposits. The hydrothermal karst with high temperature and the atmospheric freshwater karst are the main periods of the formation of deep karst fracture-cavity filling in Huanjiang area. The research results are of great significance for deep karst reservoir prediction and deep resource exploration in the later stage.