多孔介质中非等温互溶流体驱替过程的格子Boltzmann研究

利用格子Boltzmann方法,在孔隙尺度上对多孔介质内非等温混溶驱替过程进行数值研究,定量分析温度粘性膨胀系数（β_T）以及路易斯数（Le）对界面不稳定性和驱替效率的影响.结果表明：随着β_T的增大,界面不稳定性增强,驱替效率降低.当β_T>0时,随着Le的增大,界面不稳定性减弱,驱替流体与被驱替流体之间的界面趋于平缓,指尖残余率减小,驱替效率增大.当β_T<0时,Le对于驱替效率的影响相反.     

多孔介质中含溶解反应的互溶驱替过程格子Boltzmann研究

采用格子Boltzmann方法研究了孔隙尺度下多孔介质内含流固溶解反应的互溶驱替过程,重点研究了被驱替流体与驱替流体黏性差异较大的情况下,溶解反应引起的多孔介质内部结构变化对驱替过程的影响;定量分析了不同达姆科勒数及佩克莱数下多孔介质孔隙率和驱替过程驱替效率随时间的演变.研究结果表明:达姆科勒数较大时,溶解反应的发生会在多孔介质内部生成虫洞,导致一部分被驱替流体不能被波及,驱替流体沿虫洞离开多孔介质,造成驱替效率的减少.在此基础上,随着达姆科勒数的增大,孔隙率变化越大,生成的虫洞越宽,最终驱替效率变大,但仍小于无溶解反应时的驱替效率;随着佩克莱数的增大,指进增长速度越快,孔隙率变化越小,驱替效率越小.     

基于VOF方法的水驱油藏孔隙尺度剩余油分布状态研究

准确认识多孔介质内水驱剩余油的微观赋存状态对于改善高含水油田水驱开发效果,提高水驱采收率具有重要意义。在对Volume of Fluid (VOF)方法验证的基础上,充分利用其追踪两相界面动态变化、再现微观渗流物理过程的优势,开展特高含水期砂岩油藏物性条件、驱替方式对剩余油微观赋存特征和采收率的影响研究。通过分析典型孔隙结构的微观渗流特征和剩余油受力情况,揭示不同类型微观剩余油的动用机制和规律：水湿条件下驱替速度的增大和驱替方向的改变会使得微观剩余油分布较为分散且采收率得到不同程度的提高;油湿及高黏度比条件下毛管阻力、黏滞力较大,剩余油多以簇状和多孔状聚集分布,采出程度相对较低。     

基于相场方法的孔隙尺度油水两相流体流动模拟

基于真实岩心颗粒粒径分布，利用过程法构建疏松砂岩油藏的三维孔隙结构模型，利用相场方法建立两相流体流动数学模型并利用有限元方法进行求解，研究驱替速度、流体性质、润湿性对剩余油分布以及采出程度的影响.结果表明：驱替速度的增大和油水粘度比的减小会导致较大的毛管数，进而有利于采出程度的提高；就润湿性而言，水湿条件下毛管力是水驱油的动力，而在油湿条件下是阻力，因此水湿岩心采出程度更高.同时，从孔隙尺度对油水渗流机理及剩余油分布机理进行揭示，结果表明：由于多孔介质的复杂孔隙结构，流体在流经不同孔隙时呈现不同的流动特征，进而对油水两相流整体的压力分布、流速分布造成重要影响.     

含内热源多孔介质–自由流通道内非达西强制对流

基于多孔介质局部非热平衡模型,对考虑内热源条件下的多孔介质–自由流耦合通道内非达西对流换热特性进行研究。多孔介质区内流体运动方程采用Darcy-Brinkman-Forchheimer模型,利用有限差分法获得通道内各区域流体运动速度、流固相温度分布及努塞尔数,并进一步分析了相关参数对流体流动传热的影响。结果表明：在文中研究参数下,惯性参数F_f对通道内各区域流体速度及温度分布的影响仅在达西数Da大于10^-3时需要考虑;增加F_f或降低固相内热源W_s绝对值会使流固两相温差减小,且改变固相内热源换热方向会使多孔介质区内流固相发生温度分岔现象;固相内热源对通道换热效果影响较大且更为复杂,不同惯性参数F_f下,考虑W_s时可能会使Nu出现奇异点。     

扩散序谱(DOSY)实验测定缓冲体系中蛋白质表观分子量

液体核磁共振扩散序谱（DOSY）可以通过测定溶质分子的自扩散系数（D_t）来研究该分子在溶液中的表观分子量（M）.D_t与测试体系和分子本身性质相关，蛋白质体系较为复杂，从而增加了蛋白质自扩散系数（D_t-protein）测定的难度.本文以3-（三甲基硅基）丙磺酸钠（DSS）为内标，以蛋白质分子与DSS自扩散系数的比值（D_r）来表征蛋白质分子在溶液中的表观分子量（M_protein），该方法降低了缓冲体系对M_protein的影响，使得M_protein主要由分子本身的性质决定.在此基础上，测定了不同分子量蛋白质分子相对于DSS的D_r，拟合得到了D_r与M_protein的相关关系：lgM_protein =-2.6488 lgD_r-0.7863，相关系数（R²）为0.997.最后测定了通过大肠杆菌表达纯化得到的SARS冠状病毒主蛋白酶C端结构域（M^pro-C）分子相对于DSS的D_r，并计算出与文献结果一致的M_protein，进一步验证了拟合公式的准确性和实用性.     

多孔介质内混溶流体间粘性指进现象的格子Boltzmann模型

提出一个模拟多孔介质内混溶流体间粘性指进现象的格子Boltzmann模型.采用双分布函数分别求解压力场和浓度场.在浓度场平衡态分布函数中引入与浓度扩散相关的参数.通过调节参数,使粒子碰撞过程中的松弛时间保持恒定.模拟了粘度相同的流体间的混相驱替问题.不同网格下的模拟结果均与解析解吻合良好,验证了模型的可行性.进一步研究粘度比和贝克莱数(Pe)对粘性指进现象的影响.结果显示,增大粘度比会促进"手指"的增长.当粘度比不变时,存在Pe的临界值.当Pe超过临界值时,"手指"前缘会出现分裂现象.对横向平均浓度场的研究显示,混合区域的长度随时间的变化分为两个阶段,它首先随着t^1/2成线性增长,然后随着t成线性增长.     

多孔介质内泡沫渗流过程入口效应的数值研究

本文利用随机泡沫数目守恒模型对均质多孔介质内泡沫液渗流过程的入口效应进行了数值研究。随机泡沫数目守恒模型仅利用两个参数:泡沫生成速率Kg与最大泡沫个数n_∞,来描述多孔介质内非稳态驱替过程中泡沫的生成和发展。采用IMPES方法对二维介质内泡沫两相渗流过程的控制方程组进行求解,基于驱替过程中液相压力、液相饱和度及气泡密度分布等参数的分析,对泡沫驱替过程的入口效应进行了详细研究。结果发现,增大泡沫生成速率Kg,泡沫渗流过程的入口效应减小;而增大参数n_∞,虽可提高驱替压差,但对多孔介质入口段各参数的分布规律影响不大。     

CO＿2地质封存中重力对气水两相流动特性影响研究

CO₂地质封存技术是缓解温室气体排放的一种有效手段,亚岩心尺度CO₂-水在多孔介质内的流动实验对阐明CO₂地质封存中的残余捕集过程的封存潜力预测具有重要意义。本文利用高压微流控芯片在毛细管数为2.37×10^-5条件下进行了CO₂驱替纯水的实验,以研究重力对气水两相流动的影响。实时采集到了两相在非均质孔隙内的分布形态和流动特点,进一步的对CO₂饱和度进行了量化和评估了CO₂的封存容量,阐明了重力对CO₂残余捕集过程的影响。     

多孔介质BISQ模型中的慢纵波

着重研究了多孔介质BISQ模型中慢纵波的基本特性.给出了BISQ模型下慢纵波速度 和衰减的低频近似公式.与Biot理论对比，BISQ模型中慢纵波的衰减随频率降低急剧增大， 且随喷射流长度的减小而增加；相速度随喷射流长度的减小而增加，其低频极限值不是零； 孔隙流体位移与固相骨架位移之比的幅值随喷射流长度的增加而减小，其相位特点与Biot模 型预测的不同；在流体与孔隙介质的边界上可产生更大的渗流.为对比，同时也给出快纵波 的行为.依据BISQ模型可推断：非黏滞流体饱和孔隙介质中不存在喷射流机理；BISQ模型中 关键词： 多孔介质 喷射流 慢纵波 动力协调     

An analysis of the radial flow in the heterogeneous porous media based on fractal and constructal tree networks

In this paper, an analysis of the radial flow in the heterogeneous porous media based on fractal and constructal tree networks is presented. A dual-domain model is applied to simulate the heterogeneous porous media embedded with a constructal tree network based on the fractal distribution of pore space and tortuosity nature of flow paths. The analytical expressions for seepage velocity, pressure drop, local and global permeability of the network and binary system are derived, and the transport properties for the optimal branching structure are discussed. Notable is that the global permeability (K_n) of the network and the volume fraction (f_n) occupied by the network exhibit linear scaling law with the fractal dimension (D_p) of channel diameter bylogK_n∼0.46D_p and logf_n∼1.03D_p, respectively. Our analytical results are in good agreement with the available numerical results for steady-state soil vapor extraction and indicate that the fractal dimension for pore space has significant effect on the permeable properties of the media. The proposed dual-domain model may capture the characteristics of heterogeneous porous media and help understanding the transport mechanisms of the radial flow in the media.     

Effects of Monomer Size Distribution on the fractal dimensionality of diffusion-limited aggregates

Computer simulations of diffusion-limited aggregation (DLA) for monomers to investigate the effects of size and of lognormal distribution on the fractal dimensionality of the aggregates were conducted on a two-dimensional lattice. The results show the DLA clusters posses multifractal characteristics. For clusters consisting of monodisperse monomers, the bifurcation point on the graph of the pair correlation function (PCF) for each cluster is located right at the monomers size under investigation The textural dimension (D_f1) has a stable value of about 1.65, whereas the structural dimension (D_f2) decreased with increase in monomer size. For the cases with monomers in log-normal distributions, the textural dimension is around 1.67; however, the structural dimension decreases with increasing polydispersity of monomer size.     

Two—phase flow in correlated pore—throat random porous media

We have constructed a porous media model in which there are percolation clusters with varying percolation probability P and correlated site-bonds. Taking into account both the pore and the throat geometry, the viscous fingering (VF) in porous media has been investigated by using the standard over-relaxed Gauss-Seidel scheme. The simulation results show that the VF structure varies with the correlation parameter ε, the viscosity ratio M and the percolation probability P. The smaller the correlation parameter ε, the greater the deviation of the normalized size distribution of the invaded throat N_inv(r) from the truncated Rayleigh distribution. For a larger viscosity ratio M, the VF pattern looks like a diffusion-limited-aggregation structure in percolation clusters. The fractal dimension D increases with the increase of the percolation probability P and the correlation parameter ε. The velocity distribution f(α) of VF in percolation clusters is of a parabola-like curve. The tail of the distribution (large α) is longer for a larger correlation parameter ε. For a smaller ε, the distribution is very sharp. The sweep efficiency E decreases along with the decrease of the correlation parameter ε and the increase of the network size L_nz. E has a minimum as L_nz increases up to the maximum no matter what the values of P, M and ε. The E～ L_nz curve has a frozen zone and an active zone. The geometry and the topology of the porous media have strong effects on the displacement processes and the structure of VF.     

Computer simulations of two-phase flow with surface tension in the percolation cluster

Percolation clusters with varying occupation probability were constructed. Viscous fingering (VF) in the percolation cluster, based on the assumption that throat radii are Rayleigh distributed, is investigated by means of a successive over-relaxation technique. The fractal dimension and the sweep efficiency of VF in the percolation cluster when surface tension is considered are larger than when surface tension is neglected. The fractal dimension of VF will increase as the percolation probability increases or the viscous ratio decreases. VF's fractal dimension of porous media in the limit viscous ratio → ∞ is found to be identical with the DLA. The topology and the geometry of the porous medium have a strong effect on the displacement processes and the structure of the VF.     

On Different Approaches to Estimate the Mass Fractal Dimension of Coal Aggregates

Several methods to measure the structures of coal aggregates are compared. Loose and compact coal aggregates were generated through flocculation of ultrafine coal particles (mean volume diameter of 12 μm) under specific shearing conditions. Aggregate structure in terms of mass fractal dimension, D_f, was determined using various methods; namely 2D and 3D image analysis, interpretation of intensity patterns from small angle light scattering, changes in aggregation state through light obscuration, and settling behavior. In this study, the measured values of D_f ranged from 1.84–2.19 for coal aggregates with more open structures, and around 2.27–2.66 for the compact ones. All of these approaches could distinguish structural differences between aggregates, albeit with variation in D_f values estimated by the different techniques. The discrepancy in the absolute values for fractal dimension is due to the different physical properties measured by each approach, depending on the assumptions used to infer D_f from measurable parameters. In addition, image analysis and settling techniques are based on the examination of individual aggregates, such that a large number of data points are required to yield statistically representative estimations. Light scattering and obscuration measure the aggregates collectively to give average D_f values of the particulate systems; consequently ignoring any structural variation between the aggregates, and leaving possible small contaminations undetected (e.g. by dust particles or air bubbles). Appropriate utilization of a particular method is thus largely determined by system properties and required data quality.     

自仿射Sierpinski地毯中的粘滞指进

构造了自仿射Sierpinski地毯.在自仿射Sierpinski地毯中,认为喉管半径服从截断瑞利分布,采用逐次超松弛技术,模拟了自仿射Sierpinski地毯中的粘滞指进.计算了粘滞指进分维,结果表明,当粘滞比M→∞时,指进图样与在自仿射Sierpinski地毯中DLA模拟的结果类似：当M=1时,驱替流体在长标度范围内具有紧凑的结构,且具有稳定的位移. 关键词：     

Pore size distribution in porous glass: fractal dimension obtained by calorimetry

By differential Scanning Calorimetry (DSC), at low heating rate and using a technique of fractionation, we have measured the equilibrium DSC signal (heat flow) J _q ⁰ of two families of porous glass saturated with water. The shape of the DSC peak obtained by these techniques is dependent on the sizes distribution of the pores. For porous glass with large pore size distribution, obtained by sol-gel technology, we show that in the domain of ice melting, the heat flow J_q is related to the melting temperature depression of the solvent, ΔT _m , by the scaling law: J _q ⁰∼ΔT _m ^{- (1 + D)}. We suggest that the exponent D is of the order of the fractal dimension of the backbone of the pore network and we discuss the influence of the variation of the melting enthalpy with the temperature on the value of this exponent. Similar D values were obtained from small angle neutron scattering and electronic energy transfer measurements on similar porous glass. The proposed scaling law is explained if one assumes that the pore size distribution is self similar. In porous glass obtained from mesomorphic copolymers, the pore size distribution is very sharp and therefore this law is not observed. One concludes that DSC, at low heating rate ( q? 2^°C/min) is the most rapid and less expensive method for determining the pore distribution and the fractal exponent of a porous material. Received 23 July 1999 and Received in final form 16 February 2001     

FRACTAL VISCOUS FINGERING AND ITS SCALING STRUCTURE IN RANDOM SIERPINSKI CARPET

Viscous fingering (VF) in random Sierpinski carpet is investigated by means of successive over-relaxation technique and under the assumption that bond radii are of Rayleigh distribution. In the random Sierpinski network, the VF pattern of porous media in the limit M→∞ (M is the viscosity ratio and equals to η₂/η₁ where η₁ and η₂ are the viscosities of the injected and displaced fluids, respectively) is found to be similar to the diffusion-limited aggregation (DLA) pattern. The interior of the cluster of the displacing fluid is compact on long length scales when M=1, and the pores in the interior of the cluster have been completely swept by the displacing fluid. For finite values of M such as M≥10, the pores in the interior of the cluster have been only partly swept by the displacing fluid on short length scales. But for values of M in 1f(α) sites have velocites scaling as L^-α; and the scaling function f(α) is measured and its variation with M is found.     

Morphology Development in Multi‐Component Polymer Blends: I Composition Effect on Phase Morphology in PP/PET Polymer Blends

The relationship of the phase morphology of polypropylene/polyethylene‐terephthalate (PP/PET) blends and their corresponding compatibilized blends with composition was investigated using digital image analysis. A diameter, d _g , was defined and calculated to discuss the phase morphology of this polymer blend system. A figure‐estimation method was introduced to determine the width of the distribution of d _g . Based on the method, it is proven that the distribution of d _g obeys a log‐normal distribution and consequently, the distribution width, σ was calculated. Further, a fractal dimension, D _f , was introduced to describe the distribution of main sizes of the particles of the dispersed phase. The results showed that, while d _g increased with the concentration of the dispersed phase, σ and D _f show different dependence relations on composition;σ increases monotonously but D _f shows a maximum at a PET content of 30%, indicating that, even though the whole size distribution is much broader, the distribution of the main body of size becomes more uniform when the content of PET is less than 30%.     

Dynamic cracking resistance of metallic materials during rapid propagation of a self-similar crack

A model is proposed to determine the dynamic cracking resistance K _ID of metals and alloys for the case of a rapidly moving fractal or self-affine crack. The values of this characteristic correlate with the fractal dimension D _f of the future contour of a crack surface profile. K _ID is lower or higher than K _IC depending on the fractal dimension.