GEOFLUIDSgydF4y2Ba GeofluidsgydF4y2Ba 1468 - 8123gydF4y2Ba 1468 - 8115gydF4y2Ba HindawigydF4y2Ba 10.1155 / 2021/5473719gydF4y2Ba 5473719gydF4y2Ba 研究文章gydF4y2Ba 迁移的工业废水基于Thermal-Hydraulic-Mechanical破碎岩体耦合模型gydF4y2Ba https://orcid.org/0000 - 0002 - 0482 - 6235gydF4y2Ba 高gydF4y2Ba 延安gydF4y2Ba 1gydF4y2Ba 2gydF4y2Ba 3gydF4y2Ba 4gydF4y2Ba 郭gydF4y2Ba 彭gydF4y2Ba 1gydF4y2Ba 2gydF4y2Ba https://orcid.org/0000 - 0002 - 7182 - 2787gydF4y2Ba 张gydF4y2Ba ZetiangydF4y2Ba 3gydF4y2Ba https://orcid.org/0000 - 0001 - 6767 - 0234gydF4y2Ba 李gydF4y2Ba 结果表明gydF4y2Ba 4gydF4y2Ba 高gydF4y2Ba 冯gydF4y2Ba 1gydF4y2Ba 2gydF4y2Ba 王gydF4y2Ba 余gydF4y2Ba 1gydF4y2Ba 地质力学和地下工程国家重点实验室gydF4y2Ba 中国矿业大学和技术gydF4y2Ba 徐州gydF4y2Ba 江苏221116年gydF4y2Ba 中国gydF4y2Ba cumt.edu.cngydF4y2Ba 2gydF4y2Ba 力学和土木工程学院gydF4y2Ba 中国矿业大学和技术gydF4y2Ba 徐州gydF4y2Ba 江苏221116年gydF4y2Ba 中国gydF4y2Ba cumt.edu.cngydF4y2Ba 3gydF4y2Ba 地球科学与工程重点实验室的深度gydF4y2Ba 四川大学gydF4y2Ba 成都gydF4y2Ba 四川610025gydF4y2Ba 中国gydF4y2Ba scu.edu.cngydF4y2Ba 4gydF4y2Ba 广东省重点实验室的地球科学和地热能开发利用gydF4y2Ba 深圳大学gydF4y2Ba 深圳518060年gydF4y2Ba 中国gydF4y2Ba szu.edu.cngydF4y2Ba 2021年gydF4y2Ba 16gydF4y2Ba 10gydF4y2Ba 2021年gydF4y2Ba 2021年gydF4y2Ba 8gydF4y2Ba 9gydF4y2Ba 2021年gydF4y2Ba 1gydF4y2Ba 10gydF4y2Ba 2021年gydF4y2Ba 16gydF4y2Ba 10gydF4y2Ba 2021年gydF4y2Ba 2021年gydF4y2Ba 版权©2021年延安高et al。gydF4y2Ba 这是一个开放的文章在知识共享归属许可下发布的,它允许无限制的使用,分布和繁殖在任何媒介,提供最初的工作是正确的引用。gydF4y2Ba

工业废水可能长期影响环境和人类生活不断地下,造成严重的污染。好理解迁移的废水是工业废水处理的基本任务以及工业设计。研究工业废水的迁移机制在岩层中,控制方程如力学、渗流,热,和传质进行了综述,引用,并提出。热(T)液压(H)机械(M)耦合模型,建立了matrix-fault和matrix-fracture-fault多媒体。断层和裂缝的影响压力分布和污染物迁移进行了分析。断层的影响长度、宽度、倾角、渗透率和温度的废水污染物迁移参数化研究。可以得到以下结果。(1)裂缝特征影响的浓度分布,断层在浓度分布和污染物迁移。(2)污染物的迁移可以几何分为3区沿着错误的方向:饱和区,快速扩散区,浓度降低。(3)有一个峰值浓度的模型的底部。 The position of the peak is the projection of the endpoint of the fault. (4) The fault length has the most significant effect on contaminant accumulation. The temperature of the wastewater has the minimum effect on the contaminant accumulation. (5) The accumulation of concentrations can be divided into 2 stages, the slow growth stage (before 20 years) and the rapid growth stage (after 20 years). The main channel of contaminant migration in the slow growth stage is a fault. During the rapid growth stage, the contaminants penetrate through the rock matrix as well as the fault.

 深圳大学gydF4y2Ba 2020 - 3gydF4y2Ba 中华人民共和国教育部的gydF4y2Ba DESE202004gydF4y2Ba DESE202106gydF4y2Ba 四川大学gydF4y2Ba 中国国家自然科学基金gydF4y2Ba 51827901gydF4y2Ba 52078477gydF4y2Ba 1。介绍gydF4y2Ba

中国巨大的社会和经济发展带来的快速增长的行业。然而,工业的发展不可避免地随着污水排放导致环境问题不可避免。包含污染物会导致长期的废水污染和对社会生活与生态风险,因为它迁移不断地下,如土壤和岩石。受污染的土壤或岩石等因此变得有毒、有害、对环境和人类生活造成不利影响gydF4y2Ba 1gydF4y2Ba- - - - - -gydF4y2Ba 3gydF4y2Ba]。因此,研究工业废水中污染物的迁移过程岩层学者一直是一个热门话题在过去几十年。gydF4y2Ba

同时存在的不同尺度的物理结构在岩层,如矩阵,裂缝网络,和错误,导致强烈非齐次的行为(gydF4y2Ba 4gydF4y2Ba- - - - - -gydF4y2Ba 7gydF4y2Ba]。污染物在迁移过程中渗透甚至微尺度到宏观尺度。从达西线性渗流流动状态也改变岩石和天然微裂缝网络的复杂非线性的缺点。与此同时,工业废水中污染物的运输岩层是多个物理领域的结合的结果,包括温度场、液压领域,和应力场。流体压力的变化在液压领域由于流体流动行为和影响岩层的应力分布,即:,它对应力场的影响。同时,热特性是受流的速度变化的影响。岩层的孔隙度和渗透率的变化由于应力场的变化也会影响液压领域。温度的变化温度场影响水动力粘性系数和流体密度,从而影响液压领域。热应力的影响应力场的应力分布。gydF4y2Ba

如上所述,工业废水的迁移的过程是一个复杂的多重物理量耦合问题,如热工(TH),热化学(TC),胎面和(TM)和THM耦合gydF4y2Ba 8gydF4y2Ba- - - - - -gydF4y2Ba 14gydF4y2Ba]。基于节能模型,孔隙流体质量守恒方程,和弹塑性本构模型,更全面的多组分THM耦合模型建立了非饱和多孔介质由Seetharam et al。gydF4y2Ba 15gydF4y2Ba,gydF4y2Ba 16gydF4y2Ba]。然而,强的岩层的贡献没有考虑流体运动的不均匀性。等人建立了一个THMC耦合模型粘土壁垒高放射性废物存储[gydF4y2Ba 17gydF4y2Ba]。然而,孔隙流体流动和溶质扩散特征由合并后的温度梯度和浓度梯度不包括在内。研究流体渗流的数学模型,许多研究已经开展。提出了三种典型岩石渗流的基本模型,即。,the equivalent continuous medium seepage model, the fracture network discontinuous medium seepage model, and the dual medium seepage model [ 18gydF4y2Ba- - - - - -gydF4y2Ba 20.gydF4y2Ba]。雷森指出pore-fissure和interfracture渗流问题可以解决由达西定律和一个各向异性渗透率张量,提出了一种方法来描述裂隙岩体的各向异性渗流利用连续介质理论(gydF4y2Ba 21gydF4y2Ba]。这样等效模型不适用,当裂隙岩体不能被视为连续介质,和不连续介质方法应采用(gydF4y2Ba 22gydF4y2Ba]。谢长廷等人水动力学引入裂缝孔隙介质通过计算方程和连续性方程作为控制方程,提供了一种新的方法来解决这个困难的问题(gydF4y2Ba 23gydF4y2Ba]。gydF4y2Ba

大多数现有的研究不同时考虑岩体的不均匀性和多重物理量的耦合效应在污染物的迁移。因此,THM耦合模型的考虑多媒体,即。、matrix-fracture-fault提出。断层的影响和自然裂缝压力分布和污染物浓度在迁移过程中首先进行了分析。其次,断层的影响长度、宽度、倾角、渗透率,流入温度的废水污染物迁移进行了研究。本文的工业废水中污染物不指定为这是一个理论和参数研究。污染物可以被称为常见的重金属离子,如铅(铅)或汞(水银)。gydF4y2Ba

 2。THM耦合模型gydF4y2Ba

岩层中污染物的渗透的结果结合物理、化学、生物、和其他因素。因此,研究岩层中污染物迁移需要建立控制方程的耦合下multifield等化学领域,温度场、应力场和液压领域。同时,污染物的迁移包括对流、吸附、解吸、水动力弥散,化学作用,微生物分解,研究造成很大的困难。gydF4y2Ba

本文研究的重点是渗流机制下的岩层中污染物THM耦合效应。工业废水被认为是单向流动的,和吸附的影响,解吸,化学交互作用,微生物分解污染物的迁移期间不考虑。gydF4y2Ba

此外,开展数值模拟,下面的极其复杂的流程的简化和假设提出了THM耦合仍然需要建立理论模型。(1)岩石饱和,各向同性,均匀,小变形多孔弹性介质。(2)不考虑相变,即。,the transformation between gas and liquid phases due to temperature is not considered. (3) Only the heat conduction mode is considered in the rock mass, and thermal radiation is not considered. (4) The thermal strain in the rock mass due to temperature is isotropic.

2.1。等效连续介质模型理论gydF4y2Ba 2.1.1。液压领域gydF4y2Ba

废水的密度定义为温度和压力的函数(gydF4y2Ba 24gydF4y2Ba),可以表示为gydF4y2Ba (1)gydF4y2Ba ρgydF4y2Ba =gydF4y2Ba ρgydF4y2Ba 0gydF4y2Ba 1gydF4y2Ba +gydF4y2Ba ηgydF4y2Ba wgydF4y2Ba pgydF4y2Ba −gydF4y2Ba pgydF4y2Ba 0gydF4y2Ba −gydF4y2Ba βgydF4y2Ba wgydF4y2Ba TgydF4y2Ba −gydF4y2Ba TgydF4y2Ba 0gydF4y2Ba =gydF4y2Ba ρgydF4y2Ba 0gydF4y2Ba 1gydF4y2Ba +gydF4y2Ba ηgydF4y2Ba wgydF4y2Ba ΔgydF4y2Ba pgydF4y2Ba −gydF4y2Ba βgydF4y2Ba wgydF4y2Ba ΔgydF4y2Ba TgydF4y2Ba ,gydF4y2Ba 在哪里gydF4y2Ba ηgydF4y2Ba wgydF4y2Ba 是流体的压缩系数,gydF4y2Ba βgydF4y2Ba wgydF4y2Ba 的热膨胀系数是液体,gydF4y2Ba ρgydF4y2Ba 0gydF4y2Ba 液体参考密度,gydF4y2Ba pgydF4y2Ba 和gydF4y2Ba pgydF4y2Ba 0gydF4y2Ba 孔隙水压力和初始孔隙水压力,分别和gydF4y2Ba TgydF4y2Ba 和gydF4y2Ba TgydF4y2Ba 0gydF4y2Ba 分别是流体温度和初始温度。gydF4y2Ba

一般地面温度梯度的平均值是3°C / 100;因此,污水流的温度在岩层低于40°C。gydF4y2Ba

水动力粘度随着温度增加而减小,水动力粘度与温度的关系符合幂函数。废水的动态粘度与温度的表达如下(gydF4y2Ba 25gydF4y2Ba,gydF4y2Ba 26gydF4y2Ba]:gydF4y2Ba (2)gydF4y2Ba μgydF4y2Ba =gydF4y2Ba 1.787gydF4y2Ba ×gydF4y2Ba 10gydF4y2Ba −gydF4y2Ba 3gydF4y2Ba egydF4y2Ba −gydF4y2Ba 0.333gydF4y2Ba +gydF4y2Ba 1.962gydF4y2Ba ×gydF4y2Ba 10gydF4y2Ba −gydF4y2Ba 4gydF4y2Ba ×gydF4y2Ba TgydF4y2Ba ×gydF4y2Ba TgydF4y2Ba 。gydF4y2Ba

岩层中的工业废水流可以被假定为多孔介质中流体的流动。根据质量守恒定律和达西定律,液压领域方程可以表示为gydF4y2Ba (3)gydF4y2Ba ϕgydF4y2Ba ρgydF4y2Ba cgydF4y2Ba pgydF4y2Ba ∂gydF4y2Ba pgydF4y2Ba ∂gydF4y2Ba tgydF4y2Ba +gydF4y2Ba ∇gydF4y2Ba ⋅gydF4y2Ba −gydF4y2Ba ρgydF4y2Ba kgydF4y2Ba μgydF4y2Ba ∇gydF4y2Ba pgydF4y2Ba +gydF4y2Ba ρgydF4y2Ba ggydF4y2Ba ∇gydF4y2Ba HgydF4y2Ba =gydF4y2Ba 问gydF4y2Ba 米gydF4y2Ba −gydF4y2Ba ϕgydF4y2Ba ρgydF4y2Ba cgydF4y2Ba TgydF4y2Ba ∂gydF4y2Ba TgydF4y2Ba ∂gydF4y2Ba tgydF4y2Ba −gydF4y2Ba ρgydF4y2Ba ∂gydF4y2Ba ϕgydF4y2Ba ∂gydF4y2Ba tgydF4y2Ba ,gydF4y2Ba 在哪里gydF4y2Ba cgydF4y2Ba pgydF4y2Ba =gydF4y2Ba 1gydF4y2Ba /gydF4y2Ba ρgydF4y2Ba ∂gydF4y2Ba ρgydF4y2Ba /gydF4y2Ba ∂gydF4y2Ba pgydF4y2Ba ,gydF4y2Ba cgydF4y2Ba TgydF4y2Ba =gydF4y2Ba 1gydF4y2Ba /gydF4y2Ba ρgydF4y2Ba ∂gydF4y2Ba ρgydF4y2Ba /gydF4y2Ba ∂gydF4y2Ba TgydF4y2Ba ,gydF4y2Ba ϕgydF4y2Ba 孔隙度,gydF4y2Ba ρgydF4y2Ba 废水的密度,gydF4y2Ba kgydF4y2Ba 矩阵的渗透性,gydF4y2Ba ggydF4y2Ba 重力加速度,gydF4y2Ba ∇gydF4y2Ba HgydF4y2Ba =gydF4y2Ba 0gydF4y2Ba ,gydF4y2Ba 0gydF4y2Ba ,gydF4y2Ba 1gydF4y2Ba TgydF4y2Ba ,gydF4y2Ba HgydF4y2Ba 代表了水头,gydF4y2Ba 问gydF4y2Ba 米gydF4y2Ba 源项。gydF4y2Ba

 2.1.2。温度场gydF4y2Ba

岩体中的流体流动在热平衡假设。不考虑岩体之间的热交换和废水,流体流动过程的温度场方程在饱和多孔介质可以获得gydF4y2Ba 27gydF4y2Ba]gydF4y2Ba (4)gydF4y2Ba ρgydF4y2Ba CgydF4y2Ba effgydF4y2Ba ∂gydF4y2Ba TgydF4y2Ba ∂gydF4y2Ba tgydF4y2Ba +gydF4y2Ba 1gydF4y2Ba −gydF4y2Ba ϕgydF4y2Ba TgydF4y2Ba KgydF4y2Ba αgydF4y2Ba TgydF4y2Ba ∂gydF4y2Ba εgydF4y2Ba vgydF4y2Ba ∂gydF4y2Ba tgydF4y2Ba =gydF4y2Ba λgydF4y2Ba effgydF4y2Ba ∇gydF4y2Ba 2gydF4y2Ba TgydF4y2Ba +gydF4y2Ba ρgydF4y2Ba cgydF4y2Ba wgydF4y2Ba kgydF4y2Ba μgydF4y2Ba ∇gydF4y2Ba pgydF4y2Ba ∇gydF4y2Ba TgydF4y2Ba ,gydF4y2Ba 在哪里gydF4y2Ba ρgydF4y2Ba CgydF4y2Ba effgydF4y2Ba =gydF4y2Ba ϕgydF4y2Ba ρgydF4y2Ba cgydF4y2Ba wgydF4y2Ba +gydF4y2Ba 1gydF4y2Ba −gydF4y2Ba ϕgydF4y2Ba ρgydF4y2Ba 年代gydF4y2Ba cgydF4y2Ba 年代gydF4y2Ba 是有效的热容,gydF4y2Ba λgydF4y2Ba effgydF4y2Ba =gydF4y2Ba ϕgydF4y2Ba λgydF4y2Ba wgydF4y2Ba +gydF4y2Ba 1gydF4y2Ba −gydF4y2Ba ϕgydF4y2Ba λgydF4y2Ba 年代gydF4y2Ba 有效导热系数,gydF4y2Ba αgydF4y2Ba TgydF4y2Ba 是岩石的热膨胀系数矩阵,gydF4y2Ba εgydF4y2Ba vgydF4y2Ba 是岩石的体积应变矩阵,gydF4y2Ba cgydF4y2Ba wgydF4y2Ba 比热容的废水,gydF4y2Ba λgydF4y2Ba wgydF4y2Ba 是废水的导热系数,gydF4y2Ba ρgydF4y2Ba 年代gydF4y2Ba 岩体的密度,gydF4y2Ba cgydF4y2Ba 年代gydF4y2Ba 是岩石的比热容矩阵,然后呢gydF4y2Ba λgydF4y2Ba 年代gydF4y2Ba 是岩石的导热系数矩阵。gydF4y2Ba

 2.1.3。应力场gydF4y2Ba

岩体的本构关系是理想的各向同性,均匀、线性与小变形介质。岩体的变形将造成的3个因素:岩体的应力,流体的压力岩体,岩体的温度变化。因此,岩体的应力-应变关系和考虑孔隙压力和温度是获得基于线性原则:gydF4y2Ba (5)gydF4y2Ba σgydF4y2Ba 我gydF4y2Ba jgydF4y2Ba =gydF4y2Ba σgydF4y2Ba 我gydF4y2Ba jgydF4y2Ba ′gydF4y2Ba +gydF4y2Ba αgydF4y2Ba pgydF4y2Ba δgydF4y2Ba 我gydF4y2Ba jgydF4y2Ba =gydF4y2Ba 2gydF4y2Ba GgydF4y2Ba εgydF4y2Ba 我gydF4y2Ba jgydF4y2Ba +gydF4y2Ba λgydF4y2Ba δgydF4y2Ba 我gydF4y2Ba jgydF4y2Ba εgydF4y2Ba vgydF4y2Ba −gydF4y2Ba KgydF4y2Ba αgydF4y2Ba TgydF4y2Ba δgydF4y2Ba 我gydF4y2Ba jgydF4y2Ba TgydF4y2Ba −gydF4y2Ba TgydF4y2Ba 0gydF4y2Ba +gydF4y2Ba αgydF4y2Ba pgydF4y2Ba δgydF4y2Ba 我gydF4y2Ba jgydF4y2Ba ,gydF4y2Ba 在哪里gydF4y2Ba σgydF4y2Ba 我gydF4y2Ba jgydF4y2Ba 应力张量,gydF4y2Ba σgydF4y2Ba 我gydF4y2Ba jgydF4y2Ba ′gydF4y2Ba 有效应力张量,gydF4y2Ba GgydF4y2Ba 剪切模量,gydF4y2Ba GgydF4y2Ba =gydF4y2Ba EgydF4y2Ba /gydF4y2Ba 2gydF4y2Ba 1gydF4y2Ba +gydF4y2Ba νgydF4y2Ba ,gydF4y2Ba εgydF4y2Ba 我gydF4y2Ba jgydF4y2Ba 应变张量,gydF4y2Ba λgydF4y2Ba 是瘸腿的常数,gydF4y2Ba λgydF4y2Ba =gydF4y2Ba 2gydF4y2Ba GgydF4y2Ba νgydF4y2Ba /gydF4y2Ba 1gydF4y2Ba −gydF4y2Ba 2gydF4y2Ba νgydF4y2Ba ,gydF4y2Ba αgydF4y2Ba 毕奥系数:gydF4y2Ba (6)gydF4y2Ba αgydF4y2Ba =gydF4y2Ba 1gydF4y2Ba −gydF4y2Ba KgydF4y2Ba KgydF4y2Ba 年代gydF4y2Ba ,gydF4y2Ba 在哪里gydF4y2Ba KgydF4y2Ba 是岩石的体积变形模量矩阵,gydF4y2Ba KgydF4y2Ba =gydF4y2Ba EgydF4y2Ba /gydF4y2Ba 3gydF4y2Ba 1gydF4y2Ba −gydF4y2Ba 2gydF4y2Ba νgydF4y2Ba ,gydF4y2Ba KgydF4y2Ba 年代gydF4y2Ba 岩石颗粒的体积变形模量。gydF4y2Ba

基于柯西应变理论和静力平衡,平衡方程的几何关系,可以得到:gydF4y2Ba (7)gydF4y2Ba εgydF4y2Ba 我gydF4y2Ba jgydF4y2Ba =gydF4y2Ba 1gydF4y2Ba 2gydF4y2Ba ugydF4y2Ba 我gydF4y2Ba ,gydF4y2Ba jgydF4y2Ba +gydF4y2Ba ugydF4y2Ba jgydF4y2Ba ,gydF4y2Ba 我gydF4y2Ba ,gydF4y2Ba (8)gydF4y2Ba εgydF4y2Ba vgydF4y2Ba =gydF4y2Ba εgydF4y2Ba kgydF4y2Ba kgydF4y2Ba =gydF4y2Ba εgydF4y2Ba 11gydF4y2Ba +gydF4y2Ba εgydF4y2Ba 22gydF4y2Ba +gydF4y2Ba εgydF4y2Ba 33gydF4y2Ba ,gydF4y2Ba (9)gydF4y2Ba σgydF4y2Ba 我gydF4y2Ba jgydF4y2Ba ,gydF4y2Ba jgydF4y2Ba +gydF4y2Ba fgydF4y2Ba 我gydF4y2Ba =gydF4y2Ba 0gydF4y2Ba 。gydF4y2Ba

结合方程(gydF4y2Ba 5gydF4y2Ba),我们可以得到的方程下岩体的应力场应力、孔隙压力和温度:gydF4y2Ba (10)gydF4y2Ba GgydF4y2Ba ugydF4y2Ba 我gydF4y2Ba ,gydF4y2Ba 我gydF4y2Ba jgydF4y2Ba +gydF4y2Ba ugydF4y2Ba jgydF4y2Ba ,gydF4y2Ba 我gydF4y2Ba jgydF4y2Ba +gydF4y2Ba λgydF4y2Ba ugydF4y2Ba jgydF4y2Ba ,gydF4y2Ba jgydF4y2Ba 我gydF4y2Ba −gydF4y2Ba αgydF4y2Ba TgydF4y2Ba TgydF4y2Ba ,gydF4y2Ba 我gydF4y2Ba +gydF4y2Ba αgydF4y2Ba pgydF4y2Ba ,gydF4y2Ba 我gydF4y2Ba +gydF4y2Ba FgydF4y2Ba 我gydF4y2Ba =gydF4y2Ba 0gydF4y2Ba ,gydF4y2Ba 在哪里gydF4y2Ba FgydF4y2Ba 我gydF4y2Ba 岩体的体积力和吗gydF4y2Ba ugydF4y2Ba 我gydF4y2Ba 是岩体的位移分量。gydF4y2Ba

 2.1.4。孔隙度模型gydF4y2Ba

孔隙度是孔隙压力的关键因素;因此,有必要选择一个岩体的孔隙度模型。根据孔隙度的定义,推导了岩体的动态孔隙度方程定义如下(gydF4y2Ba 28gydF4y2Ba]:gydF4y2Ba (11)gydF4y2Ba ϕgydF4y2Ba =gydF4y2Ba αgydF4y2Ba −gydF4y2Ba αgydF4y2Ba −gydF4y2Ba ϕgydF4y2Ba 0gydF4y2Ba 经验值gydF4y2Ba −gydF4y2Ba 1gydF4y2Ba KgydF4y2Ba 年代gydF4y2Ba pgydF4y2Ba −gydF4y2Ba pgydF4y2Ba 0gydF4y2Ba +gydF4y2Ba αgydF4y2Ba TgydF4y2Ba TgydF4y2Ba −gydF4y2Ba TgydF4y2Ba 0gydF4y2Ba −gydF4y2Ba εgydF4y2Ba vgydF4y2Ba −gydF4y2Ba εgydF4y2Ba vgydF4y2Ba 0gydF4y2Ba 。gydF4y2Ba

2.1.5节讨论。渗透率模型gydF4y2Ba

Kozeny-Carman (KC)方程(gydF4y2Ba 29日gydF4y2Ba)是一种半经验的公式广泛应用于多孔介质的渗透率演化的研究。基于孔隙度模型、渗透率模型可以表示为gydF4y2Ba (12)gydF4y2Ba kgydF4y2Ba =gydF4y2Ba kgydF4y2Ba 0gydF4y2Ba 1gydF4y2Ba ϕgydF4y2Ba 0gydF4y2Ba −gydF4y2Ba 1gydF4y2Ba ϕgydF4y2Ba 0gydF4y2Ba −gydF4y2Ba 1gydF4y2Ba 1gydF4y2Ba −gydF4y2Ba 1gydF4y2Ba KgydF4y2Ba 年代gydF4y2Ba pgydF4y2Ba −gydF4y2Ba pgydF4y2Ba 0gydF4y2Ba +gydF4y2Ba αgydF4y2Ba TgydF4y2Ba TgydF4y2Ba −gydF4y2Ba TgydF4y2Ba 0gydF4y2Ba −gydF4y2Ba εgydF4y2Ba vgydF4y2Ba −gydF4y2Ba εgydF4y2Ba vgydF4y2Ba 0gydF4y2Ba 3gydF4y2Ba 。gydF4y2Ba

2.2。THM耦合关系gydF4y2Ba

基于定义的THM耦合模型如图以上提供的场方程gydF4y2Ba 1gydF4y2Ba。(1)液压领域影响机械领域通过改变流体的压力。(2)机械领域的渗透率和孔隙度的影响通过体积应变的变化,从而影响液压领域。(3)液压领域影响温度场的热对流项通过流体速度矢量的变化。(4)温度场影响流体密度通过改变温度,然后影响液压领域。(5)机械领域改变了热变形能量影响温度场。(6)温度场的温度变化产生的热应力,然后影响机械领域。gydF4y2Ba

THM耦合关系图。gydF4y2Ba

 2.3。Matrix-Fracture-Fault (MFF)多媒体模型gydF4y2Ba

matrix-fracture-fault (MFF)多媒体模型不同于等效连续介质模型。它分开考虑多孔介质断裂裂缝介质和离散化的每一个裂缝的影响在液压和热能领域。裂缝网络的持续改进,MFF多媒体模型趋势走向完美的建模(gydF4y2Ba 30.gydF4y2Ba]。gydF4y2Ba

2.3.1。液压领域gydF4y2Ba

控制裂缝的渗流方程可以表示如下(gydF4y2Ba 28gydF4y2Ba]:gydF4y2Ba (13)gydF4y2Ba dgydF4y2Ba fgydF4y2Ba 年代gydF4y2Ba fgydF4y2Ba ∂gydF4y2Ba pgydF4y2Ba ∂gydF4y2Ba tgydF4y2Ba +gydF4y2Ba ∇gydF4y2Ba TgydF4y2Ba ⋅gydF4y2Ba −gydF4y2Ba dgydF4y2Ba fgydF4y2Ba kgydF4y2Ba fgydF4y2Ba μgydF4y2Ba ∇gydF4y2Ba TgydF4y2Ba pgydF4y2Ba =gydF4y2Ba −gydF4y2Ba dgydF4y2Ba fgydF4y2Ba ∂gydF4y2Ba εgydF4y2Ba vgydF4y2Ba ∂gydF4y2Ba tgydF4y2Ba +gydF4y2Ba 问gydF4y2Ba fgydF4y2Ba ,gydF4y2Ba 在哪里gydF4y2Ba dgydF4y2Ba fgydF4y2Ba 裂缝宽度,gydF4y2Ba 年代gydF4y2Ba fgydF4y2Ba 裂缝的储水系数,gydF4y2Ba ∇gydF4y2Ba TgydF4y2Ba 是切向裂缝边界方向的梯度算子,即边界上的切向分化的功能gydF4y2Ba xgydF4y2Ba 方向和gydF4y2Ba ygydF4y2Ba 方向,gydF4y2Ba kgydF4y2Ba fgydF4y2Ba 裂缝渗透率。gydF4y2Ba

 2.3.2。温度场gydF4y2Ba

与等效连续介质情况下,当离散裂缝模型中,离散裂缝成为主要的流体流动通道矩阵。因此,对流的影响作用在矩阵模型的温度场可以被忽视。热传导效应被认为是矩阵中,因此,裂缝中的温度场的控制方程可以表示如下:gydF4y2Ba (14)gydF4y2Ba dgydF4y2Ba fgydF4y2Ba ρgydF4y2Ba CgydF4y2Ba effgydF4y2Ba ∂gydF4y2Ba TgydF4y2Ba ∂gydF4y2Ba tgydF4y2Ba =gydF4y2Ba ∇gydF4y2Ba TgydF4y2Ba ⋅gydF4y2Ba dgydF4y2Ba fgydF4y2Ba λgydF4y2Ba effgydF4y2Ba ∇gydF4y2Ba TgydF4y2Ba TgydF4y2Ba +gydF4y2Ba dgydF4y2Ba fgydF4y2Ba ρgydF4y2Ba cgydF4y2Ba wgydF4y2Ba kgydF4y2Ba μgydF4y2Ba ∇gydF4y2Ba pgydF4y2Ba ∇gydF4y2Ba TgydF4y2Ba TgydF4y2Ba +gydF4y2Ba 问gydF4y2Ba wgydF4y2Ba 。gydF4y2Ba

2.3.3。应力场gydF4y2Ba

一样的压力控制方程的等效连续介质,断裂的应力场控制方程可以表示为一部分gydF4y2Ba (15)gydF4y2Ba GgydF4y2Ba ugydF4y2Ba 我gydF4y2Ba ,gydF4y2Ba 我gydF4y2Ba jgydF4y2Ba +gydF4y2Ba ugydF4y2Ba jgydF4y2Ba ,gydF4y2Ba 我gydF4y2Ba jgydF4y2Ba +gydF4y2Ba λgydF4y2Ba ugydF4y2Ba jgydF4y2Ba ,gydF4y2Ba jgydF4y2Ba 我gydF4y2Ba −gydF4y2Ba αgydF4y2Ba TgydF4y2Ba TgydF4y2Ba ,gydF4y2Ba 我gydF4y2Ba +gydF4y2Ba αgydF4y2Ba pgydF4y2Ba ,gydF4y2Ba 我gydF4y2Ba +gydF4y2Ba FgydF4y2Ba 我gydF4y2Ba =gydF4y2Ba 0gydF4y2Ba 。gydF4y2Ba

由于本文研究的是一个二维扩散模型,离散骨折被视为线的几何模型。因此,骨折的变形是被忽视的。骨折的渗透率与应力之间的关系可以建立基于水的井下测试压力测量(gydF4y2Ba 31日gydF4y2Ba]:gydF4y2Ba (16)gydF4y2Ba kgydF4y2Ba fgydF4y2Ba =gydF4y2Ba kgydF4y2Ba fgydF4y2Ba 0gydF4y2Ba 经验值gydF4y2Ba −gydF4y2Ba αgydF4y2Ba fgydF4y2Ba σgydF4y2Ba ngydF4y2Ba ,gydF4y2Ba 在哪里gydF4y2Ba kgydF4y2Ba fgydF4y2Ba 0gydF4y2Ba 裂缝的渗透率是什么时候gydF4y2Ba σgydF4y2Ba ngydF4y2Ba =gydF4y2Ba 0gydF4y2Ba ,gydF4y2Ba σgydF4y2Ba ngydF4y2Ba 相对应的正常力的断裂表面gydF4y2Ba αgydF4y2Ba fgydF4y2Ba 是影响系数,它决定了断裂状态。gydF4y2Ba

 2.4。污染物迁移模型gydF4y2Ba

多孔介质中污染物的迁移机制由3部分组成:对流,即。,污染物是由下游水流迁移;扩散,污染物从高浓度区向低浓度区扩散在浓度梯度的影响;和分散,这是由于多孔介质的存在,和多孔介质导致污染物的迁移速度不同的平均流速的大小和方向。尽管扩散和机械分散的机制是不同的,他们仍然难以区分。他们通常命名为水动力弥散。污染物在多孔介质中迁移的机制可以表示如下:gydF4y2Ba (17)gydF4y2Ba ∂gydF4y2Ba cgydF4y2Ba ∂gydF4y2Ba tgydF4y2Ba =gydF4y2Ba ∇gydF4y2Ba ⋅gydF4y2Ba DgydF4y2Ba ∇gydF4y2Ba cgydF4y2Ba −gydF4y2Ba ∇gydF4y2Ba ⋅gydF4y2Ba ugydF4y2Ba cgydF4y2Ba +gydF4y2Ba 我gydF4y2Ba ,gydF4y2Ba 在哪里gydF4y2Ba cgydF4y2Ba 是溶质浓度,gydF4y2Ba DgydF4y2Ba 水动力弥散系数,gydF4y2Ba ugydF4y2Ba 是对流速度。gydF4y2Ba

 3所示。迁移在多媒体THM耦合效应gydF4y2Ba

调查THM耦合下的污染物的迁移效果的影响以及故障和骨折,2模型被使用,例如,gydF4y2Ba

matrix-fault (MF)模型gydF4y2Ba

MFF模型gydF4y2Ba

MF模型考虑了岩石矩阵和故障,和MFF模型考虑了岩石矩阵,断层和骨折。gydF4y2Ba

模型的上边界的边界流入废水。岩层的初始温度为318 K。流入污染物的温度上限是293 K。污染物的浓度是1摩尔/ mgydF4y2Ba3gydF4y2Ba。水头压力的上限是9000。参数表中列出的数值研究中使用gydF4y2Ba 1gydF4y2Ba。gydF4y2Ba

参数和值。gydF4y2Ba

参数gydF4y2Ba 价值gydF4y2Ba
岩石的杨氏模量矩阵,gydF4y2Ba EgydF4y2Ba (Pa)gydF4y2Ba 2gydF4y2Ba ×gydF4y2Ba 10gydF4y2Ba 10gydF4y2Ba
杨氏模量的岩石颗粒,gydF4y2Ba EgydF4y2Ba 年代gydF4y2Ba (Pa)gydF4y2Ba 5gydF4y2Ba ×gydF4y2Ba 10gydF4y2Ba 10gydF4y2Ba
泊松比,gydF4y2Ba νgydF4y2Ba 0.3gydF4y2Ba
岩石的热膨胀系数矩阵,gydF4y2Ba αgydF4y2Ba TgydF4y2Ba (1 / K)gydF4y2Ba 1。2gydF4y2Ba ×gydF4y2Ba 10gydF4y2Ba −gydF4y2Ba 5gydF4y2Ba
岩石颗粒的密度,gydF4y2Ba ρgydF4y2Ba 年代gydF4y2Ba (公斤/米gydF4y2Ba3gydF4y2Ba)gydF4y2Ba 2652年gydF4y2Ba
比热容的矩阵,gydF4y2Ba cgydF4y2Ba 年代gydF4y2Ba (J /(公斤·°C))gydF4y2Ba 2190年gydF4y2Ba
岩石的导热系数矩阵,gydF4y2Ba λgydF4y2Ba 年代gydF4y2Ba (W / (m·°C))gydF4y2Ba 0.56gydF4y2Ba
废水的密度,gydF4y2Ba ρgydF4y2Ba 0gydF4y2Ba (公斤/米gydF4y2Ba3gydF4y2Ba)gydF4y2Ba 1100年gydF4y2Ba
比热容的废水,gydF4y2Ba cgydF4y2Ba wgydF4y2Ba (J /(公斤·°C))gydF4y2Ba 4210年gydF4y2Ba
热导率的废水,gydF4y2Ba λgydF4y2Ba wgydF4y2Ba (W / (m·°C))gydF4y2Ba 0.65gydF4y2Ba
初始孔隙度,gydF4y2Ba φgydF4y2Ba 0gydF4y2Ba 0.001gydF4y2Ba
矩阵的初始磁导率,gydF4y2Ba kgydF4y2Ba 0gydF4y2Ba (mgydF4y2Ba2gydF4y2Ba)gydF4y2Ba 2。5gydF4y2Ba ×gydF4y2Ba 10gydF4y2Ba −gydF4y2Ba 15gydF4y2Ba
初始渗透率的错,gydF4y2Ba kgydF4y2Ba 1gydF4y2Ba (mgydF4y2Ba2gydF4y2Ba)gydF4y2Ba 1gydF4y2Ba ×gydF4y2Ba 10gydF4y2Ba −gydF4y2Ba 13gydF4y2Ba
初始渗透率的骨折,gydF4y2Ba kgydF4y2Ba 2gydF4y2Ba (mgydF4y2Ba2gydF4y2Ba)gydF4y2Ba 2。5gydF4y2Ba ×gydF4y2Ba 10gydF4y2Ba −gydF4y2Ba 14gydF4y2Ba
断层的宽度,gydF4y2Ba dgydF4y2Ba 1gydF4y2Ba (m)gydF4y2Ba 0.002gydF4y2Ba
孔径的骨折,gydF4y2Ba dgydF4y2Ba 2gydF4y2Ba (m)gydF4y2Ba 0.0005gydF4y2Ba

研究断层周围的污染物的浓度分布和浓度的积累,2监测线(gydF4y2Ba 一个gydF4y2Ba BgydF4y2Ba 和gydF4y2Ba CgydF4y2Ba DgydF4y2Ba )选择。的线gydF4y2Ba 一个gydF4y2Ba BgydF4y2Ba 沿着错误的方向。的线gydF4y2Ba CgydF4y2Ba DgydF4y2Ba 底部的几何模型。端点的坐标gydF4y2Ba 一个gydF4y2Ba ,gydF4y2Ba BgydF4y2Ba ,gydF4y2Ba CgydF4y2Ba ,gydF4y2Ba 和gydF4y2Ba DgydF4y2Ba (23.66米、30米),(6.34,0米),(0,0)和(30 m, 0米)。gydF4y2Ba

3.1。MF模型gydF4y2Ba

如图gydF4y2Ba 2gydF4y2Ba的几何模型的大小gydF4y2Ba 30.gydF4y2Ba 米gydF4y2Ba ×gydF4y2Ba 30.gydF4y2Ba 米gydF4y2Ba 建立了。故障是60°的倾斜角度,断层的宽度是0.002米,断层的中点的坐标是(15米,15米)。gydF4y2Ba

MF的几何模型。gydF4y2Ba

如图gydF4y2Ba 3gydF4y2Ba显然不同,断层周围的压力等值线集中。模型中的压力分布行为非齐次由于错。断层和岩石的渗透性的差异矩阵可能是这一现象的原因,因为它可能会导致渗流场的不均匀性。gydF4y2Ba

曼氏金融的压力分布轮廓模型。gydF4y2Ba

图gydF4y2Ba 4gydF4y2Ba是二十年后污染物浓度的分布。污染物主要集中在模型的上半部分,沿着断层。这可能意味着污染物迁移非常缓慢,和扩散是更重要的错比矩阵。gydF4y2Ba

MF的浓度分布等值线模型。gydF4y2Ba

正如上面提到的,故障可能会影响的流动大大废水和污染物的迁移。流矢量是用来说明。如图gydF4y2Ba 5gydF4y2Ba两向量的方向和密度变化显然在断层;因此,它可以是另一个方面反映出角色的断层分布的不均匀性。gydF4y2Ba

流矢量的废水MF模型。gydF4y2Ba

虽然故障有一个戏剧性的影响流场,污染物的浓度监测沿断层方向,即。的线gydF4y2Ba 一个gydF4y2Ba BgydF4y2Ba 。重要的是要注意的长度gydF4y2Ba 一个gydF4y2Ba BgydF4y2Ba 本文基于计算的起点gydF4y2Ba 一个gydF4y2Ba 。如图gydF4y2Ba 6gydF4y2Ba,污染物达到饱和浓度接近1摩尔/ mgydF4y2Ba3gydF4y2Ba在3 m向错误方向的长度。和污染物的扩散变得显著沿着断层大约3米到15米。因此,可以找到快速扩散区。从15米到模型的底部(点gydF4y2Ba BgydF4y2Ba ),浓度迅速降低。因此,污染物的迁移可分为3区沿着错误的方向,即。,饱和区,快速扩散区,浓度降低。gydF4y2Ba

污染物浓度沿行gydF4y2Ba 一个gydF4y2Ba BgydF4y2Ba 在MF模型中。gydF4y2Ba

它也可以观察到浓度随时间增加。点的浓度gydF4y2Ba BgydF4y2Ba 到达gydF4y2Ba 0.23gydF4y2Ba ×gydF4y2Ba 10gydF4y2Ba −gydF4y2Ba 3gydF4y2Ba 摩尔/ mgydF4y2Ba3gydF4y2Ba,gydF4y2Ba 1.31gydF4y2Ba ×gydF4y2Ba 10gydF4y2Ba −gydF4y2Ba 3gydF4y2Ba 摩尔/ mgydF4y2Ba3gydF4y2Ba,gydF4y2Ba 1。3gydF4y2Ba ×gydF4y2Ba 10gydF4y2Ba −gydF4y2Ba 3gydF4y2Ba 摩尔/ mgydF4y2Ba3gydF4y2Ba0.012摩尔/ mgydF4y2Ba3gydF4y2Ba0.074摩尔/ mgydF4y2Ba3gydF4y2Ba和0.138摩尔/ mgydF4y2Ba3gydF4y2Ba的5年,10年,15年,20年,25年,30年,分别。gydF4y2Ba

图gydF4y2Ba 7gydF4y2Ba的情节是浓度沿行吗gydF4y2Ba CgydF4y2Ba DgydF4y2Ba ,即,the bottom of the model. Similarly, the concentration increases with time, and it becomes the maximum after 30 years. Meanwhile, it can be found that this is a peak in each line of the plot. And the peak occurs at the point of the projection of the endpoint of the fault. The reason may be as follows. The fault has a significant effect on the migration of the contaminants as its permeability is much higher than that of the rock matrix. The projection of the endpoint of the fault is the point which is the nearest to the fault; therefore, the peak value of the concentration occurs around such a point.

污染物浓度沿行gydF4y2Ba CgydF4y2Ba DgydF4y2Ba 在MF模型中。gydF4y2Ba

 3.2。MFF模型gydF4y2Ba

在本节中,MFF模型(图gydF4y2Ba 8gydF4y2Ba)考虑岩石的矩阵,采用故障,骨折研究骨折的效果。50骨折随机生成通过使用MATLAB代码(gydF4y2Ba 32gydF4y2Ba]。骨折的长度从0.8米到4米不等。骨折不同的方向从0°- 360°。每个断裂的孔径是0.0005米。gydF4y2Ba

MFF的几何模型。gydF4y2Ba

相比之下,图gydF4y2Ba 3gydF4y2Ba有更激烈的压力变化模型(图gydF4y2Ba 9gydF4y2Ba)。MFF模型的压力分布的不均匀性更明显比MF模型。周围的压力轮廓变化自然骨折。然而,主要的压力分布模式MFF MF模型的相似模型。它也可以发现,MFF的浓度的分布模型(图gydF4y2Ba 10gydF4y2Ba)是MF模型的相似。上述现象可能表明断层具有主导作用的浓度分布和污染物迁移。gydF4y2Ba

的压力分布轮廓MFF模型。gydF4y2Ba

的浓度分布轮廓MFF模型。gydF4y2Ba

基于图gydF4y2Ba 11gydF4y2Ba,它可以发现还有3个典型区域的污染物浓度,饱和区,快速扩散区,浓度降低。3区域的几何尺度也类似于MF模型。点的浓度gydF4y2Ba BgydF4y2Ba 到达gydF4y2Ba 1.15gydF4y2Ba ×gydF4y2Ba 10gydF4y2Ba −gydF4y2Ba 3gydF4y2Ba 摩尔/ mgydF4y2Ba3gydF4y2Ba,gydF4y2Ba 1.47gydF4y2Ba ×gydF4y2Ba 10gydF4y2Ba −gydF4y2Ba 3gydF4y2Ba 摩尔/ mgydF4y2Ba3gydF4y2Ba,gydF4y2Ba 1.79gydF4y2Ba ×gydF4y2Ba 10gydF4y2Ba −gydF4y2Ba 3gydF4y2Ba 摩尔/ mgydF4y2Ba3gydF4y2Ba0.016摩尔/ mgydF4y2Ba3gydF4y2Ba0.086摩尔/ mgydF4y2Ba3gydF4y2Ba和0.157摩尔/ mgydF4y2Ba3gydF4y2Ba的5年,10年,15年,20年,25年,30年,分别。gydF4y2Ba

污染物浓度沿行gydF4y2Ba 一个gydF4y2Ba BgydF4y2Ba 在MFF模型中。gydF4y2Ba

图gydF4y2Ba 12gydF4y2Ba浓度分布的情节沿着MFF模型(底部的线吗gydF4y2Ba CgydF4y2Ba DgydF4y2Ba );还有一个峰的每一行的阴谋。每个峰值MFF高于相应的MF模型之一。峰值出现的位置是几乎一样的MFF模型。因此,2的浓度分布模型是相似的,MFF模型的大小是高于MF模型。gydF4y2Ba

污染物浓度沿行gydF4y2Ba CgydF4y2Ba DgydF4y2Ba 在MFF模型中。gydF4y2Ba

 4所示。断层的影响参数对污染物迁移gydF4y2Ba

作为讨论的部分gydF4y2Ba 3gydF4y2Ba,污染物迁移的故障具有主导作用。错误的参数研究是在这一节中进行。长度、宽度、倾斜、渗透性和宽度选择错误的参数。此外,工业废水的初始温度的影响也进行了研究。gydF4y2Ba

4.1。断层的长度gydF4y2Ba

基于部分的结果gydF4y2Ba 3所示。1gydF4y2Ba,断层的长度可能控制污染物的分布;因此影响污染物迁移的一个重要因素。在本节中,与不同的断层长度(表3模型gydF4y2Ba 2gydF4y2Ba)是用来探索断层长度对污染物的迁移的影响。gydF4y2Ba

断层的长度在不同的模型。gydF4y2Ba

模型没有。gydF4y2Ba 模型L1gydF4y2Ba 模型L2gydF4y2Ba 模型L3gydF4y2Ba
断层长度(米)gydF4y2Ba 10gydF4y2Ba 20.gydF4y2Ba 30.gydF4y2Ba

图gydF4y2Ba 13gydF4y2Ba沿着线显示了污染物的积累gydF4y2Ba CgydF4y2Ba DgydF4y2Ba 根据不同的断层的长度。线的污染物的积累gydF4y2Ba CgydF4y2Ba DgydF4y2Ba 随着时间的推移,逐渐增加。30年后,积累的污染物gydF4y2Ba CgydF4y2Ba DgydF4y2Ba 到达gydF4y2Ba 3.68gydF4y2Ba ×gydF4y2Ba 10gydF4y2Ba −gydF4y2Ba 4gydF4y2Ba 摩尔(模型L1), 0.064摩尔(L2)模型,分别和1.58摩尔(L3)模型。断层的长度对污染物的积累具有明显的影响。更长的故障将导致更多的积累的污染物。gydF4y2Ba

沿着线污染物积累gydF4y2Ba CgydF4y2Ba DgydF4y2Ba 在不同故障下的长度。gydF4y2Ba

 4.2。断层的宽度gydF4y2Ba

断层的宽度可能没有显著影响污染物的迁移模式虽然与渗透率有关。它也可能对污染物积累有重要影响。在本节中,与不同的断层宽度(表3模型gydF4y2Ba 3gydF4y2Ba)建立研究断层宽度对污染物的累积的影响。gydF4y2Ba

断层的宽度在不同的模型。gydF4y2Ba

模型没有。gydF4y2Ba 模型W1gydF4y2Ba 模型W2gydF4y2Ba 模型W3gydF4y2Ba
断层宽度(米)gydF4y2Ba 0.001gydF4y2Ba 0.002gydF4y2Ba 0.003gydF4y2Ba

沿线的污染物的积累gydF4y2Ba CgydF4y2Ba DgydF4y2Ba 根据不同的断层宽度绘制在图gydF4y2Ba 14gydF4y2Ba。30年后,线的污染物的积累gydF4y2Ba CgydF4y2Ba DgydF4y2Ba 达到0.69摩尔(W1)模型,1.58摩尔(W2)模型,分别和1.83摩尔(模型W3)。更大的断层宽度可能导致一个更大的积累。这表明断层的宽度可以增加当地的模型的渗透率,提高污染物的扩散。gydF4y2Ba

沿着线污染物积累gydF4y2Ba CgydF4y2Ba DgydF4y2Ba 根据不同的断层宽度。gydF4y2Ba

 4.3。断层倾角gydF4y2Ba

污染物的迁移过程中,断层的倾角可以确定断层的长度;因此,污染物的迁移距离以及污染物的分布可能受这种倾角的影响。三个模型(表gydF4y2Ba 4gydF4y2Ba)与断层的不同倾斜角度考虑污染物的积累。gydF4y2Ba

断层倾角不同的模型。gydF4y2Ba

模型没有。gydF4y2Ba 模型D1gydF4y2Ba 模型D2gydF4y2Ba 模型D3gydF4y2Ba
断层倾角(°)gydF4y2Ba 30.gydF4y2Ba 45gydF4y2Ba 60gydF4y2Ba

图gydF4y2Ba 15gydF4y2Ba是污染物的积累的情节沿着行吗gydF4y2Ba CgydF4y2Ba DgydF4y2Ba 根据不同的断层倾斜角度。30年后,沿线的污染物的积累gydF4y2Ba CgydF4y2Ba DgydF4y2Ba 达到0.67摩尔(模型D1), 1.25摩尔(D2)模型,分别和1.58摩尔(D3)模型。60°倾角的故障可能引起的最大污染物积累。gydF4y2Ba

沿着线污染物积累gydF4y2Ba CgydF4y2Ba DgydF4y2Ba 在不同故障下的倾斜角度。gydF4y2Ba

 4.4。断层渗透率gydF4y2Ba

渗透系数是一个指数的渗透性岩层和关键影响污染物的迁移。在本节中,初始渗透率gydF4y2Ba kgydF4y2Ba 故障设置的gydF4y2Ba 1gydF4y2Ba ×gydF4y2Ba 10gydF4y2Ba −gydF4y2Ba 13gydF4y2Ba 米gydF4y2Ba2gydF4y2Ba。三个模型(表gydF4y2Ba 5gydF4y2Ba)与断层的不同渗透率是用来调查断层渗透率的影响污染物的积累。gydF4y2Ba

渗透率在不同的模型。gydF4y2Ba

模型没有。gydF4y2Ba 模型P1gydF4y2Ba 模型P2gydF4y2Ba 模型P3gydF4y2Ba
断层渗透率(mgydF4y2Ba2gydF4y2Ba)gydF4y2Ba 0.5gydF4y2Ba kgydF4y2Ba 1。0gydF4y2Ba kgydF4y2Ba 1。5gydF4y2Ba kgydF4y2Ba

图gydF4y2Ba 16gydF4y2Ba沿着线显示了污染物的积累gydF4y2Ba CgydF4y2Ba DgydF4y2Ba 随着时间的推移,不同渗透率下。30年后,污染物的积累gydF4y2Ba CgydF4y2Ba DgydF4y2Ba 达到0.67摩尔(P1)模型,1.58摩尔(P2)模型,分别和1.84摩尔(P3)模型。更大的渗透导致一个更大的污染物的积累。gydF4y2Ba

沿着线污染物积累gydF4y2Ba CgydF4y2Ba DgydF4y2Ba 在不同渗透率。gydF4y2Ba

 4.5。工业废水流入温度gydF4y2Ba

废水流入岩层时,流入温度会影响温度场的岩层,和intereffect THM耦合从而影响迁移的废水。在本节中,污水的流入温度(表gydF4y2Ba 6gydF4y2Ba)被选中作为一个参数来研究污染物积累。gydF4y2Ba

废水流入的温度。gydF4y2Ba

模型没有。gydF4y2Ba 模型T1gydF4y2Ba 模型T2gydF4y2Ba 模型T3gydF4y2Ba
流入温度(°C)gydF4y2Ba 20.gydF4y2Ba 30.gydF4y2Ba 40gydF4y2Ba

沿线的污染物的积累gydF4y2Ba CgydF4y2Ba DgydF4y2Ba 在不同流入温度绘制在图gydF4y2Ba 17gydF4y2Ba。污染物的累积的gydF4y2Ba CgydF4y2Ba DgydF4y2Ba 随着时间的推移,逐渐增加。30年后,沿线的污染物的积累gydF4y2Ba CgydF4y2Ba DgydF4y2Ba 达到1.58摩尔(模型T1), 1.69摩尔(T2)模型,分别和1.76摩尔(T3)模型。流入水的温度有一个积极的积累与污染物的关系,虽然不影响断层参数。gydF4y2Ba

沿着线污染物积累gydF4y2Ba CgydF4y2Ba DgydF4y2Ba 在不同流入温度。gydF4y2Ba

 4.6。讨论gydF4y2Ba

从数据gydF4y2Ba 13gydF4y2Ba- - - - - -gydF4y2Ba 17gydF4y2Ba,我们可以推断,断层长度最显著影响污染物积累,而废水的温度最低影响污染物积累。gydF4y2Ba

同时,积累的污染物分为2个阶段,缓慢增长阶段,快速发展阶段。和所有的模型,缓慢增长阶段是前20年,20年后开始快速增长阶段。2阶段的原因可能是污染物运移的主要通道缓慢增长阶段是错;因此,浓度积累增加缓慢。在长期快速增长阶段,规模渗流导致整个污染物通过岩石的渗透矩阵以及断层;因此,浓度积累迅速增加。gydF4y2Ba

 5。结论gydF4y2Ba

控制方程的基础上THM耦合理论,提出了多媒体的MF和MFF。污染物的迁移进行了分析。参数研究断层的几何和废水的温度。结论可以概括如下:gydF4y2Ba

周围的压力分布强烈非齐次的错。断裂特征影响浓度分布,断层主导浓度分布和污染物迁移gydF4y2Ba

污染物的迁移可以几何分为3区沿着错误的方向,饱和区,快速扩散区,浓度降低gydF4y2Ba

有一个峰值浓度的模型的底部。这样的峰值出现在端点上的投影点的错,因为这个位置最近的故障gydF4y2Ba

污染物的积累呈正相关的断层长度、宽度、倾角、渗透率和温度的废水。断层长度对污染物积累最重要的影响。废水的温度最低影响污染物积累gydF4y2Ba

的积累浓度可分为2个阶段,缓慢增长阶段(前20年)和快速增长阶段(20年后)。污染物运移的主要通道缓慢增长阶段是一个错。在快速增长阶段,污染物穿透岩石矩阵以及断层gydF4y2Ba

符号gydF4y2Ba ηgydF4y2Ba wgydF4y2Ba :gydF4y2Ba

流体压缩系数gydF4y2Ba

 βgydF4y2Ba wgydF4y2Ba :gydF4y2Ba

液体的热膨胀系数gydF4y2Ba

 ρgydF4y2Ba :gydF4y2Ba

流体的参考密度(公斤/米gydF4y2Ba3gydF4y2Ba)gydF4y2Ba

 ρgydF4y2Ba 0gydF4y2Ba :gydF4y2Ba

废水的密度(公斤/米gydF4y2Ba3gydF4y2Ba)gydF4y2Ba

 ρgydF4y2Ba 年代gydF4y2Ba :gydF4y2Ba

岩体的密度(公斤/米gydF4y2Ba3gydF4y2Ba)gydF4y2Ba

 pgydF4y2Ba :gydF4y2Ba

孔隙压力(MPa)gydF4y2Ba

 pgydF4y2Ba 0gydF4y2Ba :gydF4y2Ba

初始孔隙压力(MPa)gydF4y2Ba

 TgydF4y2Ba :gydF4y2Ba

流体温度(°C)gydF4y2Ba

 TgydF4y2Ba 0gydF4y2Ba :gydF4y2Ba

初始温度(°C)gydF4y2Ba

 φgydF4y2Ba :gydF4y2Ba

的孔隙度gydF4y2Ba

 φgydF4y2Ba 0gydF4y2Ba :gydF4y2Ba

初始孔隙度gydF4y2Ba

 kgydF4y2Ba :gydF4y2Ba

矩阵的渗透率(mgydF4y2Ba2gydF4y2Ba)gydF4y2Ba

 kgydF4y2Ba 0gydF4y2Ba :gydF4y2Ba

的初始渗透矩阵(mgydF4y2Ba2gydF4y2Ba)gydF4y2Ba

 kgydF4y2Ba 1gydF4y2Ba :gydF4y2Ba

故障的初始渗透率(mgydF4y2Ba2gydF4y2Ba)gydF4y2Ba

 kgydF4y2Ba 2gydF4y2Ba :gydF4y2Ba

的初始渗透率骨折(mgydF4y2Ba2gydF4y2Ba)gydF4y2Ba

 ggydF4y2Ba :gydF4y2Ba

重力加速度(m / sgydF4y2Ba2gydF4y2Ba)gydF4y2Ba

 问gydF4y2Ba 米gydF4y2Ba :gydF4y2Ba

源项gydF4y2Ba

 αgydF4y2Ba TgydF4y2Ba :gydF4y2Ba

岩石的热膨胀系数矩阵gydF4y2Ba

 λgydF4y2Ba 年代gydF4y2Ba :gydF4y2Ba

岩石的导热系数矩阵(W / (m·°C))gydF4y2Ba

 λgydF4y2Ba wgydF4y2Ba :gydF4y2Ba

废水的导热系数(W / (m·°C))gydF4y2Ba

 cgydF4y2Ba wgydF4y2Ba :gydF4y2Ba

水的比热容(J /(公斤·°C))gydF4y2Ba

 cgydF4y2Ba 年代gydF4y2Ba :gydF4y2Ba

岩石的比热容矩阵(J /(公斤·°C))gydF4y2Ba

 KgydF4y2Ba :gydF4y2Ba

岩石的体积变形模量矩阵(Pa)gydF4y2Ba

 KgydF4y2Ba 年代gydF4y2Ba :gydF4y2Ba

岩石颗粒的体积变形模量(Pa)gydF4y2Ba

 EgydF4y2Ba :gydF4y2Ba

岩石的杨氏模量矩阵(Pa)gydF4y2Ba

 EgydF4y2Ba 年代gydF4y2Ba :gydF4y2Ba

杨氏模量的岩石颗粒(Pa)gydF4y2Ba

 υgydF4y2Ba :gydF4y2Ba

岩石的泊松比矩阵gydF4y2Ba

 εgydF4y2Ba vgydF4y2Ba :gydF4y2Ba

岩石的体积应变矩阵gydF4y2Ba

 dgydF4y2Ba 1gydF4y2Ba :gydF4y2Ba

断层的宽度(米)gydF4y2Ba

 dgydF4y2Ba 2gydF4y2Ba :gydF4y2Ba

骨折(m)的光圈。gydF4y2Ba

 数据可用性gydF4y2Ba

使用的数据来支持本研究的发现可以从第一作者。gydF4y2Ba

 的利益冲突gydF4y2Ba

作者宣称没有利益冲突。gydF4y2Ba

 确认gydF4y2Ba

这个项目是由中国国家自然科学基金(52078477和52078477号),深地球科学与工程重点实验室(四川大学),教育部(DESE202106和DESE202004),广东省重点实验室地球科学和地热能开发利用(深圳大学)(2020 - 3)。作者承认上述资助机构的支持。gydF4y2Ba

 陈gydF4y2Ba M。gydF4y2Ba 陈gydF4y2Ba h . Q。gydF4y2Ba 工业污水排放的时空耦合测量,在中国工业经济gydF4y2Ba 环境科学与污染研究国际gydF4y2Ba 2021年gydF4y2Ba 螳螂gydF4y2Ba 我。gydF4y2Ba VoutsagydF4y2Ba D。gydF4y2Ba 萨马拉gydF4y2Ba C。gydF4y2Ba 环境风险评估从市政和工业污水处理污泥采用化学和生物方法gydF4y2Ba 生态毒理学和环境安全gydF4y2Ba 2005年gydF4y2Ba 62年gydF4y2Ba 3gydF4y2Ba 397年gydF4y2Ba 407年gydF4y2Ba 10.1016 / j.ecoenv.2004.12.010gydF4y2Ba 2 - s2.0 - 26444581786gydF4y2Ba 陈gydF4y2Ba lgydF4y2Ba 郑gydF4y2Ba H。gydF4y2Ba 杨gydF4y2Ba H。gydF4y2Ba 规则设计和现状评价cross-industrial废水排放交易在中国典型的工业领域gydF4y2Ba 欧洲环境科学gydF4y2Ba 2020年gydF4y2Ba 32gydF4y2Ba 3gydF4y2Ba 7077年gydF4y2Ba 7084年gydF4y2Ba 梁gydF4y2Ba B。gydF4y2Ba 江gydF4y2Ba h·J。gydF4y2Ba 李gydF4y2Ba J·J。gydF4y2Ba 龚gydF4y2Ba C。gydF4y2Ba 骨折断裂参数的系统研究影响网络渗透率基于discrete-fracture模型运用有限元分析gydF4y2Ba 天然气的科学与工程》杂志上gydF4y2Ba 2016年gydF4y2Ba 28gydF4y2Ba 711年gydF4y2Ba 722年gydF4y2Ba 10.1016 / j.jngse.2015.12.011gydF4y2Ba 2 - s2.0 - 84954285036gydF4y2Ba 王gydF4y2Ba Y。gydF4y2Ba 易gydF4y2Ba y F。gydF4y2Ba 李gydF4y2Ba c . H。gydF4y2Ba 汉gydF4y2Ba j . Q。gydF4y2Ba 各向异性裂缝和能量特征的西藏大理石暴露于多层次的等幅循环加载(MLCA):一个实验室规模的测试gydF4y2Ba 工程断裂力学gydF4y2Ba 2021年gydF4y2Ba 244年gydF4y2Ba 107550年gydF4y2Ba 10.1016 / j.engfracmech.2021.107550gydF4y2Ba 王gydF4y2Ba Y。gydF4y2Ba 张gydF4y2Ba B。gydF4y2Ba 李gydF4y2Ba B。gydF4y2Ba 李gydF4y2Ba c . H。gydF4y2Ba 以疲劳损伤模型自然断裂的大理石受到冻融和单轴循环加载gydF4y2Ba 损伤力学的国际期刊gydF4y2Ba 2021年gydF4y2Ba 30.gydF4y2Ba 105678952110216gydF4y2Ba 10.1177 / 10567895211021629gydF4y2Ba 朱gydF4y2Ba h·G。gydF4y2Ba 谢gydF4y2Ba h·P。gydF4y2Ba 易gydF4y2Ba C。gydF4y2Ba 江gydF4y2Ba Y。gydF4y2Ba 刘gydF4y2Ba J。gydF4y2Ba 赖gydF4y2Ba 年代。gydF4y2Ba 越南盾gydF4y2Ba X。gydF4y2Ba 岩石裂缝中流体流动特性的分析gydF4y2Ba 中国岩石力学工程杂志》上gydF4y2Ba 2013年gydF4y2Ba 32gydF4y2Ba 4gydF4y2Ba 657年gydF4y2Ba 663年gydF4y2Ba 周gydF4y2Ba Y。gydF4y2Ba 拉贾帕克萨gydF4y2Ba r·k·n·D。gydF4y2Ba 格雷厄姆gydF4y2Ba J。gydF4y2Ba thermo-osmosis和thermal-filtration thermoporoelastic耦合模型gydF4y2Ba 国际期刊的固体和结构gydF4y2Ba 1998年gydF4y2Ba 35gydF4y2Ba 34-35gydF4y2Ba 4659年gydF4y2Ba 4683年gydF4y2Ba 10.1016 / s0020 - 7683 (98) 00089 - 4gydF4y2Ba 2 - s2.0 - 0032430718gydF4y2Ba 太阳系gydF4y2Ba j . M。gydF4y2Ba 耦合的影响对放射性核素迁移现象在Opalinus粘土和影响交通gydF4y2Ba 《污染物水文gydF4y2Ba 2001年gydF4y2Ba 53gydF4y2Ba 1 - 2gydF4y2Ba 63年gydF4y2Ba 84年gydF4y2Ba 10.1016 / s0169 - 7722 (01) 00140 - 1gydF4y2Ba 2 - s2.0 - 0035546237gydF4y2Ba 陈gydF4y2Ba x H。gydF4y2Ba PaogydF4y2Ba W。gydF4y2Ba 李gydF4y2Ba x K。gydF4y2Ba 耦合thermo-hydro-mechanical模型考虑thermal-osmosis混合物理论基础上修改gydF4y2Ba 国际工程科学杂志》上gydF4y2Ba 2013年gydF4y2Ba 64年gydF4y2Ba 1gydF4y2Ba 13gydF4y2Ba 10.1016 / j.ijengsci.2012.12.005gydF4y2Ba 2 - s2.0 - 84872965759gydF4y2Ba 蔡gydF4y2Ba g . Q。gydF4y2Ba 赵gydF4y2Ba c·G。gydF4y2Ba 田gydF4y2Ba H。gydF4y2Ba 数值模拟耦合thermo-hydro-mechanical行为设计壁垒高放射性核废料处理gydF4y2Ba 中国岩土工程杂志》上gydF4y2Ba 2013年gydF4y2Ba 35gydF4y2Ba 1gydF4y2Ba 1gydF4y2Ba 8gydF4y2Ba 陈gydF4y2Ba w Z。gydF4y2Ba 龚gydF4y2Ba Z。gydF4y2Ba 余gydF4y2Ba h . D。gydF4y2Ba Yong-shanggydF4y2Ba M。gydF4y2Ba Hong-minggydF4y2Ba T。gydF4y2Ba 回顾thermo-hydro-mechanical耦合测试和粘土的本构模型gydF4y2Ba 岩石和土力学gydF4y2Ba 2015年gydF4y2Ba 36gydF4y2Ba 5gydF4y2Ba 1217年gydF4y2Ba 1238年gydF4y2Ba 蔡gydF4y2Ba g . Q。gydF4y2Ba Yan-xingydF4y2Ba G。gydF4y2Ba 江ydF4y2Ba lgydF4y2Ba Xue-donggydF4y2Ba Z。gydF4y2Ba )gydF4y2Ba Z。gydF4y2Ba 弹塑性模型的体积变化不饱和土壤thermo-hydro-mechanical耦合条件下的行为gydF4y2Ba 岩石和土力学gydF4y2Ba 2017年gydF4y2Ba 38gydF4y2Ba 4gydF4y2Ba 1060年gydF4y2Ba 1068年gydF4y2Ba 张gydF4y2Ba z H。gydF4y2Ba 史gydF4y2Ba y . M。gydF4y2Ba 朱gydF4y2Ba M。gydF4y2Ba 耦合hydro-mechanical-chemical模型粘土衬垫gydF4y2Ba 中国岩土工程杂志》上gydF4y2Ba 2016年gydF4y2Ba 38gydF4y2Ba 7gydF4y2Ba 1283年gydF4y2Ba 1290年gydF4y2Ba SeetharamgydF4y2Ba s . C。gydF4y2Ba 托马斯。gydF4y2Ba h·R。gydF4y2Ba CleallgydF4y2Ba p . J。gydF4y2Ba 耦合热/水电/化学/力学模型不饱和土壤:数值算法gydF4y2Ba 国际期刊工程中的数值方法gydF4y2Ba 2007年gydF4y2Ba 70年gydF4y2Ba 12gydF4y2Ba 1480年gydF4y2Ba 1511年gydF4y2Ba 10.1002 / nme.1934gydF4y2Ba 2 - s2.0 - 34250357568gydF4y2Ba CleallgydF4y2Ba p . J。gydF4y2Ba SeetharamgydF4y2Ba s . C。gydF4y2Ba 托马斯。gydF4y2Ba h·R。gydF4y2Ba 包含不饱和土壤的化学行为的某些方面在热/水电/化学/力学模型gydF4y2Ba 《工程力学gydF4y2Ba 2007年gydF4y2Ba 133年gydF4y2Ba 3gydF4y2Ba 337年gydF4y2Ba 347年gydF4y2Ba 一族gydF4y2Ba 一个。gydF4y2Ba 吉马良斯gydF4y2Ba l . d . N。gydF4y2Ba OlivellagydF4y2Ba 年代。gydF4y2Ba 桑切斯gydF4y2Ba M。gydF4y2Ba 造型thermo-hydro-mechano-chemical核废物处理的交互gydF4y2Ba 《岩石力学与岩土工程gydF4y2Ba 2010年gydF4y2Ba 2gydF4y2Ba 2gydF4y2Ba 97年gydF4y2Ba 102年gydF4y2Ba 10.3724 / SP.J.1235.2010.00097gydF4y2Ba PreshogydF4y2Ba M。gydF4y2Ba 我们gydF4y2Ba 年代。gydF4y2Ba 方面能gydF4y2Ba V。gydF4y2Ba 校准双重孔隙度、裂缝储层的双重渗透率建模gydF4y2Ba 石油科学与工程》杂志上gydF4y2Ba 2011年gydF4y2Ba 77年gydF4y2Ba 3 - 4gydF4y2Ba 326年gydF4y2Ba 337年gydF4y2Ba 10.1016 / j.petrol.2011.04.007gydF4y2Ba 2 - s2.0 - 79957445946gydF4y2Ba 陆gydF4y2Ba M . M。gydF4y2Ba 谢gydF4y2Ba k . H。gydF4y2Ba 王gydF4y2Ba s Y。gydF4y2Ba 整合的垂直排水管depth-varying加载应力引起的多阶段gydF4y2Ba 电脑和土工技术gydF4y2Ba 2011年gydF4y2Ba 38gydF4y2Ba 8gydF4y2Ba 1096年gydF4y2Ba 1101年gydF4y2Ba 10.1016 / j.compgeo.2011.06.007gydF4y2Ba 2 - s2.0 - 80052930236gydF4y2Ba 雷森gydF4y2Ba r . J。gydF4y2Ba 模拟天然裂缝性储层与semi-implicit源项gydF4y2Ba 石油工程师学会》杂志上gydF4y2Ba 1977年gydF4y2Ba 17gydF4y2Ba 3gydF4y2Ba 201年gydF4y2Ba 210年gydF4y2Ba 10.2118 / 5737 - pagydF4y2Ba 雪gydF4y2Ba d . T。gydF4y2Ba 岩石裂缝间距、开口、疏密度gydF4y2Ba 土木土力学和基础部门日报》gydF4y2Ba 1968年gydF4y2Ba 94年gydF4y2Ba 1gydF4y2Ba 73年gydF4y2Ba 91年gydF4y2Ba 10.1061 / JSFEAQ.0001097gydF4y2Ba 谢长廷gydF4y2Ba P。gydF4y2Ba 纽曼gydF4y2Ba s P。gydF4y2Ba 辛普森gydF4y2Ba 大肠。gydF4y2Ba 斯泰尔斯gydF4y2Ba G。gydF4y2Ba 场的确定三维各向异性介质的渗透系数张量:1。理论gydF4y2Ba 水资源研究gydF4y2Ba 1985年gydF4y2Ba 21gydF4y2Ba 11gydF4y2Ba 1655年gydF4y2Ba 1665年gydF4y2Ba 10.1029 / WR021i011p01655gydF4y2Ba 2 - s2.0 - 0022264617gydF4y2Ba KuwaharagydF4y2Ba F。gydF4y2Ba 孝宏gydF4y2Ba U。gydF4y2Ba 中山gydF4y2Ba 一个。gydF4y2Ba 宏观动量方程多孔介质流动的双重结构gydF4y2Ba 日本化学学会gydF4y2Ba 2000年gydF4y2Ba 26gydF4y2Ba 6gydF4y2Ba 837年gydF4y2Ba 841年gydF4y2Ba 10.1252 / kakoronbunshu.26.837gydF4y2Ba 黄ydF4y2Ba x L。gydF4y2Ba 徐gydF4y2Ba G。gydF4y2Ba 张gydF4y2Ba Y。gydF4y2Ba 太阳gydF4y2Ba F。gydF4y2Ba 雪gydF4y2Ba 年代。gydF4y2Ba 研究thermo-hydro-mechanical耦合和地热井筒结构的稳定性gydF4y2Ba 能量gydF4y2Ba 2021年gydF4y2Ba 14gydF4y2Ba 3gydF4y2Ba 649年gydF4y2Ba 649年gydF4y2Ba 10.3390 / en14030649gydF4y2Ba 史gydF4y2Ba b P。gydF4y2Ba 王gydF4y2Ba y . C。gydF4y2Ba 陆gydF4y2Ba b D。gydF4y2Ba 赵gydF4y2Ba X X。gydF4y2Ba 研究新的液体的粘度和温度之间的关系gydF4y2Ba 太原大学科学和技术杂志》上gydF4y2Ba 2007年gydF4y2Ba 28gydF4y2Ba 6gydF4y2Ba 500年gydF4y2Ba 503年gydF4y2Ba 陈gydF4y2Ba b G。gydF4y2Ba 研究数值方法耦合流体流动和传热的破碎岩石的双重系统gydF4y2Ba 2014年gydF4y2Ba 清华大学gydF4y2Ba 壮族gydF4y2Ba x Y。gydF4y2Ba 黄gydF4y2Ba R。gydF4y2Ba 梁gydF4y2Ba C。gydF4y2Ba RabczukgydF4y2Ba T。gydF4y2Ba 耦合thermo-hydro-mechanical模型的压缩空气储能贴合硬摇滚gydF4y2Ba 数学问题在工程gydF4y2Ba 2014年gydF4y2Ba 2014年gydF4y2Ba 11gydF4y2Ba 179169年gydF4y2Ba 10.1155 / 2014/179169gydF4y2Ba 2 - s2.0 - 84893825113gydF4y2Ba 高gydF4y2Ba K。gydF4y2Ba THM耦合模型干热岩裂缝性储层的开采及其应用gydF4y2Ba 2020年gydF4y2Ba 中国矿业大学和技术gydF4y2Ba 徐gydF4y2Ba P。gydF4y2Ba 余gydF4y2Ba b . M。gydF4y2Ba 开发一种新形式的渗透率和Kozeny-Carman常数为均匀多孔介质分形几何的方法gydF4y2Ba 水资源的进步gydF4y2Ba 2008年gydF4y2Ba 31日gydF4y2Ba 1gydF4y2Ba 74年gydF4y2Ba 81年gydF4y2Ba 10.1016 / j.advwatres.2007.06.003gydF4y2Ba 2 - s2.0 - 36048947628gydF4y2Ba 把握现在gydF4y2Ba y D。gydF4y2Ba 周gydF4y2Ba x P。gydF4y2Ba 一个耦合的振动非平常状态peridynamics裂缝性多孔岩石gydF4y2Ba 工程分析与边界元素gydF4y2Ba 2021年gydF4y2Ba 123年gydF4y2Ba 133年gydF4y2Ba 146年gydF4y2Ba 10.1016 / j.enganabound.2020.12.001gydF4y2Ba 路易gydF4y2Ba C。gydF4y2Ba 岩石水力学在岩石的机械师gydF4y2Ba 1974年gydF4y2Ba 纽约gydF4y2Ba 施普林格-gydF4y2Ba 吴gydF4y2Ba J。gydF4y2Ba 水力断裂表面分形的特征gydF4y2Ba 2010年gydF4y2Ba 大连理工大学gydF4y2Ba