1水文地球物理测井简介杨坤彪一、有关《地下水动力学》的基本问题㈠有关基本概念1地下水受重力作用在空隙介质中的运动称为渗透2将复杂的地下水运动的水流进行宏观概化,不考虑空隙、骨架占有空间和地下水实际运动途径,这种假想水流称渗流。在渗流场中取一过水断面,面积为ω其中空隙面积为ω′。通过ω′的实际地下水流量为Q,则地下水流实际速度(断面ω′上的均值)U=Q/ω′;而渗流速度(断面ω的均值)V=Q/ωV=ω′/ω×Uω′/ω=n为空隙率3渗流通过渗流中某点单位渗流途径长度上的水头损失称该点的水力坡度J=-(dH/dL)4直线渗透达西定律将均质砂装入直圆筒中做实验,如图A得出如下结果:12HHQKLQ—流量(m3/d)K—渗透系数(m/d)ω—直圆筒横截面或水断面积(m2)H1、H2—渗流方向相距L(米)的1、2两点的渗流水头值(m)(H1-H2)/L为水力坡度,则Q=KωJV=Q/ω=KJ表明渗流速度与水力坡度呈线性关系。12HHQKL为达西定律5渗透系数与渗透率达西公式中渗透系数K数值上相当于水力坡度为1的渗流速度。渗透系数与空隙介质的结构(n,d)以及水的性质(γ,μ)有关。232dKn忽略水的重率和粘滞性,可把渗透系数作为表征空隙介质本身透水性能的参数。用K0表示232nd0K称渗透率,只反映空隙介质本身的透水性。6紊流渗透(非直线渗透)采用哲才-克拉斯诺波里斯基公式表示紊流渗透TQKJTVKJKT为地下水呈紊流时空隙介质的渗透系数。当地下水内层流与紊流并存时的混合流用斯姆列盖尔给出的公式如图A21mCQKJ1mCVKJKC为混合流时空隙介质的渗透系数。m为流态指数(1<m<2)7裘布依微分方程裘布依研究地下水运动指出,绝大多数地下水流具有缓变流的特点。为了确定通过均质各向用性含水层任一过水段面ω的流量Q,先用达西公式的微分形式确定过dω的流量dQdHdQKddL积分dHdωdLQdQK-dH/dL表示微分过水断面dω上的水力坡度。dHQKdLdHVKdLdHdL表示具有缓变流特征时过水断面的水力坡度㈡地下水向井中运动的有关特征1含水层类型:潜水含水层或承压含水层2含水层的边界条件:周界(侧向)、顶、底界。按周界影响程度分无界含水层与有界含水层。是否从相邻含水层向抽水层补给分有越流和无越流含水层。3单井、干扰井:井之间未发生明显干扰,如同单独在含水层中工作视为单井;反之孔间发生干扰称为干扰井。4完整井、非完整井:根据井所揭露含水层的程度,揭露了含水层的全部厚度并在全部厚度范围内进水的为完整井;只一部分含水层厚度进水为非完整井。二、基本方法水文地球物理测井是通过测定井液电阻率来解决水文地质问题的。井液电阻率测井是在钻井中两个金属电极间的流体以Ω·m为单位的电阻率测量。能确定离子类型时可得出矿化度;有井温资料时可得出井液传导系数。井液电阻率测井可测定含水层的位置与厚度;确定不同水质溶液的补给或漏失部位;划分井中咸淡水界面;计算单一含水层的渗透速度和渗透系数;了解多层含水层的补给关系并得出涌水量。㈠、扩散法井液的矿化度与地下水矿化度不同时,大浓度溶液的盐离子向小浓度溶液扩散,使井液电阻率随时间发生变化。例如盐化井液后,每隔一定时间测一条井液电阻率曲线,测到异常反映明显为止。称为扩散法。通常当地下水渗透速度大于1米/昼夜时,扩散法应用效果较好,渗透速度小于0.1米/昼夜时,不宜采用扩散法。可配合采用提捞法。单一含水层⑴、定含水层厚度。如图1ρ0为盐化前的井液电阻率曲线,ρ1、ρ2、ρ3…为盐化后每隔一定时间测得的曲线,该曲线簇的宽度为含水层厚度。⑵、求含水层的渗透速度3盐化后的井液浓度随时间的变化与地下水的渗透速度V的关系如下:10201.81gCCdVlTCC其中C0—盐化前的井液浓度;C1—测ρ1时刻的井液浓度;C2—测ρ2时刻的井液浓度;d—钻孔半径;T—测ρ1的ρ2时间间隔。在一定温度下盐溶液浓度;与其电阻率成反比,因此101011.81gnnnnnndVltt其中ρ0—盐化前的井液电阻率;ρn、ρn+1—为tn和tn+1时测得的井液电阻率;由此每两条相邻曲线可算出一个渗透速度V,对于同一含水层可取各值的平均值作为该层的渗透速度值。为了提高精度,应用科尔变换公式求渗透速度:10110.78511nnndVTΔT=tn+1-tn为井液电阻率ρn变到ρn+1的时间间隔;d—钻孔半径。⑶、求渗透系数K根据渗透速度V选用合适的公式求渗透系数K层流运动用达西公式计算VKI(米/秒)水力坡度—I,能得出地下水流向时,根据流向线上两孔水位差与孔距求得水力坡度。如有等水位线图可得出水力坡度。紊流运动用哲才公式计算KVI混合流运动用斯沫来盖尔公式计算1.5VKI2、两层含水层有两个含水层时,两含水层间盐水柱由压力大的向压力小的含水层移动,盐水被含水层流出的淡水冲淡,井液电阻率曲线增高纵向移动。⑴、求补给关系和补给量井液扩散曲线如图2所示根据曲线判断为两层含水层,上层补给下层。井液4电阻率曲线纵向移动。不考虑其他影响时,其纵向位移速度代表孔中纵向水流速度,按下式计算补给量Q:2QdVd—钻孔半径;V⊥—纵向水流速度⑵、确定含水层位置测量时记录时间,每条曲线测量时间已知,纵向位移也可得出,将每个计算点的纵向流速代入上式算出补给量Q,绘制沿井轴变化的流量曲线(上图b)。轴向流量曲线的斜线ab为吸水层位置,根据ρ1、ρ2变化情况估计涌水层位置。⑶、求渗透系数K在两个含水层存在补给关系,但没有分层水位和影响半径的情况下,可用近似式计算渗透系数K26.25dVKHH—含水层厚度。裘布依和集哈尔德公式连解方程可求渗透系数K5.46g2g100HsKlKlqd其中H—含水层厚度;S—抬高水位;d—钻孔半径3、多层含水层图3为成都勘探设计院实测的井液电阻率曲线。盐化后井液电阻率曲线的特征呈现两个部分:井深37米以上及52米以下,上部分曲线随时间向下产生位移,说明含水层1向2补给;下部曲线则表明底部向上部位移,说明含水层3向4补给。分析曲线特征得出如下结果:存在四个含水层;含水层1和3为出水层、含水层3是承压水;含水层2和4为吸水层;含水层2和4之间为隔水层。4、井液电阻率测井的应用条件⑴、测井前要求对可能坍塌的井段用过滤花管保护井壁。⑵、测井前必须用清水洗井,使井壁孔隙畅通。⑶、宜用碳銨(NH4HCO3)进行盐化,浓度略小于2克/升,因其比重小于食盐,以便减少曲线5混乱交叉。㈡、注入法、提捞法经过盐化以后,如果渗透速度很小,施测时间很长时,可作注入法或提捞法。1、注入法注入法是井液盐化后,连续向孔中注入淡水,测出咸、淡水界面随时间的移动,直到淡水界面移动到不再向下移动、即含水层底界面为止。(图4)注入法可确定含水层位置;计算纵向水流速度V⊥;涌水量Q;单位涌水量q;渗透系数K等的近似值。hVtΔh—相邻曲线深度差Δt—相邻曲线时间间隔(tn+1-tn)2QdVd—钻孔半径QqSS—水位升高或降低的米数(孔中水位抬高或降低一米的涌水量)2、提捞法与注入法相似,而提捞法是降低水位,以便加速含水层水流入孔中的速度,使扩散加快以便确认存在含水层,从而确定其位置和厚度。实例方法说明:扩散法:是通过测量井液电阻率来了解井液扩散状况的一种方法。先以井液电阻计测出孔中淡水井液电阻率曲线,再在孔中加入电解质(食盐),使井液均匀盐化。然后每隔一定时间测一条井液电阻率曲线。若有含水层,地下水流入钻孔,会使井液淡化,电阻率增高得出曲线异常而可划分出含水层。而且由于盐化后的井液浓度随时间变化的快慢取决于含水层的渗透速度,因此可以根据盐化后的井液浓度随时间变化的关系来确定含水层的渗透速度。通常当地下水渗透速度大于1米/昼夜时,扩散法应用效果较好,渗透速度小于0.1米/昼夜时,不宜采用扩散法。可配合采用提捞法。本次采用的井液电阻计装置系数K=0.395,曲线以欧姆·米表示。注入法:是为了求得含水层的涌水量和渗透系数而采用的方法。注入法是从孔井口注入比井液电阻率低三倍以上、长度6约5米的盐水柱,同时向钻孔中注入与井液相同的淡水,保持抬高一定的水位,使盐水柱在一定水位静压力作用下而向下移动。以井液电阻计追踪盐水柱下移位置,并记录测定的时间,直到含水层时盐水柱消失。则根据盐水柱位移速度V、孔径d计算涌水量Q22dQV根据抬高水位值S,计算单位涌水量qq=Q/S同时可由裘布依和集哈尔德公式联解方程求渗透系数K。实例1:某矿区ZK16-13孔扩散法和注入法曲线见右图(1)求补给量Qin将扩散法曲线簇过中线的深度H、时间T及深度差ΔH、时间差ΔT列入下表,分别计算出各井段的速度得出其平均值。总共16段,总速度45.49(米/时)则V=45.49/16=2.843(m/h)=7.8972×10-4(m/s)钻孔直径为0.091m半径r=0.0455m则Qin=π×r2×V=5.136×10-6(m3/s)计算数据见下表求补给量Qm线号ρi深度(m)时间(h)速度(υ)值H差△H值T差△T418915.652001118.22.64.236209.59.521.33.13.067217.5823.82.53.2822911.53.94.12.89240116.42.54.41024339.43111256.513.512.434.51327518.518.96.52.8514280.55.521.32.42.2915290.5101.74.42.27163009.56.24.52.11717310109.233.33183199122.83.2119326.57.515.23.22.34203358.518.12.92.9321338321.23.10.97总段数16总计45.49(2)、求涌水量和单位涌水量将注入法曲线过中线两相邻同类界面的深度(顶界面H1底界面H2)深度差ΔH1、ΔH1和对应的时间(顶界面T1底界面T2),时间差ΔT1、ΔT1以及计算得出的速度v1、v2列于下表求涌水量和单位涌水量顶介面底介面深度(m)时间(h)深度(m)时间(h)线号ρi值差△H1值差△T1速度(υ1)值差△H2值差△T2速度(υ2)1124.59.763164.59.78215520.59.8200.057359.6185.5219.8370.057368.43165109.9580.13872.5200.5169.9770.14107.141942910.1730.215134.9205.5510.1780.20124.75202810.3680.19541.0217.51210.3730.19561.56211910.9750.60714.8222.5510.9800.6078.272372613.9983.0238.625734.514.0013.02111.482521514.8630.86517.3261414.8680.8674.692619.515.7810.91810.3265415.7860.9184.4102708.516.6800.8999.52811616.6860.917.8112851517.4980.81810.32961517.5030.81718.4123031818.6741.17615.3314.51818.6801.17715.31332117.519.6931.09617.2330.51619.6971.01715.714330.59.520.3720.6794.03365.520.3750.6888.1153365.521.1170.7457.4338221.1230.7482.716337121.8250.7081.4340221.8280.7052.8总计742.1671.3平均值v1=742.1/15=49.47(m/h)v2=671.3/1