《环境水文地球化学课件.ppt》由会员分享，可在线阅读，更多相关《环境水文地球化学课件.ppt（54页珍藏版）》请在三一办公上搜索。
1、吉林大学环境工程系2014年5月,污染水文地质学基本原理Principles of Contaminant Hydrogeology专题讲座,水文地质学系统结构,含水层性质,水动力性质,孔隙介质,裂隙介质,岩溶介质,双重介质,水文地球化学,环境水文地球化学,机械过滤吸附交替作用溶解、沉淀作用生物作用对流、弥散作用生物地球化学作用,污染质在地下环境中的主要衰减作用,第一讲 污染物的生物地球化学降解Biogeochemical Degradation of Pollutants,O2NO3-Mn4+Fe3+SO42-,（1）地下环境中的主要氧化还原物质 地下环境氧化还原系统中主要包括气体（O2、N
2、2、CH4、CO2）、溶解性组分（NO3-、NH4+、CH2O、Fe2+、Mn2+、SO42-、HS-、H+）、固体（FeOOH、MnO2、FeCO3、MnCO3）和存在于固体中的离子交换态组分（Fe2+、Mn2+、NH4+）。,1.污染物生物地球化学降解原理,1.污染物生物地球化学降解原理,1.污染物生物地球化学降解原理,1.污染物生物地球化学降解原理,（2）氧化还原带形成的原理生物地球化学原理 微生物在利用O2、NO3-、铁锰氧化物和氢氧化物、SO42-等最终电子受体并通过好氧呼吸、硝酸盐还原、铁锰还原、硫酸盐还原和产甲烷等代谢途径对有机污染物进行降解，在此过程中，由于最终电子受体争夺最终
3、电子能力大小的不同，导致氧化还原反应发生的难易程度也不同，所以在污染晕中形成了一系列不同的氧化还原环境，或者氧化还原梯度，也就是顺序氧化还原带，即产甲烷带、硫酸盐还原带、铁还原带、锰还原带、硝酸盐还原带、氧还原带等。其实质也是由于地下环境中的各种最终电子受体争夺电子能力的不同所致。热力学原理 O2、NO3-、Mn4+、Fe3+和SO42-发生氧化还原反应所需要能量不同，造成了在不同的氧化还原环境中发生不同的氧化还原反应的难易程度和次序的不同，从而导致污染晕在时间和空间上依次出现不同的氧化还原环境，也就是顺序氧化还原带。,1.污染物生物地球化学降解原理,2.氧化还原带中物质变化特征,3.氧化还原
4、带动态发展演化过程,第二讲 多相流Multiphase Flow,2.多相流,（1）基本概念非水相液体（Nonaqueous Phase Liquids,NAPL）不与水发生混合的液体称为NAPL;重非水相液体（Dense Nonaqueous Phase Liquids，DNAPL）比水重的非水相液体；轻非水相液体（Light Nonaqueous Phase Liquids，LNAPLs）比水轻的非水相液体自由相（Free Phase）在水动力条件下可流动的连续NAPL。残余相（Residual Phase）“困”在孔隙中，在水动力条件下不易流动的NAPL。,Mobile Phase NA
5、PL or(Free-Phase NAPL)is a continuous mass of NAPL that can flow under a hydraulic gradient,NAPL包气带和含水层中相分布示意图,Residual Phase NAPL is trapped in the pore spaces between the soil particles,and cannot be easily moved hydraulically,2.多相流,Can be estimated in the lab via the following:Where:rb=soil bulk
6、density g/cm3rn=NAPL bulk densityg/cm3n=porosityTPH=Total Petroleum Hydrocarbons mgNAPL/kgdry soil,2.多相流,饱和度：指流体充填体积占总孔隙体积的比例。剩余饱和度：因毛细作用使其固定不变时的饱和度,Partitioning Tracer Test for SN,Chemicals are introduced consisting of conservative and partitioning tracers(PTT)NAPL-Water Partition Coefficient KN.KN
7、=CN/Cw where CN is the conc.in NAPL and Cw is water conc.At steady state,the effect is to delay or retard the PT transport rate according to a retardation factorR=1+(KN SN)/(1-SN)=tp/tn,2.多相流,Tp is travel time for the Partitioning tracerTn is the travel time for the conservative tracer(Br-)Sorption
8、is assumed insignificantPT compounds often used include alcohols such as ethanol，n-pentanol(1.4)，n-hexanol，2,2 dimethyl-3-pentanol，n-heptanol,Measured BTC for PTT with bromideTime difference between different chemicals allows calculation of SN,Partitioning Tracer Test for SN,C/C0,Effluent Volume(L),
9、500100015002000,Bromiden-pentanol 2,2 DMP,2.多相流,表面张力：液体与其他物质接触（如固体、不想混合的液体或气体）将会产生界面能。这种能量是由于界面与内部分子吸引力不同而引起的，这种现象称为表面张力，定义为分子离开单位面积界面需要做的功，单位为达因/cm。,2.多相流,湿润性：一种或另一种液体会优先蔓延或湿润整个固体表面。界面角和表面张力之间的关系如下：,应在密度大的液体中测出。如果小于90，液体L将优先湿润固体；如果大于90，液体G将优先湿润固体。也就是两种液体争夺一固体表面，其中一种占主导地位，覆盖固体表面。,NAPL is Non wetting
10、 fluid,2.多相流,毛细压力（毛管压力）当两种不相混合的流体接触时，接触面发生弯曲。每相靠近界面的孔隙压力不一样，这种压力差就叫毛细管压力。毛细管压力在包气带为负值，当视为张力时为正值。,球形空气-水界面的曲率半径r、毛细压力Pc、表面张力之间的关系为：,毛细管压力与表面张力成正比，和曲率半径成反比。曲率半径取决于孔径以及所填充溶液的量。,毛细管压力是两种不想混溶的流体特性函数，即使在同一种介质中，毛细管压力也会由于水和NAPL比例的不同而不同。,2.多相流,相对渗透率 当两种不相混溶的流体同时流动时，共同占有了有效孔隙。在给定饱和度的情况下，流体的固有渗透率与岩石总的固有渗透率之比。,
11、2.多相流,束缚饱和度 在没有任何附加水流动时的含水量。在油-水两相流系统，若含水量小于束缚饱和度，油能够流动，水被毛细作用力控制；若含油量小于油的残余饱和度时，水能够流动，而油不能。,2.多相流,LNAPL Release to Vadose Zone,LNAPL at the Water Table,（2）NAPL迁移的概念模型,2.多相流,2.多相流,LNAPL在含水层中的迁移,2.多相流,DNAPL在含水层中的迁移,2.多相流,透镜体对DNAPL迁移的影响,DNAPL在裂隙中的迁移,2.多相流,Composite DNAPL Site,2.多相流,2.多相流,NAPL migratio
12、n in primary vs.secondary porosity,(A)NAPL migration through primary porosity.(B)NAPL migration through secondary porosity features(fractures,root holes,etc.).,A,B,2.多相流,LNAPL Moving Upgradient,2.多相流,Dissolving NAPL,2.多相流,Apparent LNAPL Thickness,Hf：LNAPL漂浮厚度Hw：井中LNAPL厚度w：水的密度LNAPL：LNAPL的密度,Thicknes
13、s in screened well casing only indicator of actual thickness in the unitThickness in the well may be 2 to 10 times larger than actual thickness of mobile LNAPL in the surrounding aquiferNo methods exist to improve these estimates,2.多相流,DNAPL Detection Problems,2.多相流,第三讲 饱和介质中的溶质运移Mass Transport in S
14、aturated Media,（1）浓度梯度作用下的迁移 水中的溶质从浓度高的地方向浓度低的地方迁移，这个过程叫做分子扩散，或扩散。流体扩散的量与浓度梯度成正比，称为Fick第一定律。F=-Dd（dC/dx）F：单位时间通过单位面积的溶质质量；Dd：扩散系数（L2/T）C：溶质浓度（M/L3)负号表示运动沿浓度减小的方向进行。不同物质的Dd值可以查到。对于浓度随时间的系统，适用于Fick第二定律，其一维表达形式为：,1.浓度梯度作用下的迁移,：为浓度随时间的变化。,在多孔介质中，由于离子在矿物颗粒周围扩散时需要经过更长的路径，因此，引入有效扩散系数D*。D*=Dd,Ci是t时刻距离污染源x处的
15、浓度；C0是初始浓度，erfc是余误差函数。,是一个与弯曲度有关的系数，通常小于1，通过扩散试验获得。浓度分布服从正态分布或高斯分布，还可用统计特征值来描述。,Fick第二定律一维表达式的解：,1.浓度梯度作用下的迁移,地下水运动中携带溶解性固体使之发生迁移，这个过程叫做对流迁移。,Vx为平均线速度；K为渗透系数；ne为有效孔隙度；dh/dl为水力梯度。,由于对流作用，一维质量通量Fx为：Fx=vxneC,2.对流迁移,地下水以大于和小于平均线速度的速率运动着。原因有：孔隙中心处的流速比孔隙壁处大；在直线距离相同时，多孔介质中有的流体质点会比其他质点沿着更长的路径运动；孔隙的大小不一，大孔流速
16、快。,若含有溶质的地下水以相同的速度运动，那么他们就会取代不含溶质的水，并在两种水体之间产生一个突变界面。然而，由于侵入的含有溶质的水不是以相同速度运动的，所以这两种水会沿着流动路径发生混合，这种混合就称为机械弥散，他导致在水流的推进边缘被稀释。沿着流动路径方向发生的混合称为纵向弥散；溶质锋面沿着垂直于流动方向发生的混合，称为横向弥散。,3.机械弥散,在地下水流动过程中，分子扩散过程不能从机械弥散中分离出来，将二者结合起来定义一个参数称为水动力弥散系数D，即DL=Lvi+D*DT=Tvi+D*DL：平行于水流运动方向的水动力弥散系数（纵向）；DT：垂直于水流运动主方向的水动力弥散系数（横向）；
17、L：纵向弥散度；T：横向弥散度。,4.水动力弥散,（4）对流-弥散方程推导T：横向弥散度。,2.饱和介质中的物质迁移,Fx,Fy,Fz,（1）吸附过程 吸附过程包括吸附、化学吸附、吸收和离子交换。吸附：包含溶质附着到固体表面的过程。阳离子交换：阳离子被吸附到带负电的粘土矿物附近并通过静电力附着在那里；化学吸附：溶质因化学反应而与沉积层、土壤或岩石表面结合。吸收：若含水层颗粒是多孔的，且使溶质能够通过扩散进入颗粒并被吸附到其表面上。,3.溶质转化、阻滞和衰减,（2）表面反应平衡1）线性等温吸附 溶质的浓度C和固体吸附的溶质的量C*之间存在一个直接的线性关系，即：C*=KdCC*：固体单位容重所吸
18、附溶质质量（mg/kg);C：吸附平衡时溶液中溶质的浓度（mg/L）；Kd为分配系数（L/kg），等于等温吸附线的斜率。,3.溶质转化、阻滞和衰减,整理后，得：,阻滞因子Rf由下式给出：,3.溶质转化、阻滞和衰减,如果地下水平均流速为vx，则浓度为初始浓度一半的溶质锋面的平均速度vc为：,线性等温吸附模型有两个缺陷：没有限制固体吸附溶质的量；点少时，易将C*-C曲线形式误解为直线。,2）Freundlich等温吸附,3.溶质转化、阻滞和衰减,对此式进行线性化：logC*=logK+NlogC,C*=KCN,整理后得：,阻滞因子Rff由下式给出：,若N1，锋面不断扩展；若N1，锋面会自锐；若N=
19、1，就是线性等温吸附,3.溶质转化、阻滞和衰减,是与结合能有关的吸附常数（L/mg）,是固体所能吸附的溶质的最大量（mg/kg）,Langmuir等温吸附也可表示为：,整理后得：,阻滞因子Rfl由下式给出：,3）Langmuir等温吸附,第四讲 地下水污染数值模拟Numerical Simulation of Groundwater Pollution,数值模拟的目的：测试一种假设，或提高对某一给定含水层系统的认识；了解物理、化学或生物过程；设计地下水污染修复系统；预测未来状况或模拟压力对地下水的影响；资源管理,4.地下水污染数值模拟,一、建模的基本信息（一）模型范围及边界 一般以相对独立的水
20、文地质单元为模拟区范围，便于准确刻画自然边界，提高精度。实际中，由于工作范围、资料精度以及模拟所限，模拟去可能不是一个完整的地下水系统，这种情况下就要人为划定模拟边界。1、对以河流河流或其它地表水体为边界的情况，可根据河流与地下水的关系，设置为一类（有直接水力联系，已知河流水位）、二类（河流与含水层脱节，并已知河流渗入补给量）或三类边界条件（河流下有弱透水层，已知河流水位）。2、以地下水分水岭、不透水断层为边界，可以定义为隔水边界，即二类边界，边界量设为零。3、以山区-平原交界线为边界的情况，一般可定义为二类边界，山区向平原的流入量可根据地下水流场或山区汇水面积等方法确定。4、对于人为边界，可
21、设置为通用水头边界5、泉、排水沟、内河河流边界等可直接处理成含水层的补给量。,4.地下水污染数值模拟,（二）水文地质结构 根据钻孔岩性、水文地质剖面图、地下水开采情况（深度）以及研究目的，划分模拟区含水岩组、弱透水层、和隔水层，再根据地下水流动特征，确定水文地质结构是单层结构还是多层结构。各层的顶板可用其标高等值线来表示，各层厚度可用其厚度等值线表示。（三）水文地质参数 水文地质参数主要包括两类：一类是为求含水层补排量所需的参数，如降水入渗系数、河流渗漏系数、蒸发系数等；另一类表明含水层介质性质参数，如水动力参数、有关的溶质运移参数、热运动参数等。常用的水文地质参数为含水层的渗透系数、承压含水
22、层的导水系数、潜水含水层的给水度、承压含水层的储水系数和释水系数。（四）初始流场 对于模拟非稳定流，需要提供模拟初始时刻的地下水流场，地下水水位观测点太少或分布不均，无法控制模拟区地下水流场。,4.地下水污染数值模拟,（五）源汇项通过水文地质调查、测量、计算以及其它手段获得地下水补排量（源汇项）。建模需要资料概括为下表：,4.地下水污染数值模拟,二、模型的选择与空间、时间的离散（一）空间离散 研究区用有限个网格或单元组成的离散域所代替，选择何种数值方法（有限差分法或有限单元法）将决定着空间离散特性（网格或单元）。有节点法和格点法两种剖分方法，目前多数采用格点法剖分。确定网格大小的总原则是：应考
23、虑各种分区界线，如水文地质单元，参数分区，行政分区、地表水体、断层和岩性界线。加密原则：（1）重点评价区或开采地段；（2）水文地质条件变化较大单元；（3）在地下水水位变化较大地段。,（二）模拟期的确定和时间离散（依据地下水动态变化特征确定模拟期）（三）时间步长（将模拟期离散成若干时段）,4.地下水污染数值模拟,三、模型的识别与检验,模型的识别是结合水文地质条件，通过调整模型的结构、参数、各种变量等，达到模拟的地下水状态（水头、浓度）与实测地下水状态最大程度的拟合，以使模型尽可能刻画实际水文地质条件，提高模型的仿真程度。模型识别原则：1、计算的地下水流场应与实际地下水流场基本一致；2、模拟计算期地下水水位应与实际地下水水位变化趋势一致；3、实际地下水补排差应接近于计算含水层储量变化量；4、识别后的水文地质参数、含水层结构和边界条件符合实际水文地质条件。,四、模型应用1、预测地下水分布状态和变化趋势，进行地下水资源评价和地下水环境影响评价；2、用于水资源管理。,4.地下水污染数值模拟,4.地下水污染数值模拟,
