关键词： 页岩气   非理想气体效应   流固耦合   油气地质  

摘要： 页岩气是未来中国能源结构中的重要组成部分,但中国页岩气埋藏深度大,不能直接照搬美国对浅层页岩气的开采经验,需要对其输运机理有更清晰的认识.介绍了跨尺度混合模拟算法,分析了基于孔隙尺度模拟来研究场尺度的气井衰减曲线分析等问题,揭示了高努森数效应与非理想气体效应之间的耦合机理;介绍了考虑多孔介质弹性形变对渗流影响时计算表观渗透率的特征压力模型,为页岩气勘探开发提供了相关理论支持.

关键词： 数值模拟   中观尺度   裂隙   波的衰减   流体压力  

摘要： 上地壳岩石普遍为含有一定孔隙度的复杂孔隙结构,孔隙填充物可由粘土、油、气、水等成份组成,这种不均匀异质性特征使关于地震波在地下介质中传播问题的研究更为复杂。在力源的作用下,岩石颗粒会发生应力应变状态的改变,从而可以间接的反映出传输介质的多参数物理特征,因而关于地震频带范围内弹性波传播规律的研究可以更有效的为储层地球物理以及工程、环境等方面服务。关于复杂介质结构及成份的研究是讨论波传播问题的基础,对于多孔介质而言,裂隙及裂纹通常是代表性的孔隙结构,本文也依此建立了双层含裂隙孔隙模型。为了更符合实际,我们同时考虑到介质的粘弹特性和各向异性特征,通过改变上下两层介质物理参数,如裂隙密度及弹性模量等,对比波在不同物理性质中的传播规律,同时得到变量参数对试验结果的影响。地震波在不同类型介质中的传播规律具有一定的指示作用,因而被作为地震勘探的主要手段,由于地层的吸收作用,波在传播时会发生一定程度的散射和损耗,尤其可以应用到岩石饱和度分析及储层质量控制等方面。近年来,有关波在中观尺度上衰减现象的研究表明,波在中观尺度(介于波长和粒子直径之间,远小于波长)上诱导流体流(Wave Induced Fluid Flow,WIFF)的是造成地震波在复杂异质性介质中发生损耗的主要原因之一,因而使关于中观尺度地震波传播状态及相关参数影响的研究变得更有意义。在本文中,我们应用改进的高阶交错网格有限差分方法模拟了地震波在复杂异质性含裂隙孔隙介质中的传播状态,通过改变数值模拟参数及宽带雷克子波震源的应用,我们得到宏观尺度及中观尺度下的地震波场数据,作为Biot理论在孔隙介质中研究的基本结论,快、慢纵波以及横波均在试验模拟中得到实现,在模型参数的基础上,进一步计算了三种波的衰减特性,同时讨论裂隙参数及孔隙度的影响。针对介质的的不均匀性引起的横波分裂特征,我们进一步详细的讨论了裂隙等参数与横波分裂系数的关系。基于Darcy定律的条件下,同时设计含有扁平椭圆裂隙的中观尺度孔隙介质模型,并利用有限元模拟方法得到流体压力场的变化,同时对裂隙角度在其中的影响进行了讨论。在本文的研究中实现了不同尺度波场的数值模拟,并得到三类典型波的衰减特征,针对含裂隙孔隙介质中裂隙参数的影响,分别从横波分裂系数及流体压力场等几个方面进行讨论,以裂隙密度,裂隙旋转角度以及介质孔隙度等参数为主要研究对象,为进一步研究地震波在裂隙孔隙介质中的传播等问题提供基础。

关键词： 皮革和毛皮   多尺度孔隙   甲醛释放   暴露点浓度   传质模型  

摘要： 甲醛是一种无色的气体、刺激性强烈,严重的可致人头痛、失明、呼吸困难,甚至诱发癌症;作为基础的化工原料,在皮革化学品的制备和毛皮加工生产中有一定应用,且成本低廉,在短时间内不可能完全被替代。甲醛在皮革、毛皮中的存在方式比较复杂,包括:游离的、毛细管吸附的、可逆和不可逆键合的甲醛分子。前三种结合方式的甲醛在皮革、毛皮使用或存放过程中一旦释放出来并达到一定量时,会对人体和环境造成危害。因此包括我国在内的很多国家对皮革、毛皮中游离甲醛含量提出严格的限量标准。针对皮革、毛皮中甲醛的研究,目前主要集中在基于水萃取体系的检测方法、少醛合成鞣剂、甲醛捕获剂等方面,但是对于具有多尺度孔隙结构的皮革、毛皮中甲醛的释放行为的研究,无论国内还是国外均鲜有文献报道。皮革、毛皮中甲醛在微观尺度、介观尺度下的吸附/解吸附行为,随着使用、存放介质状态影响下的宏观释放行为,在孔隙中扩散方式与孔隙结构的关系等均不明确。本论文针对上述问题开展了研究,具体研究工作包括以下几个方面:（1）系统开展了我国主要特色产地皮革、毛皮中甲醛含量评估工作。收集并总结了我国皮革、毛皮重镇—桐乡地区2014～2016年度皮革、毛皮中甲醛含量的检测数据。基于现有检测数据,采用Bootstrap方法进行模拟重新采样,建立每个Bootstrap样本数据集的均值分布,并对每年样本的95%UCL（Upper Confidence Limit,置信上限）进行了计算和比较。通过Wilcoxon秩和检验,检验了2014～2016年样本数据的位置参数,获得了皮革、毛皮中甲醛含量的变化趋势。Bootstrap方法克服了非正态分布和样本量小的缺点,可以广泛应用于未知数据分布形态,为检测数据的特征分析提供了一种新的方法。研究结果表明:桐乡地区2014～2016年度的皮革、毛皮中甲醛含量分布接近正态分布,但高度偏斜。三年样本均值的单边95%UCL分别为95.45 mg/kg、93.02 mg/kg、85.64 mg/kg,逐年下降。Wilcoxon秩和检验表明,2016年的甲醛含量的总体分布与前两年相比有明显的结构性改进,产品质量在不断提升。但仍然高于GB20400-2006中规定的B类标准:直接接触皮肤的产品≤75 mg/kg。因此,皮革、毛皮中甲醛环境风险依然存在,开展甲醛释放及防治等相关问题的研究具有重要的理论研究意义和实际应用价值。（2）以铬皮粉作为皮革模拟物,开展了与甲醛吸附/解吸附行为的研究,为宏观释放行为提供理论基础。依据醛鞣工艺设计了吸附试验,研究了甲醛初始浓度、提碱pH、铬皮粉用量、吸附温度、吸附时间对吸附率和吸附量的影响;并基于试验数据,借助数学拟合方法研究了铬皮粉对甲醛的吸附平衡及动力学。并对吸附后的样品进行了解吸附试验,考察了时间、温度和样品pH对解吸行为的影响。研究结果表明:铬皮粉对甲醛的吸附等温线符合Langmuir方程,吸附动力学符合拟二级速率方程。以吸附状态来看,铬皮粉对水体溶液中甲醛的吸附是以化学吸附为主,兼有物理吸附。温度、pH对甲醛的吸附/解吸行为影响显著。温度可以促进吸附/解吸附过程。随着pH的升高,甲醛的吸附率不断提高、解吸率不断降低。（3）基于皮革、毛皮穿着时的介质条件,使用密封瓶模拟甲醛释放环境条件,通过经验模型开展了释放行为研究。对甲醛鞣制羊皮革、甲醛鞣制牛皮革以及改性戊二醛鞣制羊皮革进行甲醛释放试验,分别考察了温度、样品吸水率、吸收液体积等因素对释放行为的影响;其次基于释放模型对释放曲线方程进行了理论推导,并通过实际释放数据拟合对释放理论进行了验证;最后考察了不同样品平衡释放率。研究结果表明:密封系统中,皮革中甲醛释放可以划分为三个释放阶段:快速释放阶段、缓慢释放阶段、平衡阶段。线性释放时间只是温度的函数,温度升高,传质系数增大,线性释放时间缩短;温度降低,传质系数减小,线性释放时间延长。线性释放阶段内,温度对甲醛释放速度影响显著,温度越高甲醛释放越快,但基于孔隙差异,不同皮种中甲醛的释放对温度敏感性不同。试样吸水率随着时间线性增加,温度提升可以促进试样吸水,试样吸水率与甲醛释放量呈正线性相关,吸水率的提升可以促进甲醛的释放,温度和湿度协同作用,释放效果更明显。密封系统中,即使是在高温度（60℃）且有水存在的条件下,甲醛鞣制的样品达到释放平衡时其释放率仅为20%左右,而改性戊二醛鞣制的样品,40℃时,其平衡释放率可达43.60%。甲醛鞣制的样品中剩余的甲醛会在环境因素改变后,持续释放,其危害是长期性的。（4）基于皮革、毛皮非穿着时的介质条件,使用环境舱模拟甲醛释放环境条件,通过传质模型开展了释放行为研究。将甲醛鞣制皮革样品置于环境舱中模拟释放,以空气介质中甲醛含量为考察指标,研究了环境温度、相对湿度、换气速率、承载率等对释放行为的影响规律。基于传质理论,建立了

摘要： Implementations of three models were used to constrain permeability and vertical bulk modulus of oceanic crust using pressure data collected from CORKs at ODP/IODP Holes 1024C, 1025C, 1026B, 1027C, and 504B. The models are based on a one-dimensional half-space model, which interprets a given sinusoidal signal within formation as a combination of two signals produced by ocean tidal loading: an elastic signal at the measurement site that is in-phase with tidal loading, and a diffusive signal that travels through the formation to the instrument site. The aforementioned sites were chosen due to their thick sedimentation: the sediment layer is too thick to allow a diffusive signal through, so diffusive signals can only be generated at thinly sedimented or exposed basement and be transmitted laterally through basement. Harmonic analyses were applied to seafloor and formation pressure data gathered from the CORKs to determine the tidal loading signal at the seafloor and the tidal signal within the formation fluid; the signals analyzed were at the S2, M2, K1, and O1 frequencies. Two implementations of the one-dimensional half-space model were used: a brute force method, and a least squares method. These were used compared to each other and to previous studies; the brute force method was found to be better overall, though the least squares method agreed with the brute force method in constraining the elastic signal. A two-dimensional isotropic model was used to constrain the “bulk” or “isotropic equivalent” values of permeability; a two-dimensional anisotropic model was used to constrain the ridge-parallel and ridge-perpendicular components of permeability. All three models show a decrease in permeability with increasing crustal age, with a drop of one to two orders of magnitude ([10-11]-[10 -10] m2 to [10-12]-[10-11] m2) from roughly 1 Ma to 6 Ma. Overall, permeability showed a faster decrease with diurnal components versus semidiurnal components. The anisotropic model sugg

关键词： CO2盐水层封存   孔隙尺度   相场法   非均质性   数值模拟  

摘要： 温室效应导致自然灾害频发,严重影响人类社会的生产与生活。CO2是温室气体的主要成分,减少CO2排放成为缓解温室效应的首要目标。CO2盐水层封存是减少CO2向大气排放最有效与最具有发展潜力的方法之一。理解CO2在盐水层内的驱替过程、CO2溶于盐水引发的自然对流过程对评估该方法安全性和有效性具有重要意义。储存结构非均质性决定了CO2/盐水运移过程是一个复杂的水文地质过程。因此,研究结构非均质性对CO2盐水驱替以及对流混合过程的影响具有重要价值。多孔介质非均质性是影响CO2/盐水驱替的重要因素。相场法可以有效实现对多孔介质内CO2/盐水相界面的捕捉。本文运用相场法建立数值模型,研究孔隙结构非均质性对CO2/盐水两相驱替的影响。研究表明:1)CO2饱和度与CO2指进的形成、发育均受多孔介质非均质性影响;2)胶结物的存在改变了多孔介质的连通性,能够引起绕流,促使指进产生,增加驱替的不稳定性;3)死孔导致残余水饱和度增高,CO2饱和度降低;4)在润湿性非均质（亲水性）介质内CO2指进被抑制,导致CO2饱和度急剧增加。对流混合可以有效促进CO2的溶解速度。流体动力学与Darcy渗流描述的是不同尺度的现象,前者适用于孔隙尺度描述了孔隙内流体间的相互作用,后者适用于岩心及更大尺度描述了多孔介质内流体平均流速与多孔介质宏观物性的关系。本文在孔隙尺度下对初始时间步长及岩心尺度下对网格数量做敏感度分析,对比孔隙-岩心尺度模型,孔隙结构非均质性对对流混合的影响。研究表明:1)初始时间步长与网格数量对模拟有一定影响,都会导致对流强度与指进触底时间的改变。2)与孔隙尺度相比,岩心尺度下对流混合不稳定积累时间长,指进触底时间晚,相同时间下CO2溶解量少;3)与均质孔隙结构相比,局部孔喉减小、胶结物存在都会延长指进触底时间,局部孔喉增大会缩短指进触底时间。

关键词： 格子Boltzmann方法   孔隙尺度   渗流   传质   化学反应  

摘要： 近些年来,可再生能源开发利用的规模不断扩大,储能技术的研究与应用日益广泛。不同的储能技术都有其各自的优缺点,氧化还原液流电池被认为是将成为未来储能技术的发展趋势。氧化还原液流电池简称液流电池,根据其有无离子交换膜的结构特点可分为双液流电池和单液流电池。锌镍单液流电池就是一种无须离子交换膜的单液流电池,由防化研究院于2007年提出,它以烧结氧化镍为电池的正极,以锌为电池的负极,电解液为碱性锌酸盐溶液。锌镍单液流电池具有成本低、结构简单、循环寿命长等优点,可以为大规模储能提供一个高性价比的解决方案。本文利用格子Boltzmann方法探究了烧结镍多孔电极内渗流、传质与电化学反应的耦合问题,本文的主要研究工作如下:1)根据QSGS多孔介质重构方法的构造特点,分析其在构造一些特定结构的多孔介质时存在的不适用性原因,围绕锌镍单液流电池多孔烧结镍正极有针对性的提出了两种改进思路与方法,形成了基于改进四参数随机生长法的可控型多孔介质虚拟重构法,并利用改进的方法重构了锌镍单液流电池正极烧结镍的数字多孔介质模型。2)在锌镍单液流电池单个通道正极附近划定合理的计算区域,结合控制方程模型,利用格子Boltzmann方法求解了计算区域内的流场,通过改变计算区域的进口流速、流道宽度以及多孔介质骨架形态等条件,来探究其对烧结镍正极内部渗流的影响规律。3)提出一种描述烧结镍电极内部传质和化学反应的孔隙尺度模型,结合格子Boltzmann方法进行模拟计算,模拟了电池充放电过程中,正极荷电状态在50%时的稳态反应过程,得到了孔隙内OH-浓度、电流密度以及固相质子浓度的分布规律;探究了进口流速和放电电流密度对电极内传质和化学反应的影响。基于上述工作,作出了一些前瞻性的成果,证明了研究方法的可行性和技术路线的合理性,为进一步更加细致的探究电极电化学反应机理奠定基础。

关键词： 矿堆散体表征   双尺度孔隙结构   孔隙尺度浸矿模型   渗流一传质—反应—传热多场耦合   浸矿仿真  

摘要： 论文结合μCT扫描、计算机三维图像处理方法和浸矿过程建模仿真等手段,研究堆浸体系宏细观结构特征和水力学特性,提出矿堆结构精细化表征参数体系,基于矿堆三维图像,建立两级渗流—两级传质—固液反应—多相传热的浸矿多场耦合模型,通过氧化铜矿浸出的数值仿真,揭示孔隙尺度浸矿体系各物理场动态响应机制,主要研究工作如下:(1)通过MCT扫描获取浸柱和矿石三维图像,基于三维图像处理方法,分析矿石单颗粒和散体形态特性,形成以粒级—球度—颗粒排列—颗粒内矿物分布的矿石颗粒散体表征体系和以孔喉轴线长—断面面积—断面周长的孔隙表征体系,并以工艺矿物学分析和常规筛分实验结果部分验证图像计算的有效性。(2)提出两种基于图像计算的矿堆水力学特性预测方法:通过浸柱CT图像的颗粒散体特性分析,利用Carman-Kozeny经验方程计算介质渗透率;采用直接孔隙图像开展Navier-Stokes流动模拟计算任意(孔隙型/裂隙型)多孔介质渗透率/水力传导率。开展浸柱渗流实验,验证图像预测渗透率的可靠性。根据孔喉图像计算孔喉断面形状系数(圆度—粗糙度—径长比)和无量纲水力传导率,利用反向传播(BP)神经网络,建立形状系数与无量纲水力传导率的新关系。(3)在阐明宏细观结构和物理特性差异的基础上,根据矿堆三维CT图像构建双尺度孔隙结构矿堆,建立基于Navier-Stokes和Brinkman's双区域控制方程的矿堆渗流模型,开展两级渗流模拟,实现优势流可视化,分析细观体系重要性。(4)在两级渗流模型基础上建立包含宏观孔隙、脉石颗粒和孔雀石的双尺度孔隙结构矿堆浸矿模型。构造宏观孔隙和脉石颗粒内的两级渗流和对流—扩散传质控制方程,宏观孔隙、脉石颗粒和孔雀石中的对流—热传导多相传热控制方程,溶浸液与脉石矿物反应的区域反应控制方程,以及溶浸液与孔雀石反应的界面反应控制方程,分析渗流、传质、反应和传热过程的耦合关系,建立浸矿多物理场数学模型。(5)利用双尺度孔隙结构矿堆模型开展氧化铜矿浸出过程仿真,分析了浸矿体系离子(H+、Cu2+、Ca/Mg2+)浓度、矿堆温度场以及钙镁脉石和矿石矿物浸出程度的整体变化趋势,再现了矿堆内流场、离子浓度场、矿堆温度场以及局部钙镁脉石和矿石矿物(未)浸出率的分布;通过比较不同浸矿条件下的浸出效果,分析各因素对浸矿过程的影响。开展小型柱浸实验验证双尺度孔隙结构矿堆浸矿模型,监测浸出液H+浓度、铜及钙镁浸出率变化,并与仿真结果比较,分析实验与仿真差异及其原因。

关键词： 瓦斯流动   吸附解吸   菲克扩散   密度差驱动流   多尺度机理   双重孔隙  

摘要： 煤层瓦斯防治与利用是关乎煤矿安全开采,清洁能源高效利用以及温室气体排放三个方面的重要问题,一直是能源开采和环境保护等领域关注的热点。煤层瓦斯抽采、煤层气开采是防治瓦斯事故,捕获瓦斯能源以及控制温室气体排放的关键技术。煤体瓦斯流动理论是人们对煤中瓦斯存储、瓦斯解吸,瓦斯流动和瓦斯产出过程的理解和认识,是进行掌握瓦斯储量,评估瓦斯灾害风险,计算瓦斯涌出并制定防治措施,预测瓦斯产能,设计瓦斯资源开发等实际工程问题的基础和依据。本研究基于前人的研究成果,在煤体瓦斯流动理论的几个方面进行了以下探讨和尝试:煤层瓦斯含量计算方法,煤基质瓦斯吸附解吸机理及理论模型,煤层瓦斯多尺度流动机理及理论模型构建,数值离散及解算模拟程序开发等。本研究首先探讨了煤层瓦斯流动中关于瓦斯含量的计算问题,指出一般应采用朗格缪尔方程表示瓦斯含量,而当前仍有许多工程计算和学术文献中对瓦斯含量计算的简化(抛物线方程近似或P.M.克里钦夫斯基近似)。通过数值计算和理论计算对三种瓦斯含量计算方法进行了对比,定量分析了简化计算所造成的误差,发现:将朗格缪尔方程的瓦斯含量曲线近似为抛物线方程,由此所带来的误差较小,但采用P.M.克里钦夫斯基引起的误差是不可接受的;精确计算瓦斯流动时,应根据朗格缪尔式进行计算瓦斯含量。煤体是由裂隙和煤基质构成的典型的双重孔隙介质,煤体裂隙内的瓦斯流动机理已被广泛接受,而煤粒/煤基质内的瓦斯流动机理仍存有争议。本文通过研究当前有关煤粒/煤基质瓦斯流动理论的学术文献,认为煤中孔隙结构复杂且孔径分布广泛,气体流动应是多尺度多机理现象,仅用菲克扩散去描述多尺度孔隙中瓦斯运移是值得怀疑的。在现有理论的基础上,提出了密度差驱动流的概念,将煤体孔隙内的吸附瓦斯和游离瓦斯一并计算为瓦斯密度,将孔隙中的渗透和扩散经过传质概念的整合,统一为密度差驱动流,提出新的瓦斯吸附解吸理论模型。通过恒定环境和变压环境下的多尺度煤粒瓦斯吸附解吸实验,验证了本文提出的理论模型在描述煤粒/煤基质中瓦斯运移方面的有效性和优越性。通过研究现有的煤体瓦斯流动模型及数值模拟软件,认为有关煤基质瓦斯运移的模拟技术和数值解算方法有待改进。基于煤中瓦斯多尺度流动机理和双重孔隙模型,考虑煤基质收缩,有效应力和瓦斯多尺度流动的耦合机制,构建双重孔隙煤体瓦斯多尺度流动数学模型。该模型相较于其他双重孔隙模型的特点在于:考虑了煤基质内的渗流,扩散,将两种流动统一为密度差驱动流,并且充分考虑基质内的游离瓦斯和吸附瓦斯,与实际煤基质的瓦斯流动更为符合。通过有限差分法对双重孔隙煤体瓦斯多尺度流动数学模型进行离散,并基于Visual Basic编程平台独立开发了高效的数值模拟软件。该数值求解的优势在于:基质系统和裂隙系统不再相互重叠,流场中的任一点只有一套参数,基质参数或者裂隙参数,并且基质和裂隙之间的质量交换仅发生在基质块的表面,而不是作为在裂隙网络中均匀分布的质量源。这一处理有别于现存的其他双重孔隙模型的解算,与实际流场更为接近。通过将模拟结果与两套现场瓦斯抽放数据进行对比,验证了本章所建瓦斯流动模型及数值模拟软件的正确性。基于自主开发的煤体瓦斯流动数值模拟软件,进行了煤中气体流动和钻孔瓦斯抽采的一系列数值模拟,研究结果表明:煤体裂隙的初始渗透率越高,钻孔瓦斯流速越大。在抽采初期,改变裂隙渗透率对提高钻孔瓦斯产出效率有重要作用,但对抽采后期瓦斯产出效率影响不大。煤基质的渗透性对抽采初期的瓦斯产气效率影响较小;煤基质渗透率越大,抽采后期瓦斯产出效率越高,抽采过程中产气速率衰减越慢;煤基质的尺度越小,煤基质的释放瓦斯的速度越快,基质的累积瓦斯释放量越多。煤体中基质尺度越小,瓦斯产出速率越大,产气速率的衰减越慢。研究结果表明煤基质的尺度和渗透性是影响基质内气体运移的两个重要参数,也是影响抽采后期瓦斯产率的关键因素。同时表明瓦斯产出过程中产气机理随开采时间不断变化:在产气初期,裂隙中的达西流主导产气速率,在产气后期,基质内的密度差驱动流主导产气速率。

关键词： 孔隙尺度模拟   反常弥散   多孔介质   胶结度   分选度   死端孔隙  

摘要： 近年来地下水污染问题日益突出,且由于由固体颗粒堆积构成的多孔介质连通性好,污染物更容易进入其中造成地下水污染。长期以来,传统对流弥散模型(ADE模型)被用来模拟多孔介质中溶质运移过程。传统ADE模型将多孔介质中溶质运移过程描述为溶质羽的质量中心以平均流体流速沿主流动方向匀速迁移,而溶质羽在其质量中心周围的弥散行为可以用Fick's定律刻画。在该描述下溶质弥散系数为常数,溶质羽空间一阶矩、二阶矩均随时间线性增大。然而,越来越多的试验及模拟结果表明,多孔介质中溶质运移过程表现出很多系统偏离传统ADE模型预测结果的现象。如穿透曲线BTC与传统ADE模型预测的结果相比,常常存在首次穿透提前、完全穿透滞后的现象,即早穿透和后期长时间拖尾现象。另外,溶质弥散系数常常随运移时间增大而不是传统ADE模型预测的为常数,溶质羽空间一阶矩、二阶矩常常随运移时间非线性增大,这种现象被称为尺度效应。表现出以上早穿透、拖尾及尺度效应现象的溶质运移过程称为非菲克运移或反常运移。在很多实际问题中,需要精确预测多孔介质中溶质运移早期及后期穿透行为。另一方面,深刻理解多孔介质中的尺度效应,有助于提高多孔介质溶质运移预测精度。因此,对多孔介质溶质反常运移行为的研究具有重要的理论意义和实践价值。已有研究表明,多孔介质微观结构对溶质反常运移行为有较大影响,固体颗粒胶结度、分选度以及多孔介质中广泛存在的死端孔隙,对多孔介质孔隙结构的几何特征有显著影响。然而,针对不同胶结度和分选度的多孔介质中溶质反常运移行为特征及其内在机理、死端孔隙对溶质运移影响的研究还很少。本文运用数值模拟方法,在孔隙尺度上研究了固体颗粒胶结度、分选度不同的多孔介质中溶质运移行为特征及其内在机理。同时,研究了死端孔隙中,溶质滞留时间的概率分布形式及其影响因素,探讨了死端孔隙对溶质运移的影响。本文主要得到了以下结论:(1)不同胶结度的二维、三维多孔介质中孔隙尺度上的溶质运移模拟结果均表明,当胶结度增大时,溶质运移早穿透、后期拖尾现象均有所增强,其中拖尾现象的增强更加显著。同时,穿透曲线BTC越来越偏离传统ADE模型预测的结果,溶质运移过程更加反常。三维多孔介质中溶质羽空间矩及传递函数模拟结果表明,当胶结度不同时,溶质运移过程中溶质羽空间分布随时间的演化行为也完全不同。胶结度大的多孔介质中,空间二阶矩不再随时间线性增大,而是以幂函数的形式增大。尺度效应明显增强。另外,随着胶结度增大,溶质运移传递函数中逐渐出现停滞浓度峰,且胶结程度越大,停滞峰峰值越大、持续时间越长。(2)不同胶结度的多孔介质中流场模拟结果表明,当胶结度增大时,二维、三维多孔介质中流场特征表现出相似的变化趋势。随着胶结度增大,流场非均质性明显增强,流速概率密度函数分布更加发散,且不流动区域、优势流区域明显增多。溶质运移连通系数CT及传递函数随时间的演化行为表明,正是以上流场特征变化导致了不同的溶质运移行为。(3)在由不同粒径的两种固体颗粒掺杂、堆积构成的多孔介质中,当Rb从0.0逐渐增大到1.0的过程时,溶质运移反常程度存在一个最大值和一个最小值。且当固体颗粒粒径、流体平均流速、溶质扩散系数改变时,对应于反常程度最大和最小的Rb值即Rbmax、Rbmin基本不变。(4)在不同死端孔隙中,溶质粒子等待时间的概率密度函数具有相似的形式,即在前期随时间按幂律衰减,后期逐渐转变为指数衰减。两种衰减形式转变时刻与死端孔隙长度的平方成正比,与溶质扩散系数成反比。众所周知,等待时间随时间按幂律衰减表明溶质运移过程为反常运移过程,等待时间随时间按指数衰减表明溶质运移过程为菲克过程。因此,死端孔隙长度对溶质反常运移有很大影响。

关键词： 超临界二氧化碳   咸水层封存   孔隙尺度   数值模拟   两相流   毛细管数   接触角  

摘要： 以咸水层封存二氧化碳(CO_2)为研究背景,利用流体体积函数(Volume of Fluid,VOF)多相流模型建立孔隙尺度多孔介质计算模型,研究了超临界二氧化碳(Sc-CO_2)注入到含水多孔介质中的两相运移规律。对比分析了毛细管数、地质封存压力、Sc-CO_2注射温度、两相表面张力系数、接触角等因素对两相运移速率以及驱替效率的影响,同时将不同毛细管数下的驱替效率与实验数据进行了对比。研究结果表明,随着毛细管数的增大,驱替效率先减小然后趋于稳定,数值模拟与实验数据吻合良好;在同一孔隙率下,在壁面表现为亲水性时,壁面润湿性越好,驱替速率越快,但驱替效率有所下降;同时毛细管数越小、地质封存压力越低、注射温度越高、张力系数越小驱替速率越快,且驱替效率越高。