一种确定多层含水层混合抽水等效井半径的方法技术

本发明专利技术公开了一种确定多层含水层混合抽水等效井半径的方法，是基于现场钻孔岩心层位分析，通过真实的多层含水层地下水流通过相应的钻孔过水断面来综合确定一个等效井半径的方法，即分别求出每一段含水层的过水断面面积，累加成一个整体的有效过水总面积，此时有效厚度为含水层的累加厚度，进而在钻孔整个深度上，确定出多层含水层抽水井的等效井半径。可解决多含水层联合抽水渗透系数求取或是在已知平均渗透系数时预报涌水量的问题；与工程中常用的等效体积法相比，该方法能真实地反映地下水运动情况和机理，比等效体积法更能反映实际，更为精确有效。

【技术实现步骤摘要】
一种确定多层含水层混合抽水等效井半径的方法
本专利技术属于地下水资源开发利用领域，具体涉及一种确定多层含水层混合抽水等效井半径的方法。
技术介绍
现有含水层参数都是通过单一含水层进行抽水试验来获取水文地质参数，然而在实际生产中，如水源勘探孔，由于深度较大，常常深达几百米甚至上千米，穿越了多个含水层(段)，此时，进行抽水试验时，若一个含水层一个含水层的进行，则当主含水层抽水时，需要对其它含水层变径止水，工艺复杂，废时废钱也废力。而生产实际中为扩大涌水量常采用多含水层联合抽水，钻孔在不同的含水层根据需要采取变径方式，然后，多级离心泵下至一定深度进行大流量联合抽水，此种情况下，多含水层流量和水位关系如何，如何运用已有的地下水井流解析公式来进行水资源预测和评价，如何确定多含水层、不同井半径条件下的等效井半径则成了非常困难的一件工作。为了简化，目前生产单位常用等效体积法来求算等效井半径，就是把各含水层变径的钻孔空间体积换算成同一孔径相同钻孔深度的体积，换算得到的同一孔径即当作等效井半径，以此来预报多含水层联合抽水流量。显然，这种方法的缺点是简单的把钻孔空间都当作均一进水空间，简单的把参与计算的所有岩层都当作含水层，不能正确反映地下含水层真实层位状况。同时半径在体积表达式中为平方项从而相对降低了含水层厚度的影响，不能正确反映地下水在含水层中运动情况和动态变化规律的基本事实。事实上，其计算结果在精度上也无法保证。
技术实现思路
鉴于上述现有技术存在的不足，本专利技术提供了一种确定多层含水层混合抽水等效井半径的方法，以解决现有技术中求算等效井半径精确度低，导致抽水试验涌水量预报和水文地质参数求取精确度不高以及对地下水运动情况和机理不能正确真实反映的问题。为了解决上述技术问题，本专利技术采取如下技术方案：一种确定多层含水层混合抽水等效井半径的方法，抽水钻孔穿越多层含水层，且经过若干变径，根据各含水层与钻孔接触的面积称为过水断面，分别求出每一层含水层的过水断面面积，累加成一个整体的有效过水断面总面积，该有效过水断面总面积除以含水层整体厚度与2π的乘积，即得到一个等效的井半径，具体包括如下步骤：1)确定变径抽水钻孔的每段变径半径ri和每段岩层中含水层厚度Mi；2)确定每一层含水层与钻孔接触的面积，即每一层含水层过水断面面积Si；3)确定是否需要封堵某些或者某个含水层；4)将变径抽水钻孔通过的各个过水断面面积Si相加求得有效过水断面的总面积∑Si，若需要封堵某些或者某个含水层时，总面积Si中要扣除被封堵含水层相对应的过水断面面积Si；5)将上述参数代入公式r等＝∑Si÷[2π∑Mi]中，计算出等效井半径r等。进一步地，步骤1)中所述每段变径半径ri通过钻孔管径D来确定(为其一半)，所述含水层厚度Mi通过钻孔岩心分析或用钻孔测井方法得出。进一步地，步骤2)中每一层含水层过水断面面积Si通过通式Si＝2πriMi计算得出，其中i表示第i层含水层，i＝1，2，3，……n。进一步地，步骤3)中若需要封堵某些或者某个含水层时是通过套管或止水材料封堵。进一步地，步骤4)中若不封堵含水层，i＝1，2，3，……n，∑Si＝S1+S2+……+Sn,，若封堵第j或k层含水层，则扣除相应Sj或Sk，余下各项相加得出∑Si。进一步地，步骤5)中∑Mi为若干段岩层中含水层总厚度，i＝1，2，3，……n，∑Mi＝M1+M2+……+Mn，若需要封堵封堵某些或者某个含水层，则应扣除相应含水层的厚度Mi。进一步地，所述ri，在当i＝1，2，3，……n，r1、r2、……rn均相等时，则为不变径抽水钻孔同样适用等效井半径r等公式。与现有技术相比，本专利技术有益效果为：该方法能够根据钻孔实际穿过若干含水层的条件来综合确定较为精确的等效井半径，使抽水试验涌水量预报和水文地质参数求取较为准确真实，并能正确真实反映地下水运动情况和机理；而且该方法将为稳定流条件下已知平均渗透系数来预报多层含水层涌水量或已知涌水量条件下来计算平均渗透系数等提供精确的科学计算依据。附图说明图1为本专利技术钻孔结构示意图。具体实施方式下面结合附图和实例对本专利技术进行进一步说明。如图1所示，图中有若干个含水层，首先依钻孔岩心分析或用钻孔测井方法得到含水层的厚度，假设含水层厚度分别为M1、M2、……Mi、……Mn，其对应每一层含水层的过水断面面积分别为S1、S2、……Si、……Sn，含水层与相邻隔水层的厚度分别为M1'、M2'、……Mi'、……Mn'，其对应的体积相应为V1、V2、……Vi、……Vn，其中r1、r2、……ri、……rn分别为钻孔的半径，它们可相等，也可不等，不等则说明发生了变径。为了更好理解，采用将现有技术体积等效法和本专利技术实施例方法对比着叙述，更体现本专利技术的有益效果。现有体积等效法计算等效井半径方法：则依体积等效有，亦即从而求得采用确定多层含水层混合抽水等效井半径(有效过水断面法)的方法计算：S1＝2πr1M1，S2＝2πr2M2，……，Si＝2πriMi，……，Sn＝2πrnMn，∑Si＝S1+S2+……+Sn，∑Mi＝M1+M2+……+Mn，从而求得r等＝∑Si÷[2π∑Mi]，(i＝1,2,3，……n)由以上公式计算出等效井半径，进一步可解决多含水层联合抽水试验参数求取的问题。下面以淄博博山区谢家店XK1钻孔现场稳定流抽水试验为例，结合具体的实验数据来验证该方法的可实施性和效果。该钻孔为三次钻孔变径，其位于谢家店村东大桥北300米河东，井深353.15m，共分为四级口径：开孔口径426mm，变径深度在0～16.3m；二级口径377mm，变径深度在16.3～105.9m；三级口径325mm，变径深度在105.9～250.2m；四级口径273mm，变径深度在250.2～353.15m。揭露地层O2b、O2d、∈4O1s、其中包含12个含水层，13个隔水层，含水层总厚度为50.29m。因套管下到第二层，即抽水从第三段开始，计算有效过水断面总面积和含水层整体厚度不包括第一、二段。先采用现有技术计算方法，即生产单位常用的等效体积法计算结果如下：表1体积等效法计算结果∑V＝V3+V4＝18.0m3M总＝M3'+M4'＝247.25m注：M3'、M4'分别为第3、4段岩层的厚度(包括含水层和隔水层)；r3、r4分别为第3、4段钻孔的半径。通过等效体积法计算求得等效井半径r体等为152.23mm，再根据等效井半径和XK1特征钻孔的抽水资料运用流场稳定时的如下两个地下水动力学公式:式中：K—渗透系数(m/d)，Q—稳定涌水量(m3/d)M—含水层总厚度(m)，S—抽水孔稳定水位降深(m)r—抽水孔等效井半径(m)，R—抽水影响半径(m)运用迭代法求得渗透系数K为4.510m/d。相应地我们再采取本专利技术提供的一种确定多层含水层混合抽水等效井半径的方法，包括如下步骤：1)确定变径抽水钻孔的每段变径半径ri和每段岩层中含水层厚度Mi，所述每段变径半径ri通过钻孔的钻的半径确定，所述含水层厚度Mi通过钻孔岩心分析或用钻孔测井方法得出；2)确定每一层含水层与钻孔接触的面积，即每一层含水层过水断面面积Si，每一层含水层过水断面面积Si通过通式Si＝2πriMi计算得出，其中i表示第i层含水层,i＝1,2,3，……n本文档来自技高网...

【技术保护点】
1.一种确定多层含水层混合抽水等效井半径的方法，抽水钻孔穿越多层含水层，且经过若干变径，分别求出每一层含水层的过水断面面积，累加成一个整体的有效过水断面总面积，该有效过水断面总面积除以含水层整体厚度与2π的乘积，即得到一个等效的井半径，所述过水断面是指所述含水层与钻孔接触的面积，包括如下步骤：1）确定变径抽水钻孔的每段变径半径ri和每段岩层中含水层厚度Mi；2）确定每一层含水层与钻孔接触的面积，即每一层含水层过水断面面积Si；3）确定是否需要封堵某些或者某个含水层；4）将变径抽水钻孔通过的各个过水断面面积Si相加求得有效过水断面的总面积 ∑Si，若需要封堵某些或者某个含水层时，总面积Si中要扣除被封堵含水层相对应的过水断面面积Si；5）将上述参数代入公式r等=∑Si÷[2π∑Mi]中，计算出等效井半径r等。

【技术特征摘要】
1.一种确定多层含水层混合抽水等效井半径的方法，抽水钻孔穿越多层含水层，且经过若干变径，分别求出每一层含水层的过水断面面积，累加成一个整体的有效过水断面总面积，该有效过水断面总面积除以含水层整体厚度与2π的乘积，即得到一个等效的井半径，所述过水断面是指所述含水层与钻孔接触的面积，包括如下步骤：1）确定变径抽水钻孔的每段变径半径ri和每段岩层中含水层厚度Mi；2）确定每一层含水层与钻孔接触的面积，即每一层含水层过水断面面积Si；3）确定是否需要封堵某些或者某个含水层；4）将变径抽水钻孔通过的各个过水断面面积Si相加求得有效过水断面的总面积∑Si，若需要封堵某些或者某个含水层时，总面积Si中要扣除被封堵含水层相对应的过水断面面积Si；5）将上述参数代入公式r等=∑Si÷[2π∑Mi]中，计算出等效井半径r等。2.根据权利要求1所述的一种确定多层含水层混合抽水等效井半径的方法，其特征在于：步骤1）中所述每段变径半径ri通过所施工钻孔管径确定，所述含水层厚度Mi通过钻孔岩心分析或用钻孔测井方法得出。3.根据权利要求1所述的一种确定多层含水层混合抽水等效井半径的方法...

【专利技术属性】
技术研发人员：齐跃明，孙箐彬，许天阳，袁冬梅，杨雅琪，游京，裴毅峰，叶明达，林旺辉，蔚波，张亮，包萌，马仪鹏，
申请(专利权)人：中国矿业大学，
类型：发明
国别省市：江苏,32