含甲烷水合物的石英砂渗透率实验和分形模型对比研究

含水合物的多孔介质渗透率是影响水合物开采的关键参数,多孔介质渗透率与水合物的饱和度密切相关。定量研究多孔介质渗透率随水合物饱和度的变化,对自然界中天然气水合物藏内渗流场的研究具有重要的理论价值。本文以平均粒径为139.612 μm的石英砂为多孔介质,采用稳态注水法测量在不同甲烷水合物饱和度（0 ~ 28.56%）下的石英砂渗透率,将实验数据与两种不同水合物赋存形式（颗粒包裹、孔隙填充）下的石英砂渗透率二维分形模型进行了对比。结果表明,石英砂渗透率比K_r随甲烷水合物饱和度S_h的增大呈现指数减小的趋势。当水合物饱和度低于11.83%时,渗透率比下降缓慢。而当水合物饱和度高于11.83%时,渗透率比下降迅速;当饱和度指数n = 12时,渗透率分形模型与实验数据吻合良好。通过分形模型与实验数据对比,发现当水合物饱和度低于11.83%时,甲烷水合物的赋存形式为颗粒包裹型。在11.83% ~ 28.56%水合物饱和度范围内,甲烷水合物的赋存形式为孔隙填充型。本研究成果量化了石英砂渗透率与甲烷水合物饱和度的关系,确定了含甲烷水合物的石英砂的渗透率分形模型的参数取值。     

基于电-力-声多物理场耦合数值模型的含水合物多孔介质声速和衰减特性研究

在物理模拟实验中对水合物微观赋存模式和饱和度进行准确控制和评价尚存在技术困难,仅依赖实验技术研究含水合物沉积物声学特性、建立储层参数解释模型存在局限性。采用基于有限元的数字岩石物理技术,针对悬浮、接触和胶结三种典型的水合物微观赋存模式分别建立多孔介质的三维电-力-声多物理场耦合模型,考察了微观赋存模式和水合物饱和度对多孔介质声速和衰减的影响规律,对比了声速数值模拟与理论模型计算结果,建立了声波衰减参数与水合物饱和度之间的关系式。研究结果表明：（1）对于三种水合物赋存模式,由于水合物相比孔隙水具有更高的弹性模量,多孔介质的声速随着水合物饱和度的增大而增大;水合物的存在导致声波在传播过程中遇到更多不连续的声阻抗界面,声衰减随着水合物饱和度的增大而近似线性增大;（2）悬浮和接触赋存模式条件下,水合物饱和度对多孔介质的声速和衰减影响规律基本一致;对于相同的水合物饱和度,胶结模式条件下含水合物多孔介质具有更高的声速和更小的声衰减;（3）通过合理选择参数值,利用权重方程与Lee改进的Biot-Gassmann Theory（BGTL）模型估算的含悬浮和接触模式水合物多孔介质的声速较为准确;通过等效介质理论模型C计算的含胶结模式水合物多孔介质的声速更为准确。研究结果可为获取复杂条件下含水合物沉积物的声学特性提供数值建模方法,为基于声波测井数据的水合物储层精细评价提供理论支撑。     

多孔介质中甲烷水合物的生成特性研究进展

天然气水合物是一种清洁高效的能源,常常在自然界中的海底沉积物多孔介质孔隙中生成,同时水合物在工业上还能与多孔介质材料一起作为储存及分离气体的一种方式,因此开采利用水合物以及发挥水合物工业技术的前提都跟多孔介质有莫大的关系,对多孔介质中天然气水合物生成特性的研究进行总结与分析具有非常重要的意义。本文总结分析了国内外关于不同类型多孔介质中甲烷水合物的生成过程及特性的研究文献,将多孔介质根据其孔径大小进行划分。结果显示,在微孔介质中,甲烷水合物的生成侧重于气体的存储及运输方面;在介孔介质中,甲烷水合物的生成动力学受孔径影响较大;在大孔的沉积物中,甲烷水合物的生成及分布的机理性研究仍比较缺乏。因此,需要进一步的研究来丰富甲烷水合物在多孔介质中的生成动力学理论,本文将在文献调研的基础上为今后的研究方向提出一些展望和思路。     

多孔介质中甲烷水合物的微观分布对电阻率的影响

电阻率测井技术是天然气水合物储层资源量估算的重要手段,明确储层电阻率与水合物在沉积物孔隙中微观分布的关系对准确估算水合物饱和度有重要意义。使用自主研发的天然气水合物计算机断层扫描（CT）技术-电阻率测量装置,开展多孔介质中甲烷水合物生成实验,同时进行CT观测与电阻率测试,得到了反应过程中扫描图像和电阻率数据,讨论了孔隙中水合物在不同微观分布模式下电阻率与水合物饱和度的关系。结果表明：水合物饱和度低于10.50%时,水合物以接触分布模式为主,与水合物对孔隙的填充效应相比,排盐效应更加显著,电阻率随水合物饱和度增大略有下降;水合物饱和度为10.50% ~ 22.34%时,水合物为悬浮与接触共存分布模式,悬浮分布模式下的水合物对孔隙水连通截面的阻塞作用较大,电阻率随水合物饱和度的增大急剧升高;水合物饱和度为22.34% ~ 27.50%时,水合物仍然为悬浮与接触共存分布模式,在大部分孔隙已被堵塞的情况下水合物的填充作用对电阻率影响较小,电阻率随水合物饱和度增大而升高的趋势逐渐变缓。可见,电阻率响应特性在不同的水合物饱和度范围下明显不同,与水合物的微观分布有关。     

天然气水合物储层渗透率研究进展

天然气水合物是清洁、高效、储量巨大的未来最具潜力资源之一,而储层的渗透率是影响水合物资源开采时产气效率的重要参数。目前,国内外对天然气水合物储层的渗透率进行了大量的研究并取得一定进展,本文从数值模拟、储层现场探井、实验研究等方面全面回溯,分析总结了孔隙度、饱和度、应力应变情况等对储层渗透率的影响,讨论目前储层渗透率研究存在人工合成水合物沉积物与自然储层存在差异,不同多孔介质形成的水合物沉积物应力敏感性不同,多相渗流研究不够充分,储层渗透率改造研究不足等问题,对未来的研究方向提出了展望。     

甲烷水合物及其沉积层导热特性的研究现状

总结了近年来国内外甲烷水合物及其沉积层导热特性的研究现状,从实验测试和模拟研究两方面分析了甲烷水合物及其沉积层的导热影响因素、相关规律和导热机理。研究结果表明：纯质水合物的导热性能与外界温度、压力、客体分子数量、多孔介质、笼形结构等因素有关,其值大小主要由主体水分子形成的笼形结构决定;而水合物沉积层导热系数的大小主要依赖于实验样品的组成成分、初始水饱和度、各组分分布,与温度、压力以及垂直有效应力关系不大。最后,指出了现存研究方法中一些值得改进的地方,并对今后的研究工作进行了展望。     

多孔介质内甲烷水合物生成动力学研究

天然气水合物作为一种潜在的新型潜在替代能源,分布广泛、储量巨大。水合物生成动力学特性对于了解自然界天然气水合物形成规律与水合物相关技术应用具有重大意义。针对甲烷水合物生成动力学特性进行研究,设计开发了甲烷水合物生成测试系统,研究了不同温度和不同孔隙粒径的多孔介质对于水合物生成动力学的影响。结果表明,随着温度的降低以及多孔介质孔隙粒径的减小,甲烷水合物生成过程中的气体消耗速率越快,生成结束时消耗的气体的量越多。同时,利用水合物生成多步骤机制模型,获得了甲烷水合物生成过程各中间产物的变化规律。     

基于TDR双参数的沉积物中水合物饱和度测量

针对沉积物中水合物饱和度的测量问题,立足于时域反射技术（TDR）能够同时获得含水合物沉积物表观介电常数和电导率的优点,提出了一种基于介电常数/电导率双参数的水合物饱和度评价新方法。在分析介电常数和电导率测量原理的基础上,设计了四氢呋喃（THF）水合物模拟实验与参数测量系统以及实验方案;基于TDR测量响应分析了孔隙水电导率对水饱和的以及含水合物的模拟沉积物介电常数和电导率的影响;利用TDR获取的不同水合物饱和度条件下的介电常数测量数据对比分析了经典介电常数模型的性能,并以Lichteneker-Rother（LR）模型为原型建立了水合物饱和度与介电常数之间的关系模型,以阿尔奇公式为原型建立了基于TDR测量电导率的水合物饱和度计算模型。研究结果表明：采用LR模型和Maxwell-DeLoor模型能够较准确地描述TDR测量表观介电常数与水合物饱和度之间的关系;基于表观介电常数/电导率双参数的水合物饱和度联合评价方法为提高评价结果的准确度和可靠性提供了新途径。将来需要结合水合物储层实际特征开展模拟实验并进一步完善饱和度计算模型,将基于TDR双参数的饱和度联合评价方法推广应用到含天然气水合物沉积物。     

质子交换膜燃料电池的气体扩散层渗透率分形模型

本文提出了质子交换膜燃料电池（PEMFCs）气体扩散层（GDL）分形渗透率模型。这个模型是根据扩散层真实微观结构中的两个分形维数建立的。其中一个与毛细管流通道大小有关,另一个与通道迁曲度的描述有关。此外,气体分子的影响可以通过Adzumi方程计算。渗透率分形模型是多孔介质迁曲度分形维数、孔隙面积分形维数、孔径以及有效孔隙度的函数,模型中没有任何经验常数,可以用压汞法测量扩散层的微观结构。根据扫描电子显微镜图象,可用盒式维数法确定两个分形维数。为了检验模型的正确性,把该模型渗透率的预测数据与Toray提供的实验数据进行对比,发现该模型的渗透率预测与实验数据一致,证实了气体扩散层的分形渗透率模型的正确性。     

泡沫铜多孔材料的液体渗透特性研究

目前对波形板汽−水分离装置的研究多集中于结构参数和入口参数的分析,而从其材料特性角度探讨提高分离效果的研究鲜有报道。以去离子水为介质,探讨了泡沫铜多孔材料孔隙率、有效半径以及浸润性对其渗透特性的影响。将修正后的弯毛细管模型与达西定律相结合,建立了泡沫铜渗透率、孔隙率以及有效半径之间的定性关系,即随着泡沫铜孔隙率和有效半径的增加,其渗透性能逐渐增强,这一结论与实验结果符合良好。此外,通过对泡沫铜材料进行浸润改性,分析了浸润性对其渗透性能的影响,结果发现,随着表面能的提高,泡沫铜的渗透性能逐渐增强。     

Application of lattice Boltzmann method to a micro-scale flow simulation in the porous electrode of a PEM fuel cell

The electrode of a PEM fuel cell is a porous medium generally made of carbon cloth or paper. Such a porous electrode has been widely modeled as a homogeneous porous medium with a constant permeability in the literature of PEM fuel cell. In fact, most of gas diffusion media are not homogeneous having non-isotropic permeability. In case of carbon cloth, the porous structure consists of carbon fiber tows, the bundles of carbon fiber, and void spaces among tows. The combinational effect of the void space and tow permeability results in the effective permeability of the porous electrode. In this work, the lattice Boltzmann method is applied to the simulation of the flow in the electrode of a PEM fuel cell. The electrode is modeled as void space and porous region which has certain permeability and the Stokes and Brinkman equations are solved in the flow field using the lattice Boltzmann model. The effective permeability of the porous medium is calculated and compared to an analytical calculation showing a good agreement. It has been shown that the permeability of porous medium is strongly dependant on the fiber tow orientation in three-dimensional simulations. The lattice Boltzmann method is an efficient and effective numerical scheme to analyze the flow in a complicated geometry such as the porous medium.     

Multi-phase micro-scale flow simulation in the electrodes of a PEM fuel cell by lattice Boltzmann method

The gas diffusion layer of a polymer electrolyte membrane (PEM) fuel cell is a porous medium generally made of carbon cloth or paper. The gas diffusion layer has been modeled conventionally as a homogeneous porous medium with a constant permeability in the literature of PEM fuel cell. However, in fact, the permeability of such fibrous porous medium is strongly affected by the fiber orientation having non-isotropic permeability. In this work, the lattice Boltzmann (LB) method is applied to the multi-phase flow phenomenon in the inhomogeneous gas diffusion layer of a PEM fuel cell. The inhomogeneous porous structure of the carbon cloth and carbon paper has been modeled as void space and porous area using Stokes/Brinkman formulation and void space and impermeable fiber distributions obtained from various microscopic images. The permeability of the porous medium is calculated and compared to the experimental measurements in literature showing a good agreement. Simulation results for various fiber distributions indicate that the permeability of the medium is strongly influenced by the effect of fiber orientation. Present lattice Boltzmann flow models are applied to the multi-phase flow simulations by incorporating multi-component LB model with inter-particle interaction forces. The model successfully simulates the complicated unsteady behaviors of liquid droplet motion in the porous medium providing a useful tool to investigate the mechanism of liquid water accumulation/removal in a gas diffusion layer of a PEM fuel cell.     

A 2D analysis of stability of bioconvection in a fluid saturated porous medium — estimation of the critical permeability value

The main purpose of this paper is to perform a 2D stability analysis of bioconvection in a suspension of motile gyrotactic microorganisms in a fluid saturated porous medium and to obtain an analytical expression for the critical permeability of the porous medium. Recent numerical investigation by Kuznetsov and Jiang [1] suggests that permeability is a very important parameter for bioconvection in porous media. Their numerical results indicate that for small permeability bioconvection is stable (the microorganisms swim in the upward direction), while for large permeability it is unstable (variations of density are enhanced and macroscopic fluid circulation is induced). In the present investigation, a simple but elegant criterion of stability of the bioconvection is obtained. This criterion gives the critical permeability of the porous medium through the cell eccentricity, average swimming velocity, fluid viscosity, and other relevant parameters.     

多孔介质金属泡沫的渗透率研究

基于真实金属泡沫内部结构的特征,建立了与真实多孔介质形态相近的十四面体模型及相对应的代表性单元体。采用分形理论得到真实多孔介质的孔隙率和渗透率,由十四面体的结构特征得到了相应的孔隙率和渗透率,根据流动一致性原理,确定十四面体模型的结构参数。通过十四面体的结构特征得到渗透率的表达式,并且将模型渗透率的计算结果和现有的实验数据进行对比,结果存在良好的一致性,说明渗透率表达式的有效性。     

Two-phase modeling of mass transfer characteristics of a direct methanol fuel cell

A two-dimensional, two-phase mass transport model has been developed for a direct methanol fuel cell (DMFC). The model is numerically solved with a in-house code and validated with published experimental data in the literature. In particular, gaseous and liquid phase velocities in the anode porous structure are obtained so that the liquid–gas counter convection effect can be investigated. The numerical results show that the mass transfer of methanol is predominated by the resistance in the anode porous structure, which is affected by physical properties of the porous medium, such as porosity, permeability, and contacting angle. It is further shown that cell performance can be improved by increasing the porosity and permeability, and decreasing the contacting angle of the porous medium for a given feed methanol concentration.     

A simple mathematical model for determining the equivalent permeability of fractured porous media

A simple mathematical model has been proposed so as to determine the equivalent permeability of fractured porous media. The model consists of square blocks placed in an array with vertical and horizontal fractures between the blocks. An analytical expression valid for all macroscopic flow directions has been derived for the equivalent permeability of the fractured porous media, assuming a horizontal flow through the blocks placed in a porous medium. The analytical expression agrees well with the existing equations and also with the microscopic numerical results carried out using a unit structure with periodic boundary conditions. The foregoing two-dimensional model has been extended to a three dimensional case in which the cubic rocks are arranged in a cubic array. The resulting three-dimensional analytical expression for the equivalent permeability is found to agree very well with both existing formula and microscopic numerical simulation.     

Couette flow through a porous medium of a high prandtl number fluid with temperature-dependent viscosity

The velocity and temperature profiles for an impulsively started Couette flow have been derived for a fluid with a high and strongly temperature-dependent viscosity when the flow takes place through a porous medium. The steady as well as the transient state flows are discussed and the influence of the medium permeability is assessed. In the steady state the presence of the porous medium causes higher maximum temperatures only for sufficiently low permeabilities and in all cases significantly lower velocity profiles. The flow development times tend to be higher for the low permeabilities only for high Nahme numbers. For the transient state it is noticeable that, in all cases, the velocity develops much more rapidly than the temperature and that the presence of the porous medium accelerates this tendency. Finally the medium permeability produces skew temperature profiles and higher temperature time gradients at all times.     

The flow of non-Newtonian fluid in an inclined channel through variable permeability

A non-Newtonian fluid's Poiseuille flow in a porous medium with variable inclination and permeability is investigated. Let us assume for the sake of simplification that permeability varies as a quadratic parabolic function form. The porous medium is used by the Brinkman methodology to control the flow. The equations for velocity distribution and mass flow that result from this are evaluated using different input values. This problem describes the effect of inclination, Jeffrey parameter, and variable permeability on the classical Poiseuille flow between parallel plates. This problem can also be treated as an extension of the work of Hamdan and Kamel for non-Newtonian fluid flow in an inclined channel. Also, the effects of these variables on the variation of mass flux with Jeffrey parameter λ₁ is analyzed through graphs, and the skin friction coefficient is analyzed through table values. It is observed that the maximum permeability of the porous medium affects both the mass flow rate and the velocity, which increase with rising λ₁ and decrease with rising H_a, respectively.     

Analysis of flow and heat transfer characteristics of porous heat-storage wall in greenhouse

The flow and heat transfer characteristics of porous heat-storage wall in greenhouse are studied by using the one-dimensional steady energy two-equation model for saturated porous medium. The results show that the heat exchange between the air and the solid matrix of the porous heat-storage wall depends upon the inlet air velocity, the porosity and the permeability of porous medium, and the thermal conductivity of the solid matrix. Because the incidence of solar radiation on the porous heat-storage wall is not uniform, the new composite porous solar wall with different porosity is proposed to reduce the disadvantageous effect.     

Investigations on the formation kinetics of semiclathrate hydrate of methane in an aqueous solution of tetra-n-butyl ammonium bromide and sodium dodecyl sulfate in porous media

Natural gas hydrate is considered to be an attractive sustainable energy resource for the world. Hydrate as a technology can be of immense importance for various industrial processes, such as multicomponent natural gas separation, gas storage and transportation, and carbon dioxide capture from flue gases and sequestration. A variety of hydrate additives, which includes promoters (thermodynamics and kinetics) and porous media, are being researched to improve the hydrate formation kinetics. However, studies involving the combinations of these are rare in the open literature. In this work, the formation kinetics of methane hydrate/semiclathrate hydrate using tetra-n-butyl ammonium bromide (TBAB) and sodium dodecyl sulfate (SDS) aqueous solutions at various concentrations in a porous medium containing silica sand at initial hydrate formation pressures (7.5 and 5.5 MPa) and temperatures (273.65 and 276.15 K) have been investigated. All the experiments were conducted using 75% water saturation. Various kinetics parameters, such as gas uptake, gas-to-hydrate conversion, and induction time, have been reported. It was found that the combination of TBAB+ SDS showed favorable hydrate formation kinetics in porous media than the TBAB system. This work provides information for further studies involving semiclathrate hydrate applications for various industrial processes.