共查询到20条相似文献，搜索用时 509 毫秒
1.
对于模拟示踪试验的溶质运移,传统方法采用细网格离散小尺度表皮层,存在网格数量多和计算时间长等问题。由于表皮层外缘的参数变化剧烈,即使采用细网格,仍然存在显著误差。因此发展基于完整井的新示踪试验溶质运移模型,提出新的瞬态Robin边界条件反映表皮层吸附/解吸溶质的影响,实现表皮层的无网格。应用拉普拉斯变换导出模型的解析解,通过有限元法建立基于非完整井的数值解。结果显示新瞬态Robin边界条件精准反映负表皮层的影响。对于求解反演问题,表皮层宽度的估计范围由0.45～0.54 m缩至0.47～0.48 m,纵向弥散度由0.6～10.0 m缩至6.4～7.7 m,参数估计的可靠度显著提升。表皮层的无网格节省97%的数值解计算时间。 相似文献
2.
地下水污染问题日益严重,研究溶质运移的弥散理论开始应用于实际问题。建立地下水溶质运移模型,对地下水中污染物的运移及发展趋势进行准确预测,是对地下水进行保护、对地下水污染进行控制的基础。而弥散参数的确定则是地下水溶质运移模型建立的关键环节之一,直接影响着模型预测结果的精度和准确性。 对西宁市贵德县地下水污染的水质运移规律进行分析,在贵德县河滨公园林场采用径向收敛流水动力弥散理论方法进行了第四系含水层现场弥散试验,计算了试验场地潜水含水层的弥散度,获得纵向弥散度(aL)为0.843~0.998 cm,横向弥散度(aT)经验推断值为0.17~0.20 cm,为进一步建立该地区的地下水溶质运移模型、预测地下水污染的发展趋势和评价该地区地下水环境质量提供了数据参考。 相似文献
3.
二维分数阶对流-弥散方程的数值解 总被引：1，自引：0，他引：1
对二维时间分数阶对流-弥散方程和二维空间分数阶对流-弥散方程分别建立了差分格式,实现了对其的数值求解。针对理想算例进行计算求解,分析了时间和空间分数阶阶数取不同值时的扩散变化规律,验证了各自所描述的时间相关性与空间相关性。同时与传统的二维整数阶对流-弥散方程的求解结果作了对比。当时间和空间分数阶阶数α与γ分别取整数时,二维时间分数阶对流-弥散方程和二维空间分数阶对流-弥散方程都与传统二维整数阶对流-弥散方程的计算结果相同,说明提出的对二维分数阶对流-弥散方程的数值求解方法是可行的。其结果对地下水溶质运移的进一步研究提供了有效的手段。 相似文献
4.
地下水中污染物运移的数值模拟方法一直是学界的研究热点问题.而如何减少与消除对流-弥散方程数值解中浓度陡锋面附近的数值振荡与数值弥散,更是研究的前沿与难点.提出了一种地下水溶质运移数值模拟中减少数值弥散的新方法.该方法的核心思想是在水动力弥散系数上加上一个数值弥散估算值,得到一个修正弥散系数,用其替代方程中有明确物理意义的水动力弥散系数进行计算.并提出了一个参数——数值弥散因子(μNDF),可以根据研究需要进行参数分区并适当调节该因子的大小,从而达到控制数值振荡,减小数值弥散的目的.从一维到二维的多个数值算例的模拟计算结果表明,该方法能在消除数值振荡的基础上,较好地减少数值弥散,达到满意的精度. 相似文献
5.
为研究城市小型河流中污染物的物理迁移过程规律,分析基流条件下流动水体与暂态存储区之间的滞留交互作用,采用溴化锂(LiBr)作为保守性示踪剂进行野外现场示踪试验,结合一维溶质运移存储模型(One-dimensional Transport with Inflow and Storage model, OTIS)定量解析潜流交换特性,估算纵向弥散系数(D)、潜流交换面积(As)、主河道断面面积(A)和潜流交换系数(α).模型度量指标DaI值和均方根误差值结果表征参数模拟结果可靠性高,拟合效果理想.由泵入点O至下游1 300 m设置的A、B、C、D 4处监测点的模拟结果表明,水文参数D、As、A和α均随水文条件而变,OB河段(0~600 m)潜流交换能力较弱,主要以对流弥散过程为主;BD河段(600~1 300 m)具有较强的暂态存储能力,对溶质的滞留时间长;BC(600~1 000 m)和CD(1 000~1 300 m)河段交换系数分别为(3.42×10-6±0.65×10-6)s-1和(2.87×10-6±0.81×10-6 )s-1;河段BC存在2.2×10-5m3/(s·m)的侧向补给流量.4个河段对比发现,城市河流渠道化、河床沉积物贫瘠等特征导致潜流交换能力弱化. 相似文献
6.
基于Fredlund等提出的非饱和土轴对称固结理论,采用傅里叶级数展开法和拉普拉斯变换法提出了等应变条件下分段循环荷载作用下同时考虑径向(r方向)和竖向(z方向)渗透的非饱和土轴对称固结半解析解。基于等应变假设和径向边界条件,给出了等应变条件下的超孔隙气压力ua和超孔隙水压力uw的控制方程。然后利用级数展开与Laplace变换求解并给出了超孔隙压力及平均固结度的半解析解。此外,通过退化法以及有限差分法验证了所提解和相应结果的正确性。引入3种具体的分段循环荷载,研究了渗透系数之比ka/kw、模型尺寸N和荷载特征参数a对非饱和土轴对称固结特性的影响,所得结论更加符合实际,可为非饱和土地基处理提供理论指导。 相似文献
7.
为确定研究区含水层的弥散系数,进一步开展咸水体运移的数值模拟,根据研究区的地质背景和水文地质条件,针对含水层中咸水体分布格局设计并开展了示踪试验。利用试验获得的数据初步分析了咸水体垂向和横向运移特征,提出了依赖于单尺度示踪试验结果推求多尺度下弥散系数的内推-外延法,据此计算了试验尺度上的弥散系数。示踪试验结果直观地反映了咸水体对抽水井中氯离子质量浓度的影响,确定了试验尺度上的横向弥散度为7.42 m,垂向弥散度为41.4 m。基于弥散方程的解析解,采用半程拟合的方法进行尺度效应的分析,确立了横向及垂向弥散度与试验尺度关系分别为y=0.007 5·e0.129 3x和y=0.102 9e0.049 2x。根据上述定量关系可进行同类地区相似条件下的弥散系数计算和尺度效应分析,从而为进一步开展咸淡水交互地带水盐运移的模拟和预测提供参数基础。 相似文献
8.
采用图像分析法确定溶质浓度,以硫酸铜和EDTA二钠作为双分子反应物,在多孔介质模型中开展了不同粒径(1.52.0,2.53.0,3.54.0 mm)和流量(1.0,1.5,2.0 mL/s)下反应性溶质运移实验,探讨了应用不完全混合的对流弥散模型(IM-ADRE)对双分子反应溶质运移的模拟和预测,并进行了参数敏感性分析。结果表明:图像分析法可准确获取多孔介质中显色反应性溶质的浓度,灰度值与浓度的决定系数R2大于0.96;用IM-ADRE模型能够准确预测双分子反应性溶质硫酸铜和EDTA二钠在3种不同多孔介质中的运移过程,误差低于3.71%;实验条件的改变对IM-ADRE模型参数D、m和β0的影响显著,说明模型参数依赖于环境条件,其变化规律需要根据实际环境条件进一步率定,便于IM-ADRE模型的进一步推广应用。 相似文献
9.
我国西部多数煤矿区矿井水往往含有较多的生产污水,导致其氨氮含量较高,利用煤矿采空区充填矸石对矿井水进行预处理已成为矿井水及矸石资源化利用的重要措施。以神华神东煤炭集团公司保德矿矸石为充填介质、以矿井水中的氨氮为研究对象,通过模拟采空区的水文地质环境,开展柱模拟实验,结合CXTFIT 2.1软件,研究氨氮在充填矸石中的运移机制,这对于评价煤矿采空区充填矸石对矿井水的预处理效果具有重要意义。结果表明:达西流速3.12 cm/h、25℃时,Cl?在充填矸石中的运移可用CDE模型较好地表征(r2=0.999),较高的纵向弥散系数(D)、溶质分子扩散系数(Df)、机械弥散系数(Dh)和弥散度(λ)值与充填矸石较大粒径、较长运移距离有关,淋出液中总溶解性固体(TDS)和Cl?的质量浓度呈显著线性关系;双点位吸附溶质运移模型能够较好地表征氨氮的运移过程,氨氮的阻滞系数R为23.79,其在矸石上的平衡吸附点位仅占总吸附点位的46%,一级动力学吸附速率常数α为3.5×10?4 h?1;氨氮的孔隙流速及水动力弥散系数远低于Cl?的数值,这主要与高含量的黏土矿物高岭石对氨氮的吸附有关;在649 h的氨氮运移过程中,亚硝酸盐和硝酸盐氮以及pH值均无明显变化,实验中后期TDS稳定,这进一步证实氨氮在模拟的矸石柱运移过程中以对流、弥散迁移及吸附为主,生物转化作用可以忽略。研究结果可为评价我国煤矿采空区地下水库的水质预处理技术提供重要理论依据。 相似文献
10.
岩溶管道溶质运移的尺度效应研究对穿透曲线的正确解译非常重要,但目前针对单一弯曲管道中溶质运移尺度效应的研究仍比较缺乏。文章将岩溶管道和溶潭分别概化为透明软管和水箱,基于前期建立的水箱-管道系统(简称“管道系统”),在水箱下游设置不同长度的弯曲管道,通过示踪试验研究管道运移尺度对穿透曲线的影响,并采用暂时存储模型模拟试验曲线。结果表明：(1)随着水箱下游管道长度的增加,峰值质量浓度逐渐缓慢降低,穿透曲线上升段斜率无明显变化,穿透曲线拖尾逐渐缩短,表明运移管道长度增加对溶质运移的影响大于下游管道弯曲;(2)穿透曲线偏度系数、后段溶质羽穿透时间和溶质羽穿透时间与管道系统长度呈良好的负线性相关关系(R2>0.96);(3)当对称和不对称水箱管道系统长度分别增加至154.5 m和164.3 m时,偏度系数接近0,穿透曲线分布接近对称;(4)弥散系数、存储区截面积和交换系数与管道系统长度呈良好的负线性相关关系,当对称和不对称水箱管道系统长度分别增加至159.9 m和178.1 m时,存储区截面积接近0,水箱导致的溶质运移滞后效应基本消失。研究结果对野外岩溶管道穿透曲线的... 相似文献
11.
内源瞬态荷载作用下圆柱形孔洞的动力响应解答是土动力学的经典问题之一。已有研究大都假设孔洞周围土体为理想弹性介质或完全饱和多孔介质。然而,实际工程中不存在完全弹性和完全饱和土体。分别视衬砌结构和周围土体为弹性材料和准饱和多孔介质（饱和度 95%）,根据牛顿第二定律、达西定律和Biot波动理论推导出准饱和土体的控制方程。根据边界条件导出衬砌和土体的位移、应力和孔隙压力的Laplace变换空间的解答。利用反Laplace变换数值计算方法,将解答转换为时域解。分析了饱和度对衬砌位移、应力和孔压的影响,结果表明,当95% 99%时,饱和度对径向位移和切向应力的影响较小;99% 100%时,饱和度对径向位移和切向应力的影响较大;但饱和度对孔隙压力的影响远大于对径向位移和切向应力的影响。得出位移、应力和孔压沿径向的衰减规律,当95% 99%时,饱和度对径向位移和切向应力沿径向衰减影响较小,99% 100%时,饱和度对径向位移和切向应力沿径向衰减影响较大,但饱和度对孔压沿径向的衰减影响远大于对径向位移和切向应力沿径向的衰减。 相似文献
12.
针对非等温分布条件下由土工膜(geomembrane,简称GMB)+土工复合膨润土衬垫(geosynthetic clay liner,简称GCL)+压实黏土衬垫(compacted clay liner,简称CCL)组成的3层复合衬垫中有机污染物一维瞬态运移问题,考虑了对流、扩散、机械弥散、吸附、降解和热扩散等因素,建立了相应的数学模型。该模型可考虑扩散系数和渗透系数等参数随温度变化。利用有限差分法,获得了该模型的数值解。通过将所建模型计算结果分别与试验结果、已有解析模型计算结果和COMSOL软件计算结果进行对比,验证了所建模型的正确性。基于所定义的击穿时间tb,以甲苯作为代表性有机污染物分析和讨论了不同因素对运移行为的影响。结果表明：(1)非等温分布条件会使渗透系数、扩散系数和线性吸附系数等参数发生变化,其中渗透系数和扩散系数的变化会加快运移过程,尤其是扩散系数,而线性吸附系数的变化会减慢运移速率;(2)复合衬垫上部和底部温度差ΔT的增大会使击穿时间tb减小,底部运移通量Jb增大;当ΔT为0、10、20、30、40 ... 相似文献
13.
川西南震旦系储集层有机包裹体在油气运移研究中的应用 总被引：1，自引：0，他引：1
有机包裹体在油气成因、运移和演化研究中具有重要的意义.川西南震旦系灯影组储层中次生成因孔洞缝发育, 经历了多期白云石和硅质充填.对不同期次充填矿物中的有机包裹体特征及荧光进行了观察, 测定了各期盐水包裹体均一温度, 其主峰分布在120~150℃、160~190℃和200~210℃; 并采用真空破碎法对各期有机包裹体的化学成分及碳同位素组成进行了分析, 其化学成分以CH4和CO2为主, 各期中的w(CH4)含量分别为10.09%~22.08%、27.06%~40.40%和55.33%~74.01%, 相应的w(CO2)/w(CH4)值为3.37~5.87, 1.20~2.20和0.26~0.63.CH4和CO2的碳同位素组成分别为(-50.7~-28.4)×10-3和(-18.3~-3.6)×10-3.据此追述了研究区油气的成因性质、油气的演化和运移.结果表明, 研究区震旦系灯影组储层的天然气为油型气; CO2主要为有机成因, 有少量为无机成因; 研究区存在有3期以上的油气运移, 且各期油气演化程度不同. 相似文献
14.
为评估土工膜与土-膨润土组成的复合隔离墙对有机污染物的防污性能,建立了污染物在源区域发生降解扩散时透过复合隔离墙的一维瞬态扩散模型。通过拉普拉斯变换和Talbot数值反演获得了模型解。第3类入流边界条件下,考虑源区的扩散和降解时,第1类隔离墙(土工膜/土-膨润土)和第2类隔离墙(土-膨润土/土工膜/土-膨润土)100 a的击穿浓度分别较源浓度恒定时减少了59%和53%。由于假设土工膜渗透系数高于10–12 m/s,污染物在复合隔离墙中的对流作用不能忽略,因此,采用第2类复合隔离墙使土-膨润土能更好地发挥主导隔离作用。将土工膜的渗透系数从10–10 m/s减少到10–16 m/s,第1类复合隔离墙的击穿时间从26 a增加到188 a,第2类复合隔离墙的击穿时间从32 a增加到81 a。利用抽水井调整墙体内外为负水头差及提高源区的污染物降解能力,可使隔离墙的防污性能得到较大提升。 相似文献
15.
基于Dakshanamurthy和Fredlund提出的二维非饱和土固结理论,利用Fourier正弦级数展开、Laplace变换,分别给出了分段循环荷载作用下二维非饱和土固结问题的超孔隙气压力、超孔隙水压力和沉降的半解析解,并应用退化法验证了本文所得半解析解的正确性。然后,结合3种具体的荷载形式,分析了分段循环荷载作用下气相与液相渗透系数之比(ka/kw)、水平方向与竖直方向渗透系数之比(kx/kz)和荷载特征参数(a)对二维非饱和土固结特性的影响。结果表明:ka/kw和kx/kz的增大均会加速固结沉降进程;荷载特征参数越大,沉降发展越早,沉降值越小;二维非饱和土固结特性受分段循环荷载作用影响明显。因此,在实际施工过程中改变施工速度和设置径向排水装置可有效控制二维条件下非饱和土体的固结过程,该研究成果可为非饱和土地基的设计和施工提供重要理论依据。 相似文献
16.
合理布局抽水井和回灌井是地下水源热泵系统长期有效运行的关键因素。以郑州市郑东新区为例,利用HST3D软件建立水热运移数值模型,优化设计地下水源热泵系统抽灌水井布局,预测地下水源热泵系统长期运行后对含水层的水热影响。结果表明:介质比热容及渗透率分别对含水层温度及水位影响显著,是较敏感的参数。方案3(3个回灌井垂直天然流向分布且位于抽水井下游)为最佳布井方式。抽灌量900,1 200,1 500及2 000m3/d所对应的合理布井间距分别为50,65,70及75m。相应布井方案的水源热泵系统运行20a,对含水层温度场的影响仅限于200m×200m的范围,抽水井温度变化小于1℃。 相似文献
17.
以模拟农田尺度降雨入渗—重分布条件下阿特拉津(atrazine)在非饱和土壤中的淋溶动态为数值算例,通过在北京东郊的一块27m×27m的农田采集100个土壤样品,测定其主要理化性质,间接计算获得各土壤样品的van Genuchten型水力学参数、溶质运移的弥散度及atrazine在土壤中吸附的分配系数,进而基于柱模型假设,运用HY DRUS 1D软件得到不同时刻采样区内土壤水负压和atrazine浓度的平均分布。在此基础上,重点探讨了应用幂平均算法计算不同的等价有效参数对农田尺度下土壤水负压和atrazine浓度平均分布的数值模拟精度的影响。就所设计的降雨入渗—重分布条件下atrazine淋溶的算例而言,土壤水力学参数、溶质运移参数和吸附特性参数取几何平均或调和平均作为有效参数值,空间平面上呈正态分布的土壤水负压的初值采用算术平均值作为流场模拟的初始条件,可以获得满意的数值模拟结果。 相似文献
18.
三峡库区在高程145~175 m之间周期性水位调控,促使库区灰岩岸坡劣化损伤加剧、稳定性衰减。本文结合巫峡段顺层灰岩岸坡“滑移-弯曲”破坏实例,基于弹塑性板翘曲模型,考虑岩体动态劣化概念,结合广义H-B准则中GSI(t)岩体参数动态指标,推导得临界挠曲段平衡方程。与传统静态计算方法相比,适用性更广,可以考虑岸坡随劣化时间t动态演化过程,更加符合工程实际。结果表明:随着劣化进行临界挠曲段lcr逐步降低,岸坡稳定状态变差,下滑推力F逐步增大,抗弯刚度K逐渐减小,也正是这两方面因素导致lcr快速降低。由于岩体劣化函数GSI(t)为指数劣化形式,故lcr呈现出快速降低并趋缓直至逼近一定值,该区段为岸坡防治的“关键区段”;随着岩层倾角α增大,lcr与相对板长比RPL快速降低,随着岩层厚度hi增大,lcr与RPL缓慢变大,α对lcr与RPL的影响明显大于hi。且岸坡lcr快速下降角α大约在40°~60°,可作为类似顺向岸坡是否容易滑移-弯曲失稳的地质判据。本文推导的考虑岩体劣化的临界平衡方程,以三峡库区青石6号岸坡为例进行了计算、论证。可为其他类似顺层岸坡损伤演化计算提供借鉴,不同之处在于确定具体岸坡损伤劣化的GSI(t)函数式。 相似文献
19.
详细分析不同评价指标下的van Genuchten方程参数敏感性,对于深入认识参数物理意义、合理确定参数值、提高降水入渗数值模拟准确性具有重要意义.用HYDRUS软件建立一维变饱和水分运移模型,以实际补给量、补给过程曲线和极限蒸发深度作为评价指标,采用单因素扰动分析法对van Genuchten方程各参数敏感性进行理论分析,进而选取多组参数实例,分析同一岩性类别下,不同土样土壤水力参数差异及其对入渗补给模拟结果的影响.结果表明:形状系数 (n)、土壤进气值的倒数 (α) 与极限蒸发深度呈负相关关系,是影响极限蒸发深度的显著因素;饱和含水率 (θs)、n、饱和渗透系数 (Ks) 对于入渗补给影响的敏感性较高,三者在实际土样粒径分布、干密度和孔隙性的共同影响下发生同向变动.合理的参数敏感性分析及数值模拟调参过程需结合土壤参数的物理意义来开展. 相似文献
20.
乙烷和甲烷是深部天然气藏中的重要组成部分, 精确测量其扩散系数对于正确理解深部页岩储层中的烃类气体的分布运移及其分异具有重要意义.目前乙烷在水中的扩散系数数据局限在低压和283~333K温度范围内, 缺少天然气储层高温高压条件下的数据.本研究利用显微激光拉曼光谱, 在高压透明毛细管中原位观测了20MPa下273~393K温度范围内乙烷在纯水中的扩散, 测定了各温度下的扩散系数, 并用Speedy-Angell指数方程拟合出乙烷扩散系数D(乙烷)(m2/s)与温度T(K)之间的关系式: D(C2H6)=D0[(T/Ts)-1]γ, 式中: D0=13.8055×10-9m2/s, Ts=237.4K, γ=1.7397.相同温度压力条件下, 测得的乙烷的扩散系数小于甲烷的扩散系数.据此计算了2种气体通过低渗透盖层的扩散量的差异, 发现甲烷和乙烷溶解扩散的分异程度随盖层厚度、扩散时间而显著变化. 相似文献