裘布依地下水井流计算公式有关问题讨论

在长期的实际工作中发现,抽水井口径对抽水试验成果存在较大影响、传统的裘布依井流公式对井涌水量和水文地质参数计算结果与实际情况偏差较大。本文以不同水文地质条件下岩溶含水层中井抽水试验成果资料为依据,通过统计分析,揭示了井径与井涌水量以及含水层水文参数计算成果之间的相关关系,探讨了裘布依井流计算公式在理论和实践上存在的问题...     

岩盐井水泥固井工艺

岩盐井水泥固井工艺甘肃省地矿局第一水文地质工程地质队杨志刚我队在甘肃省漳县盐井矿区的施工中，在无专门水泥固井设备的情况下，利用自己设计加工的一套固井水泥搅拌及灌注设施共施工岩盐井５个，获得了较好的效果。１岩盐井设计要求岩盐井设计最深为４００ｍ，最浅为...     

地下水源热泵系统中回灌能力分析

近年,地下水源热泵技术在我国被广泛应用,并在节能、环保等方面取得了一定效益。但是,回灌问题仍是困扰我国地下水源热泵发展的瓶颈。以唐山市丰润区乡居假日住宅区A4区地下水源热泵系统的应用为例,从区域水文地质条件方面,对水源热泵系统中地下水回灌能力进行了分析,指出开展地下水源热泵项目时,掌握热源井所在区域水文地质条件的重要性。探讨了影响地下水回灌能力的关键因素,其中包括区域水文地质条件、热源井成井工艺、回灌井阻塞以及地下水回灌方式。     

老矿山水文地质条件及开采后变化特征——以水城观音山铁矿区为例

以水城观音山铁矿区为例,在收集资料和调查的基础上,对矿区水文地质概况;开采前后的充水因素、补径排条件、各矿体静止水位标高;开采后的水文地质变化特征等进行了分析、总结;重新划定了水文地质勘查类型。其结果为:矿山经过多年的开采,水文地质条件发生了极大变化,表现为:地下水位、排水量大幅度下降;井、泉大量干涸;河流、沟谷流量明显减少;但采空区积水、入渗补给条件、充水因素等发生变化。对老矿山开采后的水文地质特征适时动态认识与评价,为正确指导矿山安全生产提供有价值的水文地质资料。     

降水效果的数值分析和设计优化研究

降水在各类工程中使用的越来越多,已形成了一套设计和施工方法,但其降水效果分析却少有人研究。针对某露天煤矿开挖中的降水疏干工程,开展了大规模的数值模拟研究。在掌握水文地质条件的基础上,通过Modflow建立了较合理的数值模型,对不同工况的降水效果进行了详细的分析,实际计算表明,加井时间、调整单井抽水量、井的数目、井的排列以及它们的相互作用对降水效果有很大影响。分批次、分时段、优化井排列和调整单井抽水量是提高降水效果的设计措施和优化设计的重要影响因素,这些因素的作用规律将为降水设计和施工以及实际工程效果提供重要理论依据。     

地下水源热泵抽灌井优化布置及参数灵敏度

为了确定地下水源热泵的最佳布井方式,探明含水层温度场变化特征以及水文地质参数灵敏性,以石家庄地下水源热泵适宜区肖家营—东兆通为研究区,基于热储层的地热地质条件和热物性参数,利用FEFLOW软件构建了三维水-热耦合数值模拟模型。在2 000 m³/d的抽灌量下对不同井间距、不同布井模式-含水层温度场变化进行模拟,探讨发生热突破的可能性,并探讨温度场关于热导率、渗透系数、孔隙率等参数的灵敏度。研究结果表明：1)在单抽双灌模式下,方案C(两眼回灌井连线与地下水流向垂直,一眼回灌井位于抽水井正下游,另一眼回灌井位于沿地下水流向45°方向,即两回灌井与抽水井组成等腰直角三角形,抽水井正下游的回灌井为直角顶点)的布井方案对含水层温度场影响最小,为最优方案;2)单抽双灌模式下,抽灌量为2 000 m³/d时,方案C抽灌井间距设置在40～50 m较为合理;3)抽灌过程中,温度场对于渗透系数的改变灵敏度较高,而对孔隙率以及热导率的改变灵敏度较低。     

深厚强透水含水层超大基坑降水群井效应研究

基坑减压降水的幅度与群井效应密切相关。某超大面积基坑,含水层组厚40m,中下部透水性强,采用非完整井降水。对水文地质条件概化,建立了地下水三维非稳定流数值模型,对均匀布置群井、不均匀布置群井、分块开挖降水、不同的井结构、布设回灌井等工况,进行了渗流场模拟。研究表明,强透水性含水层超大面积基坑降水的群井效应极为明显;井位角密中疏布置,可实现降深调平,避免降深不足和超降;基坑分块降水,可减少坑外降深;短滤管井结构可减少基坑总涌水量和坑外降深;强透水性含水层可灌性强,回灌对减少坑外水位下降有较明显的效果。模拟结果与现场监测较为一致。研究成果可为类似工程地下水控制设计提供参考。     

Visual Basic在求承压含水层水文地质参数中的运用

在广泛了解求水文地质参数方法研究现状的基础上,结合VisualBasic程序设计理论,用解析法、标准曲线比拟法和直线图解法求承压含水层水文地质参数的原理和程序设计步骤。以鄂尔多斯盆地某地浸砂岩型铀矿床勘查的完整井和非完整井以及稳定流和非稳定流抽水试验资料为依据,分别计算了该矿床含矿含水层的水文地质参数,并对几种方法的计算原理、过程和结果进行了对比分析。     

抽水试验中不同位置自动水位计响应数据应用分析

抽水试验中,动水位数据采集记录及处理分析对水文地质参数计算具有重要意义。近年来自动水位计被广泛用于抽水试验,通过传感器压强水头变化值获取水位降深。因井管内水流动会产生水头损失,自动水位计安放位置不同会导致获取的井水位降深不同,不同于传统方法测得的井水面降深,对水文地质参数计算将产生一定影响,因此如何合理放置自动水位计以及在参数计算中如何应用其获取的水位降深都亟待开展试验研究。在黑河流域第四系大厚度含水层地区,选择典型单层试验孔和利用分层封隔技术实现的一孔同径多层抽水孔开展试验研究,在动水位以下抽水试验层段上部、中部、下部以及潜水泵上部和下部分别放置自动水位计进行了系统的数据采集分析。结果表明:抽水试验中因井管内水流沿程水头损失及速度水头差异导致不同位置自动水位计获得压强水头变化值不同,本次试验实测到井筒内不同部位井损值;井损值在潜水泵进水口处最大,随距潜水泵距离的增大而减小,为避开井筒内较大水头损失对参数计算的影响,自动水位计宜优选安放在潜水泵上部接近动水位位置;在单孔抽水试验中利用稳定流公式计算水文地质参数时,自动水位计获取水位降深含井损不可忽略,需通过多落程抽水试验数据分析扣除后使用。同时,抽水试验中自动水位计不同位置获取数据的处理分析方法为更好地理解井中水头损失提供了依据。     

空气潜孔锤钻进技术在乌蒙山地下水探采结合工程中的应用

乌蒙山地下水探采结合工程是将水文地质调查中的水文地质勘探与供水井施工相结合的一项国家重点扶贫项目,该项目50余口探采结合井的施工解决了超过5万人的生活饮用水问题,取得了良好的社会经济效益。空气潜孔锤钻进技术在该项目中发挥了重要作用。文章介绍了空气潜孔锤在水文水井钻探中的技术优势,重点阐述了该项技术的施工工艺、设备选型以及典型复杂情况的处理,并提出了该项技术在水文水井施工方面的下一步研究方向。     

露天煤矿疏干–回灌水帷幕保水采煤模式探讨

为抑制疏干排水造成的地下水水资源量衰减,并控制露天煤矿煤炭开采对地下水资源的影响范围,采用了水帷幕保水采煤技术。以某露天煤矿为背景,建立沙槽回灌平台进行物理实验,再通过数值方法对回灌区进行模拟计算,分析疏干–回灌水帷幕的形成和变化规律。物理实验方法观察到回灌水帷幕形成定水头边界,阻止漏斗扩展的过程,获得不同回灌水帷幕水头高度对注水量的影响和注水量与排泄量平衡的最佳水头。数值模拟获得物理模拟类似的结果,回灌水帷幕位置可以有效限定疏干降落漏斗范围,回灌水帷幕定水头高度与疏干区水位和疏干时间正相关;最大回灌量受帷幕定水头高度影响显著;回灌水帷幕距离采坑越远对阻隔帷幕外水源补给的作用越大。物理模拟和数值模拟表明,改变回灌水帷幕的位置或水头高度都可以使得矿坑的排水量与帷幕的注水量相等,这就可以最大限度减少疏干水的外排水量。研究表明,疏干–回灌水帷幕可望成为我国露天煤矿区保水采煤的实用模式。      

A segregated flow scheme to control numerical dispersion for multi-component flow simulations

One of the challenges for reservoir simulation is numerical dispersion. For waterflooding applications the effect is controlled due to the self-sharpening nature of a Buckley–Leverett shock. However, for multi-component flow simulations, incorrect wavespeeds can develop leading to the excessive smearing of fronts because of the coupling of compositional dispersion with the fractional flow. Rather than implementing a higher-order discretization method, we propose a simple scheme based on segregation-in-flow within a gridblock to control numerical dispersion. We extend the method originally proposed for polymer flooding to augmented waterflooding simulations in general as well as simulations of miscible or near miscible gas injection. For compositional simulations of gas injection, this is done through a coupled limited-flash/upstream-exclusion assumption. To test the scheme, an in-house streamline simulator has been modified and validated for modeling low-salinity floods as well as ternary two-phase displacements. Simulation results presented with and without segregation demonstrate the potential of the approach as a heuristic method to control numerical dispersion in multi-component flow simulations.     

Numerical and experimental studies of pressure-controlled cavity expansion in completely decomposed granite soils of Hong Kong

Compaction grouting is the injection of a viscous grout into a soil under high pressure, which then densifies the surrounding soil by reducing void space. Laboratory and field tests of compaction grouting have been carried out. In this paper, a numerical model is used to simulate the compaction grouting process with the primary purpose of investigating relationships among various control parameters, such as injection pressure, void ratio and excess pore water pressure at various radial distances from the injection point. The compaction process is treated as a cavity expansion process in the numerical simulation. The soil is modelled with an elasto-plastic Mohr–Coulomb model using the commercial finite element program ABAQUS. In addition to numerical simulations, pressure-controlled cavity expansion laboratory tests were carried out on completely decomposed granite (CDG) soil specimens. Data collected from laboratory tests are compared with the finite element simulation to validate the finite element analyses. Factors that control the compaction process, such as the coefficient of earth pressure (K), initial void ratio, number of loading cycles and effective confining pressure, are explored in the numerical simulations.     

Finite element formulation for poroelastic problem with zero effective stress boundary condition

Zero effective stress boundary condition along with constant fluid flux is commonly encountered in geotechnical applications such as uncased borehole stability, fluid injection and production at an uncased borehole, hydraulic fracturing and sand production. This complex boundary condition introduces high nonlinearity in the numerical simulation. Conventional iterative methods such as Newton–Raphson method are required to solve this nonlinear problem iteratively, which involve huge computing time and also pose numerical difficulties on the convergence. To overcome this numerical difficulty and hence reduce the computing time, a novel numerical technique is proposed in this paper. Its performance is evaluated using a numerical example simulating fluid injection around an uncased borehole. Copyright © 2009 John Wiley & Sons, Ltd.     

疏松砂岩储层窜流通道平面分布规律研究

国内大部分高含水油田,注入水窜流现象较为严重,已成为影响油田高效开发的主要矛盾,对窜流通道形成后渗透率场的计算与描述显得尤为迫切和重要。目前窜流通道参数计算方法,只能计算得到注采井点的渗透率数值,不能计算得到注采井间的渗透率分布状况,并且目前窜流通道渗透率计算方法没有将室内实验结果、数值模拟结果和油藏工程方法结果有效结合起来。以国内某高含水油田为例,通过长期注水冲刷实验,模拟了窜流通道的动态形成过程,归纳出了储层岩石注水冲刷倍数与渗透率之间的关系式。在常规数值模拟计算方法的基础上,将渗透率变化规律模型引入到数值模拟计算中,建立了计算窜流通道平面分布的新方法。以海上某油田的地质模型为基础,计算得到了该油田目前状态的窜流通道渗透率平面二维分布规律。综合对比动态数据法、概率密度模型法、经验公式法和示踪剂测试结果,证明本方法得到的结果真实可靠。     

Numerical analysis of thermal fracturing in subsurface cold water injection by finite element methods

Integration of poromechanics and fracture mechanics plays an important role in understanding a series of thermal fracturing phenomena in subsurface porous media such as cold water flooding for enhanced oil recovery, produced‐water reinjection for waste disposal, cold water injection for geothermal energy extraction, and CO₂ injection for geosequestration. Thermal fracturing modeling is important to prevent the potential risks when fractures propagate into undesired zones, and it involves the coupling of heat transfer, mass transport, and stress change as well as the fracture propagation. Analytical method, finite element method, and finite difference method as well as boundary element method have been used to perform the thermal fracturing modeling considering different degrees and combinations of coupling. In this paper, extended finite element method is employed for the thermal fracturing modeling in a fully coupled fashion with remeshing avoided, and the stabilized finite element method is employed to account for the convection‐dominated heat transfer in the fracturing process with numerical oscillation circumvented. With the thermal fracturing model, a hypothetical numerical experiment on cold water injection into a deep warm aquifer is conducted. Results show that parameters such as injection rate, injection temperature, aquifer stiffness, and permeability can affect the fracture development in different ways and extended finite element method and stabilized finite element method provide effective tools for thermal fracturing simulation. Copyright © 2012 John Wiley & Sons, Ltd.     

Numerical simulation and optimization of CO2 sequestration in saline aquifers for vertical and horizontal well injection

With heightened concerns on CO₂ emissions from pulverized-coal (PC) power plants, there has been major emphasis in recent years on the development of safe and economical geological carbon sequestration (GCS) technology. Saline aquifers are considered very attractive for GCS because of their large storage capacity in U.S. and other parts of the world for long-term sequestration. However, uncertainties about storage efficiency as well as leakage risks remain major areas of concern that need to be addressed before the saline aquifers can be fully exploited for carbon sequestration. A genetic algorithm-based optimizer has been developed and coupled with the well-known multiphase numerical solver TOUGH2 to optimally examine various injection strategies for increasing the CO₂ storage efficiency as well as for reducing its plume migration. The optimal injection strategies for CO₂ injection employing a vertical injection well and a horizontal injection well are considered. To ensure the accuracy of the results, the combined hybrid numerical solver/optimizer code was validated by conducting simulations of three widely used benchmark problems employed by carbon sequestration researchers worldwide. The validated code is then employed to optimize the proposed water-alternating-gas injection scheme for CO₂ sequestration using both the vertical and the horizontal injection wells. The results suggest the potential benefits of CO₂ migration control and dissolution. The optimization capability of the hybrid code holds a great promise in studying a host of other problems in GCS, namely how to optimally enhance capillary trapping, accelerate the dissolution of CO₂ in water or brine, and immobilize the CO₂ plume.     

井间分溶示踪估计重非水相污染物残留量的影响因素数值分析

对于重非水相污染场地,传统的钻孔取样方式因经费和取样个数的限制,往往无法准确估计含水层中重非水相污染物的残留量。井间分溶示踪法费用较低,且不会对污染源区造成破坏,可适用于推估重非水相污染物残留量,但目前该方法的准确性尚未得到验证。本文基于数值模拟方法探讨了示踪剂注入速率、示踪剂类型以及污染源区结构等多因素对井间分溶示踪估计重非水相污染物残留量准确性的影响。数值算例结果表明:井间分溶示踪实验估计重非水相污染物残留量的误差小于15%;降低示踪剂的注入速率使得分溶性示踪剂在非水相和水相流体间作用更充分,能够提高井间分溶示踪估计重非水相污染物残留量的精度,实验中示踪剂注入速率由500 m3/d降为250 m3/d后,多种情景的平均精度由91.68%提高至93%;选取低分溶系数的示踪剂有利于提高示踪剂的回收率,从而提高推估精度,实验中示踪剂由2,2-二甲基-3-戊醇改为己醇后,平均精度由87.83%提高至96.85%;复杂的重非水相液体污染源区结构易于出现示踪剂绕流等现象,导致井间分溶示踪估计重非水相污染物残留量的误差增大,实验中含水层由均质变为非均质后,平均精度由93.03%变为91.65%。重非水相污染场地调查时,建议结合数值模拟方法选择适宜的示踪剂和示踪剂注入速率,以提高污染场地的刻画精度。     

Poro‐elasto‐plastic response of an unconsolidated formation confined with stiff seal rocks under radial injection

In this article, we evaluate geomechanics of fluid injection from a fully penetrating vertical well into an unconsolidated formation confined with stiff seal rocks. The coupled behavior of an isotropic, homogeneous sand layer is studied under injection pressures that are high enough to induce plasticity yet not fracturing. Propagation of the significant influence zone surrounding the injection borehole, quantified by the extent of the plastic domain in the elasto‐plastic model, is examined for the first time. First, a new fully coupled axisymmetric numerical model is developed. A comprehensive assessment is performed on pore pressures, stresses/strains, and failure planes during the entire transient period of an injection cycle. Results anticipate existence of five distinctive zones in terms of plasticity state: liquefaction at wellbore; two inner plastic domains surrounding the wellbore, where failure occurs along two planes and major principal stress is in vertical direction; remaining of the plastic domain, where formation fails along one plane and major principal stress is in radial direction; and a non‐plastic region. Failure mechanism at the wellbore is found to be shear followed by liquefaction. Next, a novel methodology is proposed based on which new weakly coupled poro‐elasto‐plastic analytical solutions are derived for all three stress/strain components. Unlike previous studies, extension of the plastic zone is obtained as a function of injection pressure, incorporating plasticity effects on the subsequent elastic domain. Solutions, proven to be a good approximation of numerical simulations, offer a huge advantage as the run time of coupled numerical simulations is considerably long. Copyright © 2016 John Wiley & Sons, Ltd.     

Real‐scale miscible grout injection experiment and performance of advection–dispersion–filtration model

A model was developed, to describe miscible grout propagation in a saturated deformable porous medium, based on Bear's statistical model with spatial volume averaging. In a previous paper, the model was first successfully confronted to one‐dimensional laboratory experiments. In the present paper, the numerical model is used to simulate practical grouting operation in a cylindrical injection model. The cylindrical injection model lends itself to study main flow and propagation character istics for a dispersed suspension‐type grout, under axisymmetric conditions close to real scale conditions. Comparison between numerical solutions and experimental results is essential to confirm the validity and accuracy of the proposed model from a phenomenological standpoint. The numerical model performances show that the underlying mathematical model constitutes a realistic predictive model reproducing most prominent features during injection of a suspension‐type grout into a deformable porous medium. The basic mechanism by which injected miscible grout permeates a soil mass is discussed in detail. Such a tool leads to quality control criteria for grouting on a theoretical basis, which complements existing criteria acquired through engineering practice. Copyright © 2001 John Wiley & Sons, Ltd.