采空区漏风及残煤瓦斯涌出对采场气流的影响

在国内外文献调研的基础上,通过渗流力学、岩石力学、采矿工程等多学科交叉,采用理论分析、数值计算结合的方法,建立了煤矿采空区瓦斯流动数学模型。采用CFD技术对漏风流分布及瓦斯浓度场进行研究,获得了整个工作面上风流的速度和瓦斯浓度的分布情况,得到了漏入采空区和采空区返回工作面风流较为集中的范围,采空区漏风量与工作面风量之间的变化关系,漏风流速沿采空区衰减规律,确定了整个工作面上瓦斯浓度最高的地方为回风隅角,采空区内瓦斯浓度最高的区域也是整个采场瓦斯浓度最高的区域。通过对不同送风量的模拟研究,描述了瓦斯浓度的运移趋势,得到了回风隅角瓦斯浓度与送风量的关系,并且发现近风侧瓦斯浓度梯度较大,增大送风量能降低近回风隅角瓦斯浓度,但是会增加遗煤自燃的风险。揭示了采空区漏风和残煤瓦斯涌出对采场气体流动的影响规律。     

综放开采围岩活动影响下瓦斯运移规律及其控制

 博士学位论文摘要　将综放开采矿山压力新理论、覆岩结构及其活动特征与瓦斯运移聚集规律有机结合, 通过现场矿压观测、适时瓦斯监测、相似材料模拟、电液伺服试验、理论分析及有限元解算等方法, 对含瓦斯厚煤层实施综放工艺时顶煤顶板活动规律及其内瓦斯运移聚集形态进行了系统深入的定性定量分析。从分析含瓦斯煤岩体微观结构特征入手, 探讨了矿压作用下煤岩与瓦斯的宏观运动形态, 从而得出综放开采富含瓦斯煤层的安全控制理论与实践技术。论文视含瓦斯煤岩体为不连续的多相介质材料, 并深入到材料的微观机制抽象出其微结构模型, 分析了其结构形态在煤岩体变形和瓦斯运移过程中的变化形式; 定量分析了煤岩体内孔隙裂隙分布及破碎后的块度变化的分形维数值, 认为采动后支承压力等因素使煤岩体破坏, 从而影响其分维值, 煤岩的极限粒度越小,瓦斯放散解析速度越快, 瓦斯运移潜能越大, 表征其吸附特性的V an derW aals 力下降越多; 综放开采使煤岩力学介质变化显著, 支承压力作用下煤体经历初期压缩、强烈压缩、滑移与离散、流动放出4 个过程, 顶煤介质状态可由宏观连续介质转化为碎裂介质和块裂介质。基于魏家地矿和阳泉五矿的综放面及巷道矿压观测、瓦斯监测和采空区束管监测, 认为本煤层开采时支承压力使煤体孔隙率、渗透系数和瓦斯压力发生很大变化, 且支承压力变化状态与瓦斯运移变化在时间和空间上有一致性; 首次将综放面前方煤体划分为无压瓦斯自由放散区、卸压瓦斯涌出活跃区、降压瓦斯涌出变化区、升压瓦斯涌出变化区及稳压瓦斯正常涌出区; 周期来压时, 瓦斯绝对涌出量和相对涌出量分别比来压前增加了44% 和54%; 统计国内11 个矿井生产面, 覆岩关键层初次来压时, 采场瓦斯绝对涌出量是来压前的1. 21～ 3. 44 倍, 相对涌出量是初垮前的1. 21～ 6. 75 倍; 关键层失稳经历不同, 使综放采空区后方垮矸碎胀系数具有明显的分区性, 据计算碎胀系数和瓦斯浓度、氧浓度变化可将其划分为自然堆积区、载荷影响区和压实稳定区, 自然堆积区和载荷影响区是分析瓦斯流态的关键区, 其内瓦斯的紊流和过渡流态可用非线性渗流方程表示。关键层初次破断失稳的动态冲击力, 一方面克服了风流压力和漏风阻力将富集于支架上方断裂煤壁及架后采空区瓦斯压向采场, 另一方面在碰撞作用点附近使采空区垮矸空隙率降低形成模拟墙, 从而使墙体前方瓦斯迅速挤向采场; 顶板来压前工作面支承压力达到最大, 使直接顶、顶煤体剪切破坏, 煤体屈服, 卸压范围扩大, 支承压力与瓦斯压力梯度联合作用使煤体瓦斯解析并向采场涌出; 关键层作为板的破断呈“O 2X”型特征, 在其交点处首先形成导气通道, 顶板来压期间关键层及其覆岩因变形特性不同而不协调垮落, 便将离层裂隙聚集的瓦斯通过“O 2X”破断裂隙挤入采场, 因此得出结论: 综放面瓦斯大量快速涌出是矿山压力的一种显现。综放开采使垮落带和规则移动带高度增加, 为瓦斯运移聚集提供了较大活动范围。关键层与其上覆或下伏岩层间不协调变形将形成覆岩离层裂隙和破断裂隙; 在煤层采厚2. 6～ 3. 4 倍高度以下破断裂隙较发育, 其上以离层裂隙为主, 随综放开采两类裂隙的时空发展有明显的三阶段特征, 即切眼和回采面附近覆岩采动裂隙发育, 采空区中部裂隙则被重新压实; 覆岩垮落及离层高度受关键层及其结构的影响而呈动态变化, 随顶板初次来压和周期来压, 离层高度呈跳跃性变化, 从而为裂隙带内瓦斯运移聚集的剧变提供了时空条件。综放开采前期是瓦斯运移及控制的关键时期, 覆岩采动裂隙带是经破断与离层裂隙贯通后在空间形成关键层下似椭圆抛物面内外边界所包围的椭抛带(EPZ) 分布, 椭抛带层面的切割为椭圆形裂隙发育区; 关键层破断后, 裂隙带宽在初采边界处相当于初次来压步距, 在综放面上方则变化在1～ 2 倍周期来压步距之间; 采动裂隙带的发生、发展基本受制于覆岩关键层层位及其所形成砌体梁结构的变形、破断和失稳形态; 当主关键层切割椭抛带时, 采动裂隙呈椭球台状, 层面展布的椭圆形裂隙区仍将存在。表征覆岩离层特征的位移曲线主要取决于关键层断裂块度的大小, 由此可计算出关键层初次失稳前后离层裂隙的当量面积及其空隙率和渗透系数。综放开采时瓦斯涌出特点决定了其在覆岩采动裂隙带内具有升浮和扩散两种运移方式。不均衡性瓦斯涌出带, 与周围环境气体存在密度差而升浮, 在浮力作用下沿破断裂隙上升过程中不断渗入周围气体, 使涌出源瓦斯与环境气体的密度差渐减至零, 瓦斯则会漂浮在离层裂隙发育区, 瓦斯升浮高度与本煤层及邻近层瓦斯含量及涌出强度成正比; 混入矿井空气中的瓦斯(CH4) 在其浓度梯度作用下会引起气体分子的普通扩散和压强扩散, 瓦斯扩散流方向与重力压强梯度反向, 即瓦斯具有向上扩散的趋势。从理论上解释了裂隙带是瓦斯运移及聚集带, 为覆岩裂隙带内钻孔抽放、巷道排放瓦斯技术提供了科学依据。首次提出煤样全应力应变过程中渗透系数是体积应变的双值函数, 体积缩小时为2 次多项式, 体积膨胀时为5 次多项式; 渗透系数在弹塑性段急增, 峰值后仍增大, 但梯度渐缓, 最大值发生在软化段或塑性流动段, 且与最小值相差上百倍; 主应力差增大时, 渗透系数变化范围增大, 反之则小。渗透系数是影响煤层瓦斯运移的最重要指标, 而支承压力则是影响渗透系数的主导因素; 支承压力作用下综放面前方不同部位煤体渗透系数变化范围相当大, 支架上方顶煤煤体及煤壁前方5m 内渗透系数最高, 塑性变形区内, 煤层渗透系数急剧降低, 到弹性变形区则接近原始值, 两极值相差可达数10 倍甚至数百倍。认为不论原始渗透系数怎样低的煤层, 采动影响下煤层卸压后, 其内瓦斯渗流速度大增, 瓦斯涌出量也随之剧增, 为瓦斯抽排提供便利条件, 由此提出“煤层与瓦斯共采”的新概念。有限元计算表明: 均衡推进的综放面, 采用短距离循环推进则可降低煤体中因渗流场结构变化而引起的瓦斯压力较大的波动, 一定程度上可减弱综放面前方煤体中瓦斯挤压和抛出煤体(动力异常) 的危险程度。提出以采场矿压监测为主的连续危险源非接触式法预测采场瓦斯大量涌出或涌出异常, 并成功预测了打通一矿工作面突出; 魏家地矿的应用实践表明, 三巷型布置较适宜富含瓦斯倾斜厚煤层的综放开采, 高抽巷应开掘在采动裂隙带内; 预采顶分层或开采解放层即可预释放大量瓦斯又可减缓综放开采的矿压显现程度;代替采空区井的采动区井(孔) 底处于裂隙带内能够充分抽排瓦斯, 淮北局的应用充分证明了该论点; 充分监测综放面顶板来压, 可有效防止综放采空区瓦斯爆炸; 充分利用矿压显现特点且有利于综放开采防治瓦斯、煤层自燃和煤尘等的煤体注水技术在魏家地矿取得了成功。     

煤体双重孔隙瓦斯双渗流模型及无因次解算

 考虑瓦斯在煤层中的解吸、放散与渗流,利用达西定律分别描述煤基质与裂隙内的瓦斯运移,以煤基质与裂隙之间的传质通量为桥梁,发展煤体双重孔隙瓦斯双渗流模型,推导无因次模型,并运用有限差分法进行编程解算。结果表明：瓦斯压力、含量在裂隙内的下降速度要远大于煤基质;基质空间内瓦斯压力及含量的分布具有非均匀性及非稳态性;增大裂隙渗透性或煤层瓦斯压力,或减小煤壁表面瓦斯压力,均能导致瓦斯涌出速度的增大;煤体游离瓦斯含量对瓦斯涌出速度影响较小。结合潘一矿煤层瓦斯参数,对比模拟结果和实测数据,验证了煤体双重孔隙瓦斯双渗流模型的正确性。     

综放开采J型通风采空区渗流场数值分析

以潞安矿区王庄煤矿5201综放面通风和瓦斯抽放为背景，分析J型和U型通风时采空区破碎岩体内的渗流和瓦斯浓度分布规律。J型通风采空区渗流为“一源两汇”，并有明显取近道的趋势，J型通风比U型通风更有利于瓦斯排放。为建立破碎岩体渗流数值模型，还兼顾分析采空区上方的支承压力，得出该处支承压力的分区特征。     

页岩中气体流动规律数值模拟研究

为了更好地对页岩气进行开发,亟需对其渗流规律进行深入研究。常规气藏的气体渗流理论没有考虑气体在页岩纳米尺度孔隙中流动规律的特殊性。本文考虑气体在页岩微观孔隙中流动时的Knudsen扩散和壁面滑脱效应,建立了针对页岩一维流动的非线性渗流数学模型及求解方法,研究了页岩中气体的渗流特征和规律。通过实验,对本文的数学模型和常规气藏数学模型的结果进行了验证,可以看出:本文的模型结果能更好地符合实验规律,而常规气藏模型计算的流动平衡时间和实验结果有50%左右的差别,因此考虑气体在页岩中的特殊渗流规律的影响,对于页岩气渗流规律的研究是必要的。     

U+L型通风系统采空区瓦斯运移规律数值模拟

文章采用FLUENT软件对高瓦斯煤层"U+L"型通风系统采空区气体渗流进行数值模拟,得出采空区瓦斯浓度分布状况及其运移规律。对于高瓦斯矿井"U+L"型通风系统瓦斯综合治理与利用具有显著的实践指导意义。     

破碎岩体的渗流理论及其应用研究

以国家自然基金重点项目“煤矿瓦斯灾害预防及煤层气开采中的应用基础研究”和教育部博士点基金“峰后应力状态下的岩石渗透性研究”为背景，通过理论、试验、数值模拟和结合潞安王庄煤矿采空区通风渗流问题的具体应用，进行了如下具有创新意义的工作：(1)自行研制开发出目前已获国家专利的破碎岩石压实渗透仪，并给出了一套进行破碎岩石压实水渗透测试的试验方法。利用开发的渗透仪和试验方法，系统测定并分析了砂岩、页岩和煤3种岩样在不同破碎程度下压实的渗透系数变化规律。(2)利用置信区间方法、随机过程理论深入分析了破碎岩体渗透性的随机性特性，并建立破碎岩体渗流的随机微分控制方程。(3)用控制体欧拉法建立了描述破碎岩体塑性大变形的基本方法，开展了试验背景下破碎岩体渗流场的数值模拟。推导出求解破碎岩体渗流方程的一、二阶摄动随机有限元列式。(4)计算并给出采空区破碎岩体的支承压力分布规律，建立了采空区破碎岩体渗流分析的数值模型，结合潞安王庄煤矿综放面采空区的破碎岩体通风渗流，系统分析并得出了与理论和实测相吻合的流场分布和近似瓦斯浓度的分布规律。     

考虑吸附膨胀应力影响的煤层瓦斯流-固耦合渗流数学模型及数值模拟

已有研究表明,煤吸附瓦斯后要发生相应的膨胀变形,而当这种膨胀变形受到一定的限制时,将会产生膨胀应力。为此,在考虑煤吸附瓦斯产生膨胀应力的前提下,根据煤体受力平衡条件建立考虑吸附膨胀应力的煤体有效应力表达式,并根据流–固耦合渗流理论的基本思想,建立煤层瓦斯流–固耦合数学物理模型。采用有限元和有限差分相结合的方法对模型进行离散化处理,并对其进行模拟求解。和未考虑吸附膨胀应力的流–固耦合模型的模拟结果进行比较分析,研究成果表明,未考虑吸附膨胀应力的流–固耦合模型计算的压力下降较该模型要小,两种模型的模拟结果存在一定的差别,从而表明在煤层瓦斯流–固耦合数学模型中应考虑吸附膨胀应力的影响。     

低渗透油藏CO2混相驱过程中考虑对流扩散的流固耦合模型研究

为研究低渗透油藏中对流扩散和流固耦合对CO2混相驱油的影响,推导低渗透油藏CO2混相驱过程中多相多组分渗流场、对流扩散场和应变场三场全耦合控制方程,并给出求解方法。模型考虑对流扩散、混相和岩石骨架的应力–应变对整个渗流过程的影响,并对边界条件进行约束。应用建立的考虑对流扩散的流固耦合模型进行实例分析,结果表明:所建模型综合考虑对流扩散与流固耦合现象对渗流的影响,体现出模型与计算方法的实用性,并提高CO2混相驱的模拟精度,本模型计算结果与常规数值模拟结果具有一定的可比性。同时,通过界面张力、产能、混相区域等参数定量描述对流扩散现象对CO2混相驱的影响。     

煤层瓦斯渗流与煤体变形的 耦合数学模型及数值解法

根据固流耦合作用的基本理论，考虑瓦斯渗流对煤体本构关系的影响，提出了煤和瓦斯耦合作用的数学模型，给出了用有限元求其数值解的方法及程序。在此基础上，对采动影响下考虑煤层变形影响时瓦斯在采空区的流动规律进行了数值模拟分析，为采空区瓦斯抽放提供了理论基础和科学依据     

预测压气法隧道施工中的漏气量的模型研究

在阐述压气法隧道施工原理的基础上,建立了水–气二相流模型来模拟压气法隧道施工过程中超压空气驱替隧道周围土体的孔隙水的水–气二相流过程。该模型同时考虑了土体中气相和水相的流动以及水相和气相相互作用的影响,比较符合水–气二相渗流过程的物理意义。通过工程算例,计算出不同施工进度下开挖面和隧道混凝土衬砌上的空气损失速率以及空气入流对隧道周围土体渗流场的影响,结果表明:在均质的地质条件下,开挖面上的漏气速率一般为常数,而衬砌上的漏气速率将随着隧道长度的增加而线性增大。最后,提出了利用水–气二相流模型进一步提高压气法隧道施工的漏气量的数值模拟精度,需要深入研究的两点意见。     

Characterization of two-phase flow in a single rock joint

An experimental study using a triaxial apparatus was used to analyze the two-phase flow patterns in jointed rock specimens. Rock specimens having a single natural fracture were tested for two-phase flow of water and air. Triaxial tests were conducted to characterize the two-phase flow through fractured granite specimens at low confining pressures. It was found that for a relatively smooth joint (JRC<6), bubble flow pattern occurred within the rock joint when the gas velocity is below 15 m/s. The average velocity of water usually varied between 0.1 and 0.5 m/s for bubble flow patterns. In this velocity range, air bubbles were able to form along the joint walls or to be randomly displaced within the water phase. When the gas velocity inside the rock joint exceeded 22 m/s, the flow patterns took annular form for non-zero capillary pressures (i.e., injected gas pressure is not equal to injected water pressure). At elevated (>0.25 MPa) gas injection pressures, the gas occupied the main part of the fracture and the liquid was able to flow as an unstable film forming an annular flow along the joint. When the annular flow developed, the mixture flow pattern was independent of the air flow velocity. This was due to the fact that once the injected air velocity reached a critical value (i.e., 20 m/s), water velocity inside the joint was negligible for a given confining pressure and injected water pressure. Further increase in inlet air pressures developed a single-phase air flow with no water flow.     

水位下降过程中气相对土坡稳定性的影响

基于饱和-非饱和水-气二相流模型,对土坡在水位下降过程中的孔隙水压力、孔隙气压力、毛细压力及水相饱和度的变化过程进行模拟;基于模拟得到的不同时刻的渗流状态,采用Bishop简化方法推导出考虑孔隙气压力及坡外水压力的土坡安全系数计算公式,分析了水位下降过程中不同渗透性土坡安全系数的变化规律,及基质吸力和气相的存在对边坡稳定性的影响。结果表明:水位下降过程中土坡的安全系数逐渐减小,水位下降后土坡的安全系数略有增大;考虑气相影响后土坡的安全系数降低,考虑基质吸力影响后土坡的安全系数增大;岸坡土体的渗透率越大,水位下降过程中安全系数较大,而水位下降后安全系数较小。     

砂土渗流过程的细观数值模拟

基于散体介质理论,利用PFC2D内置FISH语言定义的流固之间的作用力方程和压力梯度方程,求解不可压缩流体中两相介质的连续方程和Navier-Stoke方程。并尝试用该理论模拟不同水压下渗流引起砂土特性变化的全过程,数值试验得到了流速、渗透系数、孔隙率和砂的流失量等参量的定性变化规律,表明水土相互作用贯穿于渗流的全过程。模拟结果与试验结果相符,且流速与压力梯度满足达西定律,表明应用该理论来模拟复杂的水土相互作用是可行的,它具有一定的理论意义,为进一步研究多孔介质中考虑流固耦合的渗透破坏研究如流土、管涌等奠定了一定的理论基础。     

基于TOUGH2和FLAC3D的水-气二相流-固耦合模型

要实现饱和-非饱和土体的水-气二相流-固耦合作用的数值模拟,关键是如何实现渗流场与应力场的耦合,目前,松弛耦合方法不仅精度能够满足一般的工程要求,而且求解方法易于实现,是多相流-固耦合问题研究的主要方向。本研究根据松弛耦合原理,利用TOUGH2/EOS3来模拟水-气二相渗流过程;利用FLAC3D来模拟土体变形;基于水-气二相渗流过程与力学过程的相互影响,给出饱和度-密度、孔隙水压力和气压力-有效应力、体积应变-孔隙率、孔隙率-固有渗透率,以及孔隙率-毛细压力等状态变量之间的耦合关系式;利用C++控制TOUGH2、FLAC3D在流固耦合过程中的有序运行并利用耦合关系式传递相关状态变量,以实现对水-气二相流-固耦合作用的模拟。利用该模型对一现场压气试验进行模拟,通过模型计算值与实测值的对比,验证了模型的准确性和有效性。     

Prediction of coal mine goaf self-heating with fluid dynamics in porous media

Coal mine goaf self-heating due to exothermic coal oxidation has been recognized as a major threat to coal mine safety. To evaluate the risk of the coal mine goaf self-heating hazard, both gas distribution characteristics in goaf and heat transfer mechanisms in porous media must be studied. However, due to the difficulty of determining goaf permeability and the complexity of the overlaying strata caving characteristics, it is a considerable challenge to determine the thermal-fluid field coupling. Based on the volumetric average method in porous media, this work develops three numerical models for solving fluid dynamics and heat transfer in both longwall face and goaf. The Brinkman-Forchheimer-extended Darcy model is introduced to describe inertia and frictional forces in fluid phase exerted by solid phase. Temperature profile is highly dependent on and affects coal oxidation rate; therefore, governing equations for energy and oxygen mass equilibrium must be coupled. A two-dimensional goaf permeability distribution model is established based on different caving conditions in goafs. Three scenarios are simulated and validated by field and experimental data, and it is observed that these models are capable of predicting gas flow pattern and temperature distribution.     

Untersuchungen zum Strömungsverhalten und zur Lösungskinetik von Gasen im Mehrphasensystem “Aquifer”

In situ remediation techniques like oxygen sparging, air sparging or bio sparging are based on mass transfer reactions between multi-component phases and use the ability of the aquifer to capture gas phases of significant amount (～ 10 % of the pore space). The governing equations of a new kinetic model (variable-volume model), and the results of column experiments on dissolution kinetics of a residual gas phase (one-component gas) are presented. A consistent picture of the vertical gas flow, of the stability and extension of the capillary network used by the gas flow is derived from large-scale lab experiments on gas injection: The comparison between the measured capillary pressure and the calculated one due to the pore size distribution, which was derived from sieve analysis (single-sphere-model), shows that the formation of the capillary network is accompanied by particle rearrangement. Otherwise the small measured capillary pressure are not explainable. The linearity between the gas flux and the dynamical gas phase volume proven by experiments means that an increase in the flux leads to an extension of the capillary network.     

A CFD approach for modeling the rime-ice accretion process on a horizontal-axis wind turbine

The main objective of this paper is to model the rime-ice accretion process on a horizontal-axis wind turbine (HAWT) operating under icing conditions. The model calculation procedure was divided into two stages. The first stage consists of computing the two-phase flow that is composed of air and water phases. The two-phase flow was solved using the Eulerian–Eulerian approach with the help of Fluent, a commercial CFD code. At the end of the calculations, it is possible to obtain the local volume fraction and local water incidence velocity, which can be used to estimate the local collision coefficient. In the second stage, the thickness of the newly accreted ice layer was calculated and the geometry so obtained was processed with an in-house smoothing algorithm. From the model simulation, it is possible to obtain the 3-D ice shape, as well as the ice load on the HAWT blade, which can then be used to analyze the performance of the HAWT under icing conditions, or to help develop a feasible de-icing solution. The important model input parameters include: air speed, median volume diameter (MVD) of water droplets, liquid water content (LWC) and air temperature. The turbine under consideration is a NREL VI turbine that rotates along its axis at a constant angular speed. In the flow domain, a rotating reference frame was defined to account for the angular movement of the wind turbine. The results reveal, among other things, ice geometry information that can otherwise hardly be obtained by on-site observations.     

页岩气藏流固耦合渗流模型及有限元求解

 页岩气渗流模型是页岩气藏动态分析和数值模拟的基础。将裂缝性页岩气藏视为基质孔隙–裂缝双重介质,同时考虑岩石骨架变形对气体渗流场的影响,建立页岩气藏流固耦合渗流模型。模型假设基质孔隙内作克努森流动,裂缝中作达西渗流,综合考虑页岩气壁面滑脱流动与孔内扩散作用、吸附与脱附、应力敏感性等渗流机制。采用有限元法离散控制方程及全隐式耦合求解方法,编制计算机程序。考虑真实页岩参数取值,利用该模型进行算例分析。结果表明,页岩气藏压力下降速率小于常规裂缝性气藏压力下降速率;裂缝渗透率是影响裂缝渗流压力衰减的主要因素,需考虑页岩裂缝导流能力与基质产气速率的匹配关系;原始地层压力越小,裂缝渗流压力衰减越慢。所建模型可为页岩气藏模拟器开发及动态分析提供理论基础。     

弱胶结油藏大孔道出砂的渗流与管流耦合模型

将渗流力学、岩石力学、流体力学知识相结合,建立了弱胶结油藏大孔道出砂的渗流与管流耦合模型,并进行了实例计算。油藏渗流模型主要基于达西定律,并且考虑了储层的可压缩性。将大孔道中的油藏流体与砂粒看作单分散的、稀释的悬浮液。考虑砂粒在大孔道中的沉降及液固两相流之间的相互作用,基于Navier-Stokes方程推导了油藏流体-砂粒两相运动的控制方程。其中,出砂流速门槛值和砂粒剥落本构方程通过实验方法确定,油藏流体渗流与大孔道管流之间通过压力和流量关系进行耦合。最后,通过计算研究了矩形区域内-注-采情况下大孔道的延伸情况,得到了开发过程中流场压力分布图和累计出砂量变化曲线,分析了油藏流体粘度、注采速度对出砂的影响。