迂曲裂縫內(nèi)支撐劑運(yùn)移數(shù)值模擬

胡曉東 張鵬天 周福建 李昕桐 白亞超 陳超

摘要：水力壓裂后，儲(chǔ)層形成較多的迂曲裂縫，支撐劑在迂曲裂縫中的運(yùn)移和鋪置規(guī)律尚不明確，在稠密顆粒流模型的基礎(chǔ)上耦合顆粒物料仿真軟件，模擬迂曲角度、支撐劑尺寸、注入速率、攜砂液密度和支撐劑密度對(duì)迂曲裂縫中支撐劑運(yùn)移和砂堵的影響。結(jié)果表明：大多支撐劑會(huì)堆積在裂縫第一轉(zhuǎn)折角處，迂曲角度的變化會(huì)顯著影響支撐劑的運(yùn)移和砂堵的發(fā)生，迂曲角度大，容易發(fā)生砂堵且砂堵程度嚴(yán)重；密度大和粒徑小的支撐劑發(fā)生砂堵時(shí)會(huì)形成致密的支撐劑堆積，加劇砂堵程度；注入速率大，支撐劑在裂縫中運(yùn)移距離較遠(yuǎn)，但注入速度和攜砂液密度的增加，會(huì)加劇砂堵程度。

關(guān)鍵詞：迂曲裂縫； 支撐劑運(yùn)移； 砂堵預(yù)測(cè)； 數(shù)值模擬

中圖分類(lèi)號(hào)：TE 357.12?? 文獻(xiàn)標(biāo)志碼：A

文章編號(hào)：1673-5005（2024）03-0111-08?? doi：10.3969/j.issn.1673-5005.2024.03.012

Numerical simulation of proppant transport in tortuous fractures

HU Xiaodong1，2， ZHANG Pengtian1，2， ZHOU Fujian1，2， LI Xintong1，2， BAI Yachao1，2， CHEN Chao1，2

（1.State Key Laboratory of Petroleum Resources and Engineering， Beijing 102249， China;2.Unconventional Petroleum Research Institute， China University of Petroleum（Beijing）， Beijing 102249， China）

Abstract： After hydraulic fracturing， numerous tortuous fractures can be generated and formed in the reservoir， and the migration and distribution behavior of proppants in these tortuous fractures is not clear. In this study， in coupling with a simulation software of granular materials， a dense particle flow model was proposed to simulate the effects of tortuosity angle， proppant size， injection rate， carrying fluid density and proppant density on the migration of proppants and sand plugging in tortuous fractures. The results show that most proppants can accumulate at the location of first turning angle of the fracture. The variation in tortuosity angle can significantly affect the migration of proppants and the occurrence of sand plugging， with a higher tortuosity angle leading to more severe sand plugging. Proppants with higher density and smaller particle size tend to form dense proppant accumulation during sand plugging， exacerbating the severity of sand plugging. Higher injection rates can result in longer migration distances of proppants in the fracture， but increasing injection velocity and carrying fluid density can exacerbate the severity of sand plugging.

Keywords： tortuous fractures; proppant transport; sand plugging predictions; numerical simulation

近年來(lái)水力壓裂已經(jīng)成為非常規(guī)增產(chǎn)改造的重要手段，支撐劑在裂縫中的運(yùn)移和分布是水力壓裂的一個(gè)重要部分，由于地層的非均質(zhì)性，水力壓裂形成許多迂曲裂縫，在支撐劑注入過(guò)程不能參照平面平行裂縫中支撐劑運(yùn)移規(guī)律，因此研究支撐劑在迂曲裂縫中的運(yùn)移十分關(guān)鍵［1-3］。水力壓裂形成的裂縫形態(tài)復(fù)雜，如果泵注設(shè)計(jì)不合理，壓裂施工過(guò)程中就會(huì)存在砂堵風(fēng)險(xiǎn)，支撐劑在裂縫中的堆積狀態(tài)決定了裂縫導(dǎo)流能力［4］，最終會(huì)影響采收率。通過(guò)調(diào)研國(guó)內(nèi)外相關(guān)研究，發(fā)現(xiàn)目前對(duì)迂曲裂縫中支撐劑運(yùn)移的研究較少，對(duì)迂曲裂縫內(nèi)運(yùn)移過(guò)程砂堵風(fēng)險(xiǎn)的討論鮮有提及，僅有的一些研究主要采用試驗(yàn)的方法和Fluent數(shù)值模擬兩種途徑。然而，由于試驗(yàn)場(chǎng)地限制和試驗(yàn)條件難以控制且大多數(shù)的試驗(yàn)適用于平直裂縫，試驗(yàn)結(jié)果的準(zhǔn)確性并不理想；采用Fluent進(jìn)行模擬，無(wú)法精確模擬支撐劑顆粒的運(yùn)移，無(wú)法進(jìn)一步進(jìn)行砂堵的模擬。針對(duì)上述問(wèn)題，筆者采用Fluent和EDEM耦合［5］的方法進(jìn)行模擬研究。稠密離散相模型（DDPM）是專(zhuān)門(mén)用來(lái)模擬稠密顆粒流的模型，可以理解為簡(jiǎn)化了的離散元方法（DEM）［6］，該模型考慮了顆粒的空隙率以及碰撞，但是對(duì)于碰撞的計(jì)算進(jìn)行了簡(jiǎn)化，不是計(jì)算真實(shí)的碰撞過(guò)程，因此為了更加真實(shí)地模擬支撐劑的碰撞過(guò)程，筆者在DDPM模型的基礎(chǔ)上耦合EDEM軟件來(lái)模擬支撐劑的運(yùn)移過(guò)程，直觀地觀察支撐劑的運(yùn)移過(guò)程，分析迂曲裂縫內(nèi)支撐劑運(yùn)移的影響因素以及對(duì)應(yīng)工況下裂縫入口壓力特征，為砂堵預(yù)測(cè)預(yù)警提供有效建議。

1 模型的開(kāi)發(fā)與驗(yàn)證

通過(guò)調(diào)研以前的支撐劑運(yùn)移數(shù)值模擬研究，發(fā)現(xiàn)數(shù)值模型主要分為歐拉-拉格朗日模型和歐拉-歐拉模型兩種類(lèi)型。歐拉-拉格朗日模型［7］將流體看作連續(xù)相，顆?？醋麟x散相，用歐拉法處理流體運(yùn)動(dòng)，拉格朗日方法計(jì)算顆粒運(yùn)動(dòng)軌跡。歐拉-拉格朗日模型計(jì)算結(jié)果精確，但耗費(fèi)時(shí)間較長(zhǎng)。歐拉-歐拉模型把流體和顆?？闯梢环N混合的連續(xù)相進(jìn)行計(jì)算，通過(guò)歐拉-歐拉模型可以快速獲得模擬結(jié)果。但是歐拉-歐拉模型把固液作為一個(gè)整體計(jì)算，無(wú)法得到顆粒的運(yùn)動(dòng)軌跡。

Larry等［8］使用透明測(cè)試儀研究了水基凍膠壓裂液中的支撐劑運(yùn)移規(guī)律，但由于試驗(yàn)方案設(shè)計(jì)不足，試驗(yàn)沒(méi)有取得成功。Nolte［9］提出了一種廣泛用于支撐劑運(yùn)移研究的方法。在這個(gè)方法中定義了流體效率，它是用裂縫體積除以注入流體體積，其目的是準(zhǔn)確且較為簡(jiǎn)單地研究支撐劑運(yùn)移過(guò)程。然而，上述兩個(gè)模型都只適用于使用高黏度的攜砂液且不考慮支撐劑沉降過(guò)程。Siddhamshetty等［10］提出了一個(gè)用于頁(yè)巖中低黏度流體的支撐劑運(yùn)移模型，該模型中的支撐劑運(yùn)移是由歐拉-歐拉模型模擬的。由于支撐劑被認(rèn)為是連續(xù)相，其間的相互作用如碰撞和分散不能被準(zhǔn)確模擬。因此使用歐拉-歐拉模型進(jìn)行支撐劑運(yùn)移設(shè)計(jì)不準(zhǔn)確。

在歐拉-拉格朗日模型的基礎(chǔ)上設(shè)計(jì)一個(gè)實(shí)驗(yàn)室規(guī)模的支撐劑運(yùn)移模型。假設(shè)：攜砂液是牛頓流體；攜砂液在裂縫中的濾失忽略不計(jì)。

1.1 數(shù)學(xué)模型

本文中的模型是基于稠密離散相法DDPM［11］建立的。DDPM模型是歐拉-拉格朗日模型的一個(gè)子類(lèi)型模型，被廣泛用于石油工程中的支撐劑運(yùn)移模擬。為了更加真實(shí)地描述和模擬支撐劑在迂曲裂縫的運(yùn)移過(guò)程，通過(guò)DDPM模型接口建立了CFD-DEM耦合模型。CFD-DEM耦合模型可以直觀準(zhǔn)確地模擬支撐劑運(yùn)移問(wèn)題。

流場(chǎng)采用CFD方法求解，流體流動(dòng)過(guò)程方程式［12］為

u=0，（1）

ρut+ρ（u）u=-p+μ2.（2）

式中，ρ為流體密度，kg/m3；u為流體速度，m/s；p為流體壓力，Pa；

μ為流體黏度，Pa·s。

式（1）、（2）稱(chēng)為納維-斯托克斯方程，滿足了流體的質(zhì)量和動(dòng)量守恒。顆粒運(yùn)動(dòng)采用DEM方法求解。該方法將牛頓第二定律應(yīng)用于單個(gè)粒子軌跡的計(jì)算，該軌跡由作用于每個(gè)粒子的力和力矩的總和產(chǎn)生。粒子運(yùn)動(dòng)方程［13］為

mp=dvpdt=Fp，f+Fp，p+Fp，w+Fb，（3）

Ip=dωpdt=rpFt.（4）

式中，mp為粒子質(zhì)量，kg；vp為粒子速度，m/s；t為時(shí)間，s；

Fp，f為周?chē)黧w施加在粒子上的力，N；Fp，p為粒子與粒子間接觸力，N；

Fp，w為粒子與壁面間接觸力，N；Fb為粒子體積力，N；rp為粒子半徑，m；Ft為作用在粒子上的切向力，N；Ip為粒子的轉(zhuǎn)動(dòng)慣量，對(duì)于球形粒子， Ip為（2/5）mpr2p，kg·m2 ；wp為粒子角速度，s-1。

根據(jù)Hertz的法向接觸理論［14］和Mindin等的切向接觸理論［15］，F(xiàn)p，p由法向力和切向力組合表示為

Fp，p=（knδn-γnΔup）+（ktδt-γtΔut）.（5）

式中，k為彈力系數(shù)，N/m；δ為重疊高度，m；γ為阻尼系數(shù)，N·s/m；Δu為接觸點(diǎn)處的相對(duì)速度，m/s;下標(biāo)n和t分別表示法向和切向。彈力系數(shù)和阻尼系數(shù)可以根據(jù)粒子材料特性得出［16］。

流體作用在粒子上的力［17］表示為

Fp，f=∫（-p+μ2u）dΩs.（6）

式中，Ωs為粒子域。

由于攜砂液攜帶支撐劑在裂縫中的流動(dòng)形態(tài)為湍流流態(tài)，因此本文中全部采用標(biāo)準(zhǔn)K-E模型［18］，其中K為紊流脈動(dòng)動(dòng)能，K越大表明紊流脈動(dòng)長(zhǎng)度和時(shí)間尺度越大；E為紊流脈動(dòng)動(dòng)量的耗散率，E越大代表紊流脈動(dòng)長(zhǎng)度和時(shí)間尺度越小。

1.2 物理模型

采用三維迂曲裂縫模型來(lái)模擬迂曲裂縫中的支撐劑運(yùn)移過(guò)程。建立一條長(zhǎng)0.4 m、高0.1 m和寬0.002 m、迂曲角度為45°的裂縫中固液兩相流模型，利用Fluent和EDEM耦合的方法對(duì)迂曲裂縫中的支撐劑運(yùn)移和砂堵預(yù)測(cè)進(jìn)行模擬。圖1為基本案例模型。

1.3 模型驗(yàn)證

首先建立一個(gè)三維垂直平面裂縫模型，用于模型驗(yàn)證。裂縫的長(zhǎng)度、高度和寬度分別為0.381、0.0762和0.002 m，網(wǎng)格總數(shù)為10000，支撐劑粒徑、質(zhì)量分?jǐn)?shù)和密度分別為0.6 mm、1.3%和2650 kg/m3，攜砂液密度和黏度分別為1000 kg/m3和0.001 Pa·s，注入速度為0.3 m/s。圖2為裂縫長(zhǎng)度和高度方向的網(wǎng)格和幾何形狀。在模擬中使用了速度入口和壓力出口。假設(shè)為湍流流態(tài)，壓力-速度耦合與校正采用SIMPLE算法［19］。動(dòng)量和體積分?jǐn)?shù)的空間離散化處理采用QUICK方案 ［20］。為了便于收斂，湍流動(dòng)能方程和湍流耗散率系數(shù)方程采用一階迎風(fēng)格式［21］。

通過(guò)Tong等［22］的試驗(yàn)數(shù)據(jù)和Hu等［23］的模擬結(jié)果來(lái)驗(yàn)證本文的模型。比較注入速度為0.3 m/s時(shí)的支撐劑分布情況，結(jié)果如圖3所示，可以看到該模型中支撐劑床的分布與試驗(yàn)數(shù)據(jù)相似。為了更加清楚地比較結(jié)果，采用無(wú)量綱支撐劑運(yùn)移長(zhǎng)度（支撐劑運(yùn)移長(zhǎng)度與裂縫總長(zhǎng)度之比，DPBL）和無(wú)量綱支撐劑床層高度（最大支撐劑床層高度和總裂縫高度的比，DPBH）［23］。相同條件下，Tong等、Hu等［23］和本文中無(wú)量綱支撐劑運(yùn)移長(zhǎng)度分別為86.27%、95.10%和90.20%，無(wú)量綱支撐劑床層高度分別為50%、45%和57.14%。本文模型結(jié)果和文獻(xiàn)［22］試驗(yàn)結(jié)果及文獻(xiàn)［23］模擬結(jié)果之間的DPBL的最大誤差為4.9%，DPBH的最大誤差為12.14%，說(shuō)明本文的模型與試驗(yàn)結(jié)果有合理的匹配。本文中模型與Tong等的試驗(yàn)誤差主要是由于試驗(yàn)中的分支縫，試驗(yàn)中存在一個(gè)旁路，這導(dǎo)致支撐劑運(yùn)移到二級(jí)裂縫中；與Hu等的模擬結(jié)果存在誤差的原因在于本文的顆粒運(yùn)移計(jì)算采用EDEM軟件，而在文獻(xiàn)［23］中的模擬，顆粒和流體都是用ANSYSFLUENT軟件進(jìn)行計(jì)算的。

2 數(shù)值模擬結(jié)果分析

通過(guò)數(shù)值模擬的方法研究迂曲角度、支撐劑尺寸、注入速率、支撐劑密度、裂縫寬度和支撐劑密度對(duì)迂曲裂縫中支撐劑沉降和砂堵的影響。

2.1 基本案例

基本案例，裂縫的長(zhǎng)度、高度和寬度分別為0.4、0.1和0.002 m，支撐劑粒徑、質(zhì)量分?jǐn)?shù)和密度分別為0.6 mm、1.3%和2650 kg/m3，攜砂液密度和黏度分別為1000 kg/m3和0.001 Pa·s，注入速度為0.3 m/s。圖4為支撐劑運(yùn)移與砂堵情況。由圖4可以看出，攜砂液攜帶的支撐劑最初傾向于在第一段裂縫的下部區(qū)域沉淀并形成支撐劑床。隨著支撐劑的不斷注入，支撐劑在裂縫寬度變窄時(shí)，首先在第一轉(zhuǎn)折角形成支撐劑堆積，隨支撐劑不斷增多最終形成密實(shí)的支撐劑堆積，阻礙后進(jìn)入的支撐劑運(yùn)移，這時(shí)入口壓力激增形成砂堵。

2.2 迂曲角度對(duì)支撐劑運(yùn)移和砂堵的影響

圖5為不同迂曲角度下注入時(shí)間為2 s時(shí)支撐劑的運(yùn)移和砂堵?tīng)顩r。在迂曲角度為0°時(shí)，支撐劑運(yùn)移距離較遠(yuǎn)，由于整個(gè)裂縫寬度相同，支撐劑的堆積面是一個(gè)平穩(wěn)的曲線面，同時(shí)在壓力曲線上，裂縫的入口壓力沒(méi)有劇烈變化，因而沒(méi)有發(fā)生砂堵現(xiàn)象；當(dāng)迂曲角度增到15°時(shí)，可以明顯地看到在兩個(gè)裂縫轉(zhuǎn)折處支撐劑堆積呈現(xiàn)翹角的狀態(tài)，這是由于在裂縫轉(zhuǎn)折處裂縫變窄，流體與支撐劑的流動(dòng)方向和速度發(fā)生改變，但其裂縫入口壓力也沒(méi)有劇烈變化，未出現(xiàn)砂堵現(xiàn)象；當(dāng)迂曲角度增到30°和45°時(shí)，支撐劑大部分堆積在第一段裂縫內(nèi)且在第一個(gè)轉(zhuǎn)折處出現(xiàn)翹角，翹角的出現(xiàn)阻礙了支撐劑的運(yùn)移，同時(shí)裂縫入口壓力都在2 s內(nèi)的某一時(shí)刻發(fā)生激增，出現(xiàn)了砂堵現(xiàn)象。

在迂曲角度為30°和45°的裂縫中，2 s內(nèi)發(fā)生了砂堵現(xiàn)象。通過(guò)計(jì)算可得45°迂曲裂縫入口壓力突變時(shí)間和入口壓力峰值時(shí)間比30°迂曲裂縫的分別提前了34%和29%，而入口壓力峰值卻只比30°迂曲裂縫的高了1.7%。這說(shuō)明迂曲角度會(huì)影響砂堵形成的時(shí)間，對(duì)支撐劑封堵裂縫的嚴(yán)重程度影響較小。迂曲角度越小，支撐劑運(yùn)移長(zhǎng)度越長(zhǎng)，且越不容易形成砂堵。

2.3 注入速率對(duì)支撐劑運(yùn)移和砂堵的影響

在壓裂施工過(guò)程中，注入速率是影響支撐劑運(yùn)移的重要參數(shù)。不同注入速率的支撐劑運(yùn)移與砂堵情況如圖6所示。

當(dāng)注入速率為0.1 m/s時(shí)，可以看出在2 s時(shí)，支撐劑的運(yùn)移尚未到達(dá)裂縫第一轉(zhuǎn)折角處，支撐劑主要堆積在注入口附近，其裂縫入口壓力在2 s內(nèi)沒(méi)有激增變化；注入速率為0.2 m/s時(shí)，支撐劑已經(jīng)運(yùn)移到裂縫第一個(gè)轉(zhuǎn)折角處，且在第一個(gè)轉(zhuǎn)折角前形成了大量沉積，入口壓力在2 s內(nèi)也沒(méi)有發(fā)生激增，說(shuō)明2 s時(shí)雖然支撐劑大量堆積在第一段裂縫中，轉(zhuǎn)折角處幾乎被支撐劑堆滿，但其入口處還可以進(jìn)入流體和支撐劑，不會(huì)出現(xiàn)砂堵現(xiàn)象；注入速率為0.3和0.4 m/s時(shí)，支撐劑雖然有運(yùn)移更遠(yuǎn)距離的趨勢(shì)，但是裂縫發(fā)生轉(zhuǎn)折，使其寬度以及流體、支撐劑運(yùn)移方向和速度改變，只有極少量的支撐劑可以通過(guò)第一個(gè)裂縫轉(zhuǎn)折處進(jìn)入第二段裂縫內(nèi)，同時(shí)堆積在轉(zhuǎn)折處的支撐劑墻阻礙著后進(jìn)入的流體和支撐劑。在入口壓力曲線上，也可以看到0.3 和0.4 m/s的入口壓力均在2 s內(nèi)發(fā)生激增，即在該時(shí)刻發(fā)生了砂堵現(xiàn)象。注入速率為0.4 m/s的裂縫入口壓力開(kāi)始突變時(shí)間和入口壓力峰值時(shí)間比0.3 m/s的提前了71%和73%，入口壓力峰值也比0.3 m/s的提高了38%。這說(shuō)明注入速率的增加在一定程度上確實(shí)可以使支撐劑運(yùn)移的更遠(yuǎn)，但是當(dāng)支撐劑在迂曲角前已形成堆積時(shí)，增加注入速率會(huì)提前砂堵發(fā)生的時(shí)間也會(huì)增加支撐劑封堵裂縫的嚴(yán)重程度。

2.4 支撐劑粒徑對(duì)支撐劑運(yùn)移和砂堵的影響

迂曲角度為45°，支撐劑以0.3 m/s的固定速率注入，入口支撐劑質(zhì)量分?jǐn)?shù)恒定為10%，不同粒徑支撐劑運(yùn)移和砂堵情況如圖7所示。

從圖7中可以看出：粒徑為0.4 mm的支撐劑在通過(guò)裂縫時(shí)，只要有很少的一部分在裂縫第一個(gè)轉(zhuǎn)折角處堆積就會(huì)形成砂堵；粒徑為0.5 和0.6 mm的支撐劑通過(guò)裂縫時(shí)，支撐劑在2 s內(nèi)也會(huì)在裂縫的第一個(gè)轉(zhuǎn)折角度堆積，最終造成砂堵；粒徑為0.7 mm的支撐劑通過(guò)裂縫時(shí)，支撐劑也會(huì)在第一個(gè)轉(zhuǎn)折角處形成堆積，但0.7 mm的支撐劑在2 s內(nèi)雖有堆積，但其未造成砂堵。同時(shí)，隨著支撐劑粒徑的增大，裂縫第一個(gè)轉(zhuǎn)折角前形成的支撐劑堆積增多。說(shuō)明支撐劑粒徑越大，其形成的支撐劑堆積越不致密，其中的縫隙越容易使流體通過(guò)。粒徑為0.4、0.5和0.6 mm的支撐劑入口壓力都在不到2 s發(fā)生了激增，而粒徑0.7 mm的支撐劑入口壓力在2 s內(nèi)幾乎沒(méi)有變化。這也說(shuō)明0.4、0.5 和0.6 mm粒徑的支撐劑在運(yùn)移過(guò)程中都造成了砂堵，而0.7 mm粒徑的支撐劑運(yùn)移過(guò)程中，2 s內(nèi)并未發(fā)生砂堵。同時(shí)支撐劑顆粒越小，入口壓力激增的時(shí)間越靠前，迂曲裂縫越快形成砂堵，同時(shí)支撐劑粒徑越小，封堵裂縫的程度越嚴(yán)重。

2.5 攜砂液密度對(duì)支撐劑運(yùn)移和砂堵的影響

圖8為不同攜砂液密度下發(fā)生砂堵時(shí)的情況。從圖8中可以看出，當(dāng)發(fā)生砂堵前一時(shí)刻，不同攜砂液密度下，支撐劑沉積情況和在第一轉(zhuǎn)折角前形成的支撐劑堆積幾乎相同。攜砂液密度為

1000、950 和900 kg/m3時(shí)入口壓力先發(fā)生突變，且發(fā)生突變的時(shí)間也幾乎相同，而攜砂液密度為850 kg/m3時(shí)入口壓力突變時(shí)間稍滯后。這說(shuō)明在密度較大的攜砂液里支撐劑沉降較慢，支撐劑可以更快地隨攜砂液到達(dá)裂縫第一轉(zhuǎn)折角處，從而造成砂堵。裂縫入口處最大壓力對(duì)應(yīng)的攜砂液密度依次為1000、950、900、850 kg/m3。在迂曲裂縫中，攜砂液密度越大越容易形成砂堵且支撐劑封堵裂縫程度嚴(yán)重。

2.6 支撐劑密度對(duì)支撐劑運(yùn)移和砂堵的影響

裂縫迂曲角度為45°，支撐劑以0.3 m/s的固定速率注入，入口支撐劑質(zhì)量分?jǐn)?shù)恒定為10%，不同支撐劑密度下支撐劑運(yùn)移和砂堵情況如圖9所示。

角處形成支撐劑堆積。隨著支撐劑密度增大，發(fā)生砂堵時(shí)支撐劑在轉(zhuǎn)折角度形成的砂堵堆積越多。這就說(shuō)明支撐劑密度越小，支撐劑在運(yùn)移過(guò)程中越不容易沉降，其運(yùn)移到裂縫轉(zhuǎn)折角處的時(shí)間也就越短，越快速地在裂縫轉(zhuǎn)折角處形成支撐劑堆積，從而造成砂堵。同時(shí)，隨著支撐劑密度的增大，迂曲裂縫入口壓力激增的時(shí)間也在不斷增大。這也說(shuō)明支撐劑密度越大，形成砂堵的時(shí)間越晚，但是迂曲裂縫入口壓力峰值在整體規(guī)律上卻在不斷增大，砂堵裂縫的程度越嚴(yán)重。

3 結(jié) 論

（1）支撐劑在迂曲裂縫中的運(yùn)移過(guò)程中，大多數(shù)支撐劑會(huì)在第一個(gè)轉(zhuǎn)折角處形成沉積，最終形成砂堵。

（2）迂曲角度會(huì)顯著影響迂曲裂縫中支撐劑的運(yùn)移和砂堵現(xiàn)象的發(fā)生，在現(xiàn)場(chǎng)壓裂施工的過(guò)程中應(yīng)盡量避免大角度裂縫的形成。

（3）注入速率越大，支撐劑運(yùn)移的越遠(yuǎn)。但是當(dāng)裂縫轉(zhuǎn)折角處形成支撐劑堆積時(shí)，增大注入速率只會(huì)更快更嚴(yán)重地形成砂堵。

（4）支撐劑的粒徑并不是越小越好。較小的支撐劑形成支撐劑堆積比較致密，不容易使流體和支撐劑通過(guò)，會(huì)加重支撐劑封堵裂縫的程度。

（5）攜砂液的密度對(duì)支撐劑的運(yùn)移影響較小，但密度較大的攜砂液會(huì)加快砂堵的形成，且加重支撐劑封堵裂縫的程度。

（6）支撐劑密度越大，在迂曲裂縫中形成砂堵的時(shí)間越慢。但是支撐劑密度并不是越大越好，較重的支撐劑會(huì)加重裂縫的封堵程度。

參考文獻(xiàn)：

［1］ 張皎生，楊煥英，王晶，等.基于復(fù)雜縫網(wǎng)精細(xì)刻畫(huà)的致密油氣藏水平井多段壓裂數(shù)值模擬技術(shù)［J］.中國(guó)海上油氣，2023，35（4）：103-111.

ZHANG Jiaosheng， YANG Huanying， WANG Jing， et al. A method for numerical simulation of multi-stage fracturing of horizontal well in tight oil and gas reservoirs based on fine characterization of complex fracture network［J］.China Offshore Oil and Gas，2023，35（4）：103-111.

［2］ 敬季昀，蘭夕堂，鄒劍，等.裂縫型儲(chǔ)層復(fù)雜人工裂縫壓后定量診斷方法研究［J］.中國(guó)海上油氣，2023，35（5）：185-192.

JING Jiyun， LAN Xitang， ZOU Jian， et al. Research of quantitative diagnosis of complex hydraulic fractures for fractured reservoirs［J］.China Offshore Oil and Gas，2023，35（5）：185-192.

［3］ 陳遲，郭建春，路千里，等.致密氣藏多尺度支撐機(jī)理研究與應(yīng)用［J］.西南石油大學(xué)學(xué)報(bào)（自然科學(xué)版），2022，44（3）：131-138.

CHEN Chi， GUO Jianchun， LU Qianli， et al. Research and application of multi-scale support mechanism in tight gas reservoir［J］. Journal of Southwest Petroleum University （Science & Technology Edition）， 2022，44（3）：131-138.

［4］ 高根英，張靜嫻，陳金峰，等.支撐裂縫導(dǎo)流能力的影響因素分析［J］.石化技術(shù)，2019，26（2）：91.

GAO Genying， ZHANG Jingxian， CHEN Jinfeng， et al. Analysis of factors influencing the inflow capacity of supporting fractures［J］. Petrochemical Technology， 2019，26（2）：91.

［5］ 王建明，邱欽宇，何訊超.攪拌罐內(nèi)基于EDEM-FLUENT耦合的多相流混合數(shù)值模擬［J］.鄭州大學(xué)學(xué)報(bào)（工學(xué)版），2018，39（5）：79-84.

WANG Jianming， QIU Qinyu， HE Xunchao. Numerical simulation of multiphase flow mixing based on EDEM-FLUENT coupling in a stirred tank［J］. Journal of Zhengzhou University （Engineering Edition）， 2018，39（5）：79-84.

［6］ 張超.基于CFD-DEM耦合的粗糙裂縫輸砂機(jī)理及數(shù)值模擬研究［D］.青島：中國(guó)石油大學(xué)（華東），2021.

ZHANG Chao. Research on sand transport mechanism and numerical simulation of rough fracture based on CFD-DEM coupling ［D］. Qingdao：China University of Petroleum （East China）， 2021.

［7］ 趙小薇.基于歐拉-拉格朗日方法的倒置液固流化床及噴動(dòng)床流動(dòng)特性數(shù)值模擬研究［D］.大慶：東北石油大學(xué)，2021.

ZHAO Xiaowei. Numerical simulation of flow characteristics of inverted liquid-solid fluidized bed and jet bed based on Euler-Lagrangian method［D］. Daqing：Northeastern Petroleum University， 2021.

［8］ LARRY H， HANNAH J， ROBERT R， et al. Dynamic experiments on proppant settling in cross-linked fracturing fluids［R］. SPE 8342， 1979.

［9］ NOLTE K G. Determination of proppant and fluid schedules from fracturing-pressure decline［J］. SPE Production Engineering， 1986，1（4）：255-265.

［10］ SIDDHAMSHETTY P， KWON J S. Modeling of hydraulic fracturing and designing of online pumping schedules to achieve uniform proppant concentration oil reservoirs［J］. Computer and Chemical Engineering， 2018，114（9）：306-317.

［11］ 戴天舒.基于DDPM的發(fā)動(dòng)機(jī)復(fù)式油氣分離器優(yōu)化設(shè)計(jì)［D］.武漢：華中科技大學(xué)，2018.

DAI Tianshu. Optimized design of engine compound oil-gas separator based on DDPM［D］.Wuhan： Huazhong University of Science and Technology， 2018.

［12］ WELLER G， TABOR H. A tensorial approach to computational continuum mechanics using object-oriented techniques［J］. Computers in Physics，1998，12（6）：620-631.

［13］ ZENG Junsheng， LI Heng， ZHANG Dongxiao. Numerical simulation of proppant transport in hydraulic fracture with the upscaling CFD-DEM method［J］. Journal of Natural Gas Science and Engineering，2016，33：264-277.

［14］ HERTZ H. Ueber die beruehrung elastischer koerper［J］. Journal Für die Reine und Angewandte Mathematik， 1882，92：156-171.

［15］ MINDLIN R D， DERIESWICZ H. Elastic spheres in contact under varying oblique forces［J］. Journal of Applied Mechanics， 1953，20：327-344.

［16］ CHRISTOPH K， CHRISTOPH G， ALICE H， et al. Models， algorithms and validation for opensource DEM and CFD-DEM［J］. Progress in Computational Fluid Dynamics， 2012，12（2/3）：140-152.

［17］ ZHOUZ Y， KUANG S B， CHU K W， et al. Discrete particle simulation of particle-fluid flow： model formulations and their applicability［J］. Journal of Fluid Mechanics， 2010，661：482-510.

［18］ 任志安，郝點(diǎn)，謝紅杰.幾種湍流模型及其在FLUENT中的應(yīng)用［J］.化工裝備技術(shù)，2009，30（2）：38-40，44.

REN Zhian， HAO Dian， XIE Hongjie. Several turbulence models and their applications in FLUENT［J］. Chemical Equipment Technology， 2009，30（2）：38-40，44.

［19］ 王為術(shù)，徐維暉，翟肇江，等.PISO算法的實(shí)現(xiàn)及與SIMPLE， SIMPLER， SIMPLEC算法收斂性的比較［J］.華北水利水電學(xué)院學(xué)報(bào)， 2007，28（4）：33-36.

WANG Weishu， XU Weihui， ZHAI Zhaojiang， et al. Implementation of PISO algorithm and comparison of convergence with SIMPLE， SIMPLER and SIMPLEC algorithms［J］. Journal of North China Institute of Water Resources and Hydropower， 2007，28（4）：33-36.

［20］ LEONARD B P， MOKHTARI S. ULTRA-SHARP nonoscillatory convection schemes for high-speed steady multidimensional flow［R］. NASA Lewis Research Center 102568，1990.

［21］ HIRSCH C. Numerical computation of internal and external flows［J］. International Journal of Heat and Fluid Flow， 2007，10（4）：371.

［22］ TONG Songyang， MOHANTY K K. Proppant transport study in fractures with intersections［J］. Fuel， 2016，181：463-477.

［23］ HU Xiaodong， WU Kan， LI Gensheng， et al. Effect of proppant addition schedule on the proppant distribution in a straight fracture for slickwater treatment［J］. Journal of Petroleum Science and Engineering， 2018，167：110-119.

（編輯 李志芬）

基金項(xiàng)目：國(guó)家自然科學(xué)基金聯(lián)合基金項(xiàng)目（U23B2084）

第一作者及通信作者：胡曉東（1990-），男，副研究員，博士，研究方向?yàn)樗毫蚜芽p擴(kuò)展/支撐劑運(yùn)移機(jī)制、水力壓裂診斷與評(píng)估、智能壓裂。E-mail： huxiaodong@cup.edu.cn。

引用格式：胡曉東，張鵬天，周福建，等.迂曲裂縫內(nèi)支撐劑運(yùn)移數(shù)值模擬［J］.中國(guó)石油大學(xué)學(xué)報(bào)（自然科學(xué)版），2024，48（3）：111-118.

HU Xiaodong， ZHANG Pengtian， ZHOU Fujian， et al.