• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    鹽酸對(duì)碳酸鹽巖儲(chǔ)層滲透率的影響

    2021-11-04 07:37:32尹思雨李催芳魏凱董琳琳王晨陽嚴(yán)梁柱
    遼寧化工 2021年10期

    尹思雨,李催芳,魏凱,董琳琳,王晨陽,嚴(yán)梁柱

    鹽酸對(duì)碳酸鹽巖儲(chǔ)層滲透率的影響

    尹思雨,李催芳,魏凱,董琳琳,王晨陽,嚴(yán)梁柱

    (長江大學(xué) 石油工程學(xué)院,湖北 武漢 430100)

    酸化是油氣田開發(fā)過程中有效的增產(chǎn)措施之一,可通過酸巖反應(yīng)改善儲(chǔ)層滲透性,但目前酸巖反應(yīng)的模擬理論較少。為了直觀反映酸化對(duì)儲(chǔ)層滲透性的影響,以鹽酸為主體酸,碳酸鹽巖為儲(chǔ)層,基于多孔介質(zhì)滲流、傳質(zhì)模型與酸化滲流耦合模型,利用有限元分析的方法,建立酸巖反應(yīng)的物理-化學(xué)耦合模型,模擬了鹽酸進(jìn)入儲(chǔ)層后的濃度變化,直觀地反映了酸化對(duì)儲(chǔ)層滲透性的改善程度。結(jié)果表明:由于射孔位置滲流速度較大,對(duì)流作用使鹽酸快速進(jìn)入地層并發(fā)生酸巖反應(yīng),致使地層孔隙度增大,滲透率由原始 8.1×10-4μm2增加到2.053×10-2μm2,儲(chǔ)層滲透性得以改善。該方法證實(shí)了酸化可以有效增滲,能夠?yàn)樗峄O(shè)計(jì)和現(xiàn)場施工提供理論指導(dǎo)。

    酸巖反應(yīng);達(dá)西定律;耦合;滲透率;鹽酸;仿真模擬

    在油田開發(fā)生產(chǎn)的過程中,隨著雜質(zhì)、礦物沉積、反應(yīng)生成的沉淀等各種不利因素的影響,油層孔隙裂縫因污染而造成堵塞,使得原油采收率降 低[1]。常用的解堵技術(shù)包括壓裂、酸化以及近年來發(fā)展起來的壓裂充填防砂、水力振蕩解堵、聲波和超聲波防蠟防垢等技術(shù)[2]。采用哪種解堵技術(shù)主要取決于儲(chǔ)層特征及巖礦特征。碳酸鹽巖儲(chǔ)層一般由天然裂縫和溶洞發(fā)育,且非均質(zhì)性嚴(yán)重,縫洞系統(tǒng)是油氣儲(chǔ)集和滲流的主要場所[4],在進(jìn)行壓裂時(shí)會(huì)導(dǎo)致壓裂液大量濾失。一方面,壓裂液濾失進(jìn)入地層會(huì)導(dǎo)致儲(chǔ)層傷害;另一方面,減產(chǎn)效應(yīng)產(chǎn)生的水力裂縫變短,增產(chǎn)效果減弱,酸化解堵工藝技術(shù)可以有效解決這些問題。其次,碳酸鹽巖油藏礦物組分簡單,主要為灰?guī)r及白云巖,極易與鹽酸反應(yīng),且反應(yīng)物無沉淀[5]。因此,常用鹽酸為主體酸來酸化改造從而解除碳酸鹽巖地層的污染問題。酸化施工時(shí),將酸液通過井筒以不壓裂巖石的壓力注入地層,與地層中的巖石發(fā)生溶蝕反應(yīng),從而提高地層的孔隙度及滲透率,達(dá)到解除污染的效果。酸巖反應(yīng)過程中,一方面其自身與巖石發(fā)生化學(xué)反應(yīng),另一方面,酸液在巖石空隙中流動(dòng)遵循達(dá)西定律,因此研究酸巖反應(yīng)和酸液在儲(chǔ)層巖石孔隙中流動(dòng)時(shí)遵循的達(dá)西定律之間的聯(lián)系以及相互的影響,為提高酸巖反應(yīng)速率以及裂縫導(dǎo)流能力奠定基礎(chǔ),在酸化施工增產(chǎn)中起重要作用。

    1 酸巖反應(yīng)動(dòng)力學(xué)方程的建立

    1.1 酸液的選擇

    酸化解堵一般不使用支撐劑,而是依靠酸液對(duì)裂縫壁面的不均勻刻蝕產(chǎn)生一定的導(dǎo)流能力。因此,酸化壓裂的應(yīng)用通常局限于碳酸鹽巖地層。一般選用鹽酸作為溶蝕酸,一是因?yàn)辂}酸成本低,對(duì)地層的溶蝕能力強(qiáng),二是其反應(yīng)生成物(氯化鈣、氯化鎂及二氧化碳)可溶,不產(chǎn)生沉淀,不會(huì)對(duì)地層造成二次傷害。而其他酸液體系,特別HF類酸液,其中F-的濃度一般都較高,在酸化時(shí),易產(chǎn)生簡單氟化物(膠狀的白色沉淀CaF2)、氟硅酸鹽與氟鋁酸鹽等沉淀。酸壓時(shí)對(duì)裂縫壁面的不均勻刻蝕程度高。通常高濃度鹽酸處理效果往往較好,原因?yàn)椋孩偎嵋河行ё饔脮r(shí)間較長,有效作用范圍大;②單位體積鹽酸可產(chǎn)生的二氧化碳?xì)怏w較多,利于殘酸的返排;③單位體積鹽酸可產(chǎn)生較多的氯化鈣、氯化鎂,提高殘酸的黏度,控制了酸巖反應(yīng)速度,并有利于懸浮、攜帶固體顆粒從地層排出。

    1.2 反應(yīng)速率方程

    碳酸鹽巖主要由灰?guī)r組成,灰?guī)r與鹽酸反應(yīng)的化學(xué)方程為:

    2HCl+CaCO3=CaCl2+CO2↑+H2O。 (1)

    根據(jù)質(zhì)量作用定律,在一定的溫度、壓力下化學(xué)反應(yīng)速率與反應(yīng)物濃度的一定方次的乘積成正比。在酸巖反應(yīng)中巖石的濃度可視為定值,反應(yīng)系統(tǒng)的酸巖反應(yīng)速率可表示為:

    式中:—酸巖反應(yīng)速率,mol·(cm2·s)-1;

    —反應(yīng)速率常數(shù);

    —酸液濃度,mol·L-1;

    —反應(yīng)級(jí)數(shù)。

    對(duì)上式兩邊取對(duì)數(shù),得:

    反應(yīng)速率常數(shù)和反應(yīng)級(jí)數(shù)在一定條件下為常數(shù),對(duì)lg和lg進(jìn)行線性回歸處理,求得和值,從而確定酸巖反應(yīng)動(dòng)力學(xué)方程。

    1.3 H+有效傳質(zhì)系數(shù)

    酸壓時(shí),酸液沿裂縫流動(dòng)反應(yīng),濃度逐漸降低,H+有效傳質(zhì)系數(shù)將發(fā)生變化。多孔介質(zhì)的傳質(zhì)模型為:

    基于奈維-斯托克斯方程和連續(xù)性方程,求解定常條件下酸液旋轉(zhuǎn)流動(dòng)反應(yīng)時(shí)的對(duì)流擴(kuò)散偏微分方程[6],可得到:

    式中:—?dú)潆x子有效傳質(zhì)系數(shù);

    —酸液平均運(yùn)動(dòng)黏度,cm2·s-1;

    —旋轉(zhuǎn)角速度,s-1;

    C—時(shí)刻酸液濃度,mol·L-1。

    設(shè)定鹽酸的濃度為 4.5mol·L-1,得到酸巖反應(yīng)各物質(zhì)濃度變化如圖 1所示。

    圖1 酸巖反應(yīng)速率

    2 建立酸化的物理-化學(xué)耦合模型

    2.1 數(shù)學(xué)模型的建立

    對(duì)于鹽酸對(duì)碳酸鹽巖儲(chǔ)層滲透率的影響的研究,酸液在多孔介質(zhì)中運(yùn)移時(shí)會(huì)使孔隙堵塞物或者骨架巖石發(fā)生溶蝕,也會(huì)使地層孔隙度發(fā)生變化,而孔隙性的改變反過來又會(huì)影響酸液的流動(dòng)及濃度分布。因此,鹽酸對(duì)碳酸鹽巖儲(chǔ)層的酸化過程屬于流-固-化間的耦合問題,各種介質(zhì)間的變化相互影響,變化規(guī)律受滲流場-應(yīng)力場-化學(xué)場耦合關(guān)系的控制。流-固-化耦合模型(HMC模型)間的相互作用關(guān)系如圖 2 所示。

    圖2 流-固-化耦合模型(HMC模型)間的相互作用關(guān)系

    2.1.1 多孔介質(zhì)滲流模型

    小變形條件下,巖體有效應(yīng)力與等效密度、流速、流量間的關(guān)系為:

    式中:m—液相源相,m3·s-1;

    p—巖體體積應(yīng)變量;

    —液相密度,kg·m-3;

    ▽—拉普拉斯算子;

    —速度,m·s-1。

    巖石骨架的連續(xù)性方程為:

    式中:—儲(chǔ)水系數(shù),1/Pa。

    當(dāng)滲流速度較低時(shí),屬層流流動(dòng),流體運(yùn)動(dòng)速度與壓力梯度滿足達(dá)西定律[7],即:

    式中:—滲透率,D;

    —液相黏度,Pa·s;

    —壓力,Pa。

    根據(jù)流體動(dòng)學(xué)理論,聯(lián)立式(6)、式(7)、式(8)得多孔介質(zhì)滲流模型如下:

    2.1.2 多孔介質(zhì)傳質(zhì)模型(對(duì)流、擴(kuò)散)

    幾種基礎(chǔ)物理化學(xué)過程中,分子擴(kuò)散作為最主要的傳輸方式,它的驅(qū)動(dòng)力是離子在多孔介質(zhì)內(nèi)的濃度梯度,當(dāng)液體在多孔介質(zhì)中運(yùn)動(dòng)時(shí),機(jī)械彌散與分子擴(kuò)散是同時(shí)存在的,但當(dāng)流速較小時(shí),分子擴(kuò)散占據(jù)主導(dǎo)作用[8],還有對(duì)流運(yùn)動(dòng)。因此,在簡化了的多孔介質(zhì)傳質(zhì)模型中,可忽略機(jī)械彌散的影響;而電遷移是基于離子電勢的帶電粒子定向移動(dòng)過程[9]。 地下水溶質(zhì)運(yùn)移過程中,不涉及外加電場,故不考慮電遷移過程。只考慮一維對(duì)流擴(kuò)散的地下水溶質(zhì)運(yùn)移方程[10]:

    式中:—濃度,mol·L-1;

    —速度,m·s-1;

    D—擴(kuò)散系數(shù),L2·T-1。

    方程左邊為對(duì)流項(xiàng),它表示濃度分布作不變形輸移;方程右邊為擴(kuò)散項(xiàng),濃度峰值被削減,分布被變形。

    2.1.3 酸化滲流耦合模型

    采用應(yīng)力敏感系數(shù)來比較有效反應(yīng)與滲透率的關(guān)系[11]:

    式中:—有效應(yīng)力,MPa;

    —滲透率,mD。

    巖石的空隙壓縮系數(shù)為[12]:

    對(duì)式(11)、式(12)分別積分得:

    聯(lián)立式(12)、式(13)、式(14),并通過積分得:

    如(15)式所示,孔滲指數(shù)由裂隙系統(tǒng)主導(dǎo)[13],酸化對(duì)地層滲透性的影響:

    式中:—有效應(yīng)力,MPa;

    —滲透率,mD;

    0—初始滲透率,mD;

    —孔隙度;

    ε—初始孔隙度;

    —毛細(xì)管半徑,μm。

    2.2 假設(shè)條件

    考慮水力耦合的儲(chǔ)層頁巖滲透性模擬是建立在水力耦合本構(gòu)方程之上。模擬過程中,頁巖固體顆粒設(shè)為彈性材料,且流體在固體顆粒間的流動(dòng)遵循達(dá)西定律。本節(jié)以彈塑性力學(xué)和流體力學(xué)為基礎(chǔ),結(jié)合有效應(yīng)力原理,建立模擬儲(chǔ)層頁巖滲流特性的水力耦合本構(gòu)方程。

    2.2.1 只考慮水在地層中的流動(dòng)

    綜合考慮流體以及介質(zhì)做出如下假設(shè):

    1)流體在多孔介質(zhì)中的流動(dòng)符合達(dá)西線性滲流規(guī)律流體;

    2)流體進(jìn)入多孔介質(zhì)時(shí)滲流速度在適當(dāng)范圍;

    3)忽略重力的影響,不考慮流固耦合。

    2.2.2 鹽酸在地層中的流動(dòng)與反應(yīng)

    鹽酸與碳酸鹽巖儲(chǔ)層中碳酸鈣會(huì)發(fā)生反應(yīng),考慮流固耦合,在模擬過程中做出如下假設(shè):

    1)考慮氫離子的傳質(zhì)擴(kuò)散;

    2)考慮儲(chǔ)層的各向異性,且認(rèn)為鹽酸體系在地層中的流動(dòng)為平面徑向流;

    3)鹽酸與灰?guī)r反應(yīng)服從與其對(duì)應(yīng)的酸巖反應(yīng)動(dòng)力學(xué)方程;

    4)滲流過程等溫。

    2.3 參數(shù)設(shè)置

    相關(guān)參數(shù)見表1、表2、表3。

    表1 模型尺寸

    表2 地層參數(shù)

    表3 流體參數(shù)

    2.4 幾何模型

    酸液經(jīng)射孔進(jìn)入地層,由于井筒對(duì)稱性,可將酸液進(jìn)入地層簡化為平面研究,相較于空間研究更為便捷,且不影響結(jié)果,射孔兩邊為封閉邊界,其余三邊為地層,是開放邊界,假設(shè)地層為多孔介質(zhì),流體的流動(dòng)不符合達(dá)西定理,流體與巖石會(huì)發(fā)生物理化學(xué)反應(yīng),建立有限元分析模型, 見圖 3 。

    圖3 簡化后幾何模型

    2.5 有限元分析

    本文運(yùn)用有限元分析軟件建立酸化的物理-化學(xué)耦合模型。有限元分析軟件優(yōu)點(diǎn)在于高度的靈活性、強(qiáng)大的求解器和較高的計(jì)算精度,它含有一些內(nèi)嵌的經(jīng)典物理模型,包括單物理場和多個(gè)物理場耦合模型,可以直接調(diào)用物理場進(jìn)行建模[14]。

    1)物理場選取。選取達(dá)西定律物理場和稀物質(zhì)傳遞化學(xué)場。達(dá)西定律物理場用于描述儲(chǔ)層中酸液的流動(dòng);稀物質(zhì)傳遞化學(xué)場用于描述酸液和儲(chǔ)層的反應(yīng)過程;通過流動(dòng)耦合接口描述儲(chǔ)層中酸液的流動(dòng)和反應(yīng)之間的相互作用。

    2)參數(shù)選取原則。根據(jù)調(diào)研得到的實(shí)際儲(chǔ)層和酸液參數(shù)性質(zhì),在有限元軟件中進(jìn)行設(shè)置。

    3)網(wǎng)格劃分原則。采用物理場控制網(wǎng)格,網(wǎng)格樣式為自由三角形,考慮到實(shí)際應(yīng)用和電腦計(jì)算時(shí)間,網(wǎng)格單元大小采用較細(xì)化,完整網(wǎng)格包含 1 274 個(gè)域單元和 90 個(gè)邊界元。

    3 結(jié)果及處理

    有限元分析模擬實(shí)驗(yàn)通過速度與時(shí)間、壓力、濃度的函數(shù)達(dá)到壓力與濃度漸變的目的,整個(gè)實(shí)驗(yàn)中變量是時(shí)間 1 000 s,通過時(shí)間的變化計(jì)算出壓力、濃度和滲透率的變化。由于條件限制,忽略模擬過程,僅展示最終實(shí)驗(yàn)結(jié)果[15]。

    3.1 水在在多孔介質(zhì)中的滲流模擬結(jié)果

    設(shè)定多孔介質(zhì)的孔隙率為 0.2,滲透率為 8.1×10-4μm2,得到水在多孔介質(zhì)中的滲流模型,見圖 4。

    圖 4 水在多孔介質(zhì)中的滲流模型

    由圖 4 可知,水在多孔介質(zhì)內(nèi)的滲流情況以射孔為中心向四周輻射,距離越遠(yuǎn)速率越低。

    3.2 鹽酸在多孔介質(zhì)中的滲流結(jié)果

    設(shè)碳酸鹽巖儲(chǔ)層的初始滲透率為 8.1×10-4μm2,注入儲(chǔ)層的鹽酸的濃度為4.5 mol·m-3,將鹽酸溶液從射孔處注入,得到儲(chǔ)層中HCl濃度分布和CaCl2濃度分布分別見圖 5和圖 6。

    圖 5 HCl濃度分布

    由圖 5可知,在射孔位置處,HCl濃度最高,越遠(yuǎn)離射孔的位置,HCl濃度越低。

    由圖 6 可知,靠近射孔位置生成CaCl2濃度較低,是由于考慮了滲流作用引起的。

    3.3 酸化對(duì)滲透性的影響

    設(shè)碳酸鹽巖儲(chǔ)層的初始滲透率為 8.1×10-4μm2,注入鹽酸的濃度為 4.5 mol·m-3,反應(yīng)時(shí)間為 1 000 s,得到酸化作用下滲透率分布見圖 7。

    由圖 7 可知,由于射孔位置滲流速度較大,對(duì)流作用使HCl快速進(jìn)入地層并發(fā)生酸巖反應(yīng),致使地層孔隙度增大,從而改善其滲透特性。

    圖6 CaCl2濃度分布

    圖7 酸化作用下滲透率分布

    碳酸鹽巖儲(chǔ)層滲透率隨時(shí)間的變化關(guān)系如 圖 8 所示。

    圖8 碳酸鹽巖儲(chǔ)層滲透率隨時(shí)間的變化關(guān)系

    由圖 8 可知,隨著反應(yīng)時(shí)間的增加,碳酸鹽巖的滲透率也隨之增大,最后趨向于穩(wěn)定。

    4 結(jié) 論

    以鹽酸作為主體酸,碳酸鹽巖作為儲(chǔ)層,利用有限元分析的方法,直觀地反映了酸液對(duì)儲(chǔ)層滲透性的改變,此次模擬基于多孔介質(zhì)滲流、傳質(zhì)模型與酸化滲流耦合模型,以數(shù)學(xué)模型為基礎(chǔ),表明了該方法的可行性。得到的結(jié)論如下:

    1)研究分析了酸化對(duì)儲(chǔ)層滲透性的影響,得到了酸化作用下滲透率分布圖。通過有限元的方法仿真模擬了流體在多孔介質(zhì)中的滲流,得到了水在多孔介質(zhì)中的滲流模型。可知水在多孔介質(zhì)內(nèi)的滲流情況,以射孔為中心向四周輻射,距離越遠(yuǎn)速率越低。

    2)結(jié)合酸化對(duì)儲(chǔ)層滲透性的影響,將水換成鹽酸,最終建立了鹽酸與碳酸鹽巖的物理-化學(xué)耦合模型,可知儲(chǔ)層中鹽酸濃度分布和CaCl2濃度分布,靠近射孔位置生成物濃度較低,是由于考慮了滲流作用引起的。

    3)由于射孔位置滲流速度較大,對(duì)流作用使HCl快速進(jìn)入地層并發(fā)生酸巖反應(yīng),致使地層孔隙度增大,從而改善其滲透特性。隨著反應(yīng)時(shí)間的增加,碳酸鹽巖的滲透率也隨之增大,最后趨向于穩(wěn)定。

    基于碳酸鹽巖儲(chǔ)層,相較于砂巖的成分較為簡單,本方法并不適用于砂巖和其他巖類。

    [1]梁春,許小龍. 水力脈沖解堵技術(shù)的研究與應(yīng)用[C]. 吉林省科學(xué)技術(shù)協(xié)會(huì).增強(qiáng)自主創(chuàng)新能力 促進(jìn)吉林經(jīng)濟(jì)發(fā)展——啟明杯?吉林省第四屆科學(xué)技術(shù)學(xué)術(shù)年會(huì)論文集(上冊).吉林省科學(xué)技術(shù)協(xié)會(huì):吉林省科學(xué)技術(shù)協(xié)會(huì)學(xué)會(huì)學(xué)術(shù)部,2006.

    [2]王飛航,李力,張建坤,等.淺談?dòng)途霎a(chǎn)技術(shù)[J].石化技術(shù),2018,25(4):248.

    [3]孫淑娟,張群正,劉金磊,等.碳酸鹽儲(chǔ)層酸化解堵配方優(yōu)化及解堵機(jī)理[J].油田化學(xué),2017,34(4):585-589.

    [4]魏星,王永剛,鄭淑媛.碳酸鹽巖儲(chǔ)層酸化工作液現(xiàn)狀研究[J].西部探礦工程,2008(8):63-65.

    [5]姚軍,劉丕養(yǎng),黃朝琴,等.碳酸鹽巖油藏酸巖反應(yīng)流動(dòng)模擬研究進(jìn)展[J].地球科學(xué),2017,42(8):1263-1272.

    [6]咼義,劉平禮,吳雙,等.酸巖反應(yīng)動(dòng)力學(xué)實(shí)驗(yàn)影響因素研究[J].重慶科技學(xué)院學(xué)報(bào)(自然科學(xué)版),2014,16(5):17-19.

    [7]茹忠亮,簡闊,馬國勝.考慮Klinkenberg效應(yīng)的多孔介質(zhì)氣體滲流模型[J].中國科技論文,2018,13(3):310-313.

    [8]徐玉璐. 多孔介質(zhì)中污染物運(yùn)移及彌散系數(shù)對(duì)流速依賴性實(shí)驗(yàn)研究[D].合肥:合肥工業(yè)大學(xué),2017.

    [9]吳靜新. 混凝土內(nèi)氯離子傳輸模型及其數(shù)值研究[D].哈爾濱:哈爾濱工業(yè)大學(xué),2012.

    [10]楊國錄.對(duì)流擴(kuò)散問題的新解法[J].水利學(xué)報(bào),1993(3):56-63.

    [11]張睿,寧正福,楊峰,等.頁巖應(yīng)力敏感實(shí)驗(yàn)與機(jī)理[J].石油學(xué)報(bào),2015,36(2):224-231.

    [12]DAVID C, WONG T F, ZHU W L,et al. Laboratory measurement of compaction-induced permeability change in porous rocks: Implications for the generation and maintenance of pore pressure excess in the crust[J]., 1994, 143(1-3) : 425-456.

    [13]劉永茜. 循環(huán)載荷作用下煤體滲透率演化的實(shí)驗(yàn)分析[J]. 煤炭學(xué)報(bào),2019,44(8):2579-2588.

    [14]李宏遠(yuǎn),唐海,呂棟梁,等. 基于COMSOL軟件的巖心竄流系數(shù)研究[J]. 石油化工應(yīng)用,2013,32(11):9-11.

    [15] 柳崎,索麗敏.地下巖體內(nèi)多孔介質(zhì)中裂隙流運(yùn)移過程的COMSOL Multiphysics仿真模擬[J]. 系統(tǒng)仿真技術(shù),2019,15(3): 184-187.

    Influence of Hydrochloric Acid on the Permeability of Carbonate Reservoirs

    ,,,,,

    (School of Petroleum Engineering, Yangtze University, Wuhan Hubei 430100, China)

    Acidification is one of the effective production stimulation measures in the development of oil and gas fields. It can improve the permeability of the reservoir through acid-rock reaction, but there are few simulation theories for acid-rock reaction at present. In order to intuitively reflect the impact of acidification on reservoir permeability, hydrochloric acid was used as the main acid and carbonate rock was used as the reservoir, based on the coupling model of porous media seepage, mass transfer model and acidification seepage, the physical-chemical coupling model of acid-rock reaction was established by finite element analysis. The physical-chemical coupling model of the reaction simulated the concentration change of hydrochloric acid after entering the reservoir, and intuitively reflected the degree of improvement of the permeability of the reservoir by acidification. The example simulation results showed that, due to the high seepage velocity at the perforation position, the convection made the hydrochloric acid quickly enter the formation and the acid-rock reaction occurred, which increased the porosity of the formation, increased the permeability from 8.1×10-4μm2to 2.053×10-2μm2. This method proves that acidification can effectively increase permeability and can provide theoretical guidance for acidification design and on-site construction.

    Acid-rock reaction; Darcy's law; Coupling; Permeability; Hydrochloric acid; Simulation

    中國石油集團(tuán)重大專項(xiàng),陸相中高成熟度頁巖油勘探開發(fā)關(guān)鍵技術(shù)研究與應(yīng)用(項(xiàng)目編號(hào):2019E-2605);湖北省大學(xué)生創(chuàng)新創(chuàng)業(yè)訓(xùn)練計(jì)劃(項(xiàng)目編號(hào):S202010489034);長江大學(xué)大學(xué)生創(chuàng)新創(chuàng)業(yè)訓(xùn)練計(jì)劃(項(xiàng)目編號(hào):2019307)。

    2021-03-20

    尹思雨(2000-),女,湖北省宜昌市人,研究方向:石油工程。

    魏凱(1983-),男,副教授,研究方向:石油工程數(shù)值模擬技術(shù)、智能鉆井風(fēng)險(xiǎn)評(píng)價(jià)理論與方法。

    TE344

    A

    1004-0935(2021)10-1523-06

    99re6热这里在线精品视频| 日韩欧美精品免费久久| 中文字幕最新亚洲高清| 黄色怎么调成土黄色| 青春草亚洲视频在线观看| 亚洲欧洲日产国产| 国产精品久久久久久久电影| 日本-黄色视频高清免费观看| 国产精品国产三级专区第一集| 九九在线视频观看精品| 国产成人精品在线电影| 精品国产露脸久久av麻豆| 九九在线视频观看精品| 国产福利在线免费观看视频| 日韩熟女老妇一区二区性免费视频| 一级毛片黄色毛片免费观看视频| 男女下面插进去视频免费观看 | 国精品久久久久久国模美| 一本久久精品| 精品亚洲成a人片在线观看| 黄色配什么色好看| 精品国产乱码久久久久久小说| 欧美另类一区| 国产日韩欧美视频二区| 蜜桃在线观看..| 纵有疾风起免费观看全集完整版| 免费久久久久久久精品成人欧美视频 | 夜夜骑夜夜射夜夜干| 男人舔女人的私密视频| 十八禁网站网址无遮挡| 五月伊人婷婷丁香| 一区二区三区四区激情视频| 亚洲精品国产av蜜桃| 免费在线观看完整版高清| 天天操日日干夜夜撸| 十分钟在线观看高清视频www| 91精品国产国语对白视频| 青春草国产在线视频| 另类亚洲欧美激情| 2022亚洲国产成人精品| 黄片播放在线免费| 黄色 视频免费看| 欧美激情 高清一区二区三区| 一级片免费观看大全| 精品久久蜜臀av无| 最近中文字幕2019免费版| 18在线观看网站| 国产精品三级大全| 五月天丁香电影| 亚洲国产毛片av蜜桃av| 男女下面插进去视频免费观看 | 亚洲,欧美,日韩| 黄片无遮挡物在线观看| 久久 成人 亚洲| 在线亚洲精品国产二区图片欧美| 色网站视频免费| 七月丁香在线播放| 丝袜脚勾引网站| 成年动漫av网址| 中文天堂在线官网| 免费大片18禁| 欧美成人精品欧美一级黄| 欧美激情极品国产一区二区三区 | av黄色大香蕉| 边亲边吃奶的免费视频| 亚洲国产精品一区三区| 在线观看人妻少妇| 国产一区二区三区综合在线观看 | 久久国产精品男人的天堂亚洲 | 欧美成人精品欧美一级黄| 大香蕉久久网| 国产欧美另类精品又又久久亚洲欧美| 亚洲人与动物交配视频| 国产在线视频一区二区| 视频在线观看一区二区三区| 爱豆传媒免费全集在线观看| 亚洲精品aⅴ在线观看| 人体艺术视频欧美日本| 国产免费福利视频在线观看| 黄片无遮挡物在线观看| 免费播放大片免费观看视频在线观看| 2018国产大陆天天弄谢| 久久精品国产综合久久久 | 18+在线观看网站| 成人毛片a级毛片在线播放| 免费播放大片免费观看视频在线观看| 亚洲国产精品一区三区| 狂野欧美激情性bbbbbb| 精品人妻一区二区三区麻豆| 在线精品无人区一区二区三| 美女大奶头黄色视频| 亚洲色图 男人天堂 中文字幕 | 午夜91福利影院| 老熟女久久久| 亚洲av电影在线进入| 日本av免费视频播放| 久久精品国产自在天天线| 久久久国产一区二区| 亚洲欧洲日产国产| 国产在视频线精品| 女的被弄到高潮叫床怎么办| 久久精品aⅴ一区二区三区四区 | 一本久久精品| 欧美丝袜亚洲另类| 国产又色又爽无遮挡免| 亚洲av日韩在线播放| 国产一区二区三区av在线| 哪个播放器可以免费观看大片| 少妇人妻久久综合中文| 亚洲欧洲精品一区二区精品久久久 | 9热在线视频观看99| 男人舔女人的私密视频| 亚洲国产毛片av蜜桃av| 亚洲性久久影院| 日本-黄色视频高清免费观看| 午夜av观看不卡| 亚洲av福利一区| 亚洲婷婷狠狠爱综合网| 亚洲精品中文字幕在线视频| 夜夜爽夜夜爽视频| 欧美老熟妇乱子伦牲交| 成人毛片60女人毛片免费| 极品人妻少妇av视频| 久久99热这里只频精品6学生| 国产 精品1| 成人无遮挡网站| 欧美性感艳星| 国产黄色视频一区二区在线观看| 美女中出高潮动态图| 免费av不卡在线播放| 久久99蜜桃精品久久| 亚洲av电影在线观看一区二区三区| 韩国av在线不卡| 久久久欧美国产精品| 男的添女的下面高潮视频| 欧美精品国产亚洲| 在现免费观看毛片| 亚洲伊人色综图| 免费大片黄手机在线观看| 国产黄色免费在线视频| 日日撸夜夜添| 激情五月婷婷亚洲| 人妻人人澡人人爽人人| 国产爽快片一区二区三区| 久热这里只有精品99| 少妇 在线观看| 国产乱人偷精品视频| av有码第一页| 国产一区二区三区av在线| 成人漫画全彩无遮挡| av不卡在线播放| 国产精品久久久av美女十八| 成年美女黄网站色视频大全免费| 国产麻豆69| av国产精品久久久久影院| 99热6这里只有精品| 亚洲欧美成人精品一区二区| 午夜日本视频在线| 男女边吃奶边做爰视频| 久久久久久久久久人人人人人人| 在线天堂中文资源库| 亚洲一区二区三区欧美精品| 欧美精品国产亚洲| 曰老女人黄片| 欧美精品高潮呻吟av久久| 狠狠婷婷综合久久久久久88av| 精品一区二区免费观看| 久久韩国三级中文字幕| 狠狠精品人妻久久久久久综合| 成年女人在线观看亚洲视频| 一级毛片黄色毛片免费观看视频| 一区二区av电影网| 欧美日韩视频高清一区二区三区二| 国产精品一区www在线观看| 午夜精品国产一区二区电影| 欧美激情 高清一区二区三区| 人体艺术视频欧美日本| 亚洲av成人精品一二三区| 一二三四中文在线观看免费高清| 韩国高清视频一区二区三区| 80岁老熟妇乱子伦牲交| 亚洲第一区二区三区不卡| 国产精品成人在线| 久久久国产一区二区| 男人爽女人下面视频在线观看| 涩涩av久久男人的天堂| 99re6热这里在线精品视频| 高清黄色对白视频在线免费看| 日本色播在线视频| 亚洲一级一片aⅴ在线观看| 国产成人精品久久久久久| 免费观看a级毛片全部| 男女国产视频网站| 三级国产精品片| 黄色配什么色好看| 亚洲美女视频黄频| 最近手机中文字幕大全| 国产亚洲精品第一综合不卡 | 精品第一国产精品| 一级黄片播放器| 精品少妇黑人巨大在线播放| av免费观看日本| 亚洲成人手机| 亚洲美女搞黄在线观看| 午夜福利视频在线观看免费| 久久久久久久久久人人人人人人| 亚洲精品色激情综合| 美国免费a级毛片| 日本91视频免费播放| 国产av精品麻豆| 一本色道久久久久久精品综合| 久久热在线av| 天美传媒精品一区二区| 在现免费观看毛片| 天堂俺去俺来也www色官网| 亚洲少妇的诱惑av| 久久韩国三级中文字幕| 十八禁高潮呻吟视频| 99国产综合亚洲精品| 爱豆传媒免费全集在线观看| 精品第一国产精品| 亚洲成人一二三区av| 午夜免费男女啪啪视频观看| 只有这里有精品99| 在线观看美女被高潮喷水网站| 免费黄色在线免费观看| 在线观看人妻少妇| 国产成人一区二区在线| 国产免费一级a男人的天堂| 两性夫妻黄色片 | 午夜久久久在线观看| 亚洲精品久久午夜乱码| 99久久精品国产国产毛片| 丰满饥渴人妻一区二区三| 亚洲国产毛片av蜜桃av| 高清视频免费观看一区二区| 国产熟女欧美一区二区| 亚洲国产av影院在线观看| 有码 亚洲区| 久久韩国三级中文字幕| 国产乱来视频区| 亚洲国产欧美日韩在线播放| 亚洲伊人色综图| 久久毛片免费看一区二区三区| 九草在线视频观看| 亚洲国产精品一区二区三区在线| 制服诱惑二区| 久久久久久伊人网av| 成人18禁高潮啪啪吃奶动态图| 久久午夜福利片| 秋霞在线观看毛片| 久久久a久久爽久久v久久| 久久亚洲国产成人精品v| 久久国产精品大桥未久av| 久久人妻熟女aⅴ| 69精品国产乱码久久久| 欧美精品国产亚洲| 久久精品久久精品一区二区三区| 久久久久久人人人人人| 日韩在线高清观看一区二区三区| 免费观看av网站的网址| 久久久亚洲精品成人影院| 亚洲熟女精品中文字幕| 欧美精品国产亚洲| 国产成人a∨麻豆精品| 免费日韩欧美在线观看| 热re99久久国产66热| 18禁裸乳无遮挡动漫免费视频| 日韩大片免费观看网站| av天堂久久9| 高清欧美精品videossex| 亚洲av免费高清在线观看| 国产片特级美女逼逼视频| 秋霞在线观看毛片| 国产色爽女视频免费观看| 捣出白浆h1v1| 久久国内精品自在自线图片| 下体分泌物呈黄色| 18禁观看日本| 另类亚洲欧美激情| 免费女性裸体啪啪无遮挡网站| 日韩中字成人| 伦理电影免费视频| av不卡在线播放| 妹子高潮喷水视频| 免费少妇av软件| 五月开心婷婷网| 免费人妻精品一区二区三区视频| 日韩制服骚丝袜av| 色视频在线一区二区三区| 最后的刺客免费高清国语| 在线观看www视频免费| 国产淫语在线视频| 免费女性裸体啪啪无遮挡网站| 中文欧美无线码| 又粗又硬又长又爽又黄的视频| 久久 成人 亚洲| 性色av一级| 国产xxxxx性猛交| 国产男女超爽视频在线观看| 国产精品成人在线| 男女国产视频网站| 99久国产av精品国产电影| 熟妇人妻不卡中文字幕| 日韩,欧美,国产一区二区三区| 男的添女的下面高潮视频| 少妇人妻久久综合中文| 天美传媒精品一区二区| 97精品久久久久久久久久精品| 国产一区有黄有色的免费视频| 国产一区二区三区综合在线观看 | 国精品久久久久久国模美| av女优亚洲男人天堂| av黄色大香蕉| 久久毛片免费看一区二区三区| 一级毛片黄色毛片免费观看视频| 国产成人91sexporn| 国产精品女同一区二区软件| 国产精品久久久久久精品电影小说| av在线app专区| 日韩熟女老妇一区二区性免费视频| av在线观看视频网站免费| 国产精品免费大片| 春色校园在线视频观看| 欧美日韩亚洲高清精品| 宅男免费午夜| 亚洲av日韩在线播放| 各种免费的搞黄视频| 91精品伊人久久大香线蕉| 免费看av在线观看网站| 久久久精品免费免费高清| 国产有黄有色有爽视频| 欧美日韩视频精品一区| 成人影院久久| 久久久久精品性色| 精品亚洲成国产av| 国产高清不卡午夜福利| 久久久久精品久久久久真实原创| 水蜜桃什么品种好| 看十八女毛片水多多多| 国产视频首页在线观看| 飞空精品影院首页| 精品国产露脸久久av麻豆| 久久ye,这里只有精品| 新久久久久国产一级毛片| 亚洲美女黄色视频免费看| 亚洲四区av| a级片在线免费高清观看视频| 日韩av不卡免费在线播放| 在线观看美女被高潮喷水网站| 国产色爽女视频免费观看| 91aial.com中文字幕在线观看| 国产亚洲精品久久久com| 亚洲久久久国产精品| 日韩中文字幕视频在线看片| 国内精品宾馆在线| 午夜福利,免费看| 亚洲伊人久久精品综合| 久久国内精品自在自线图片| 亚洲,欧美精品.| 99热网站在线观看| 妹子高潮喷水视频| 永久免费av网站大全| 午夜免费鲁丝| 国产淫语在线视频| 日韩中字成人| 亚洲欧美清纯卡通| av视频免费观看在线观看| 日韩精品有码人妻一区| 午夜激情av网站| 十八禁网站网址无遮挡| 日韩中文字幕视频在线看片| 秋霞在线观看毛片| 亚洲第一区二区三区不卡| 日韩成人av中文字幕在线观看| 少妇精品久久久久久久| 一区二区三区精品91| 亚洲欧美清纯卡通| 亚洲国产精品999| 日韩三级伦理在线观看| 激情视频va一区二区三区| 插逼视频在线观看| 两个人看的免费小视频| 精品亚洲乱码少妇综合久久| 欧美精品人与动牲交sv欧美| 黄色怎么调成土黄色| 国产一区有黄有色的免费视频| 午夜av观看不卡| 久久精品国产综合久久久 | 满18在线观看网站| 欧美精品av麻豆av| 免费黄网站久久成人精品| 亚洲国产精品一区二区三区在线| 青春草亚洲视频在线观看| 国产日韩欧美亚洲二区| 国产极品天堂在线| 好男人视频免费观看在线| 日韩一区二区视频免费看| 国产极品粉嫩免费观看在线| 国产精品久久久久久精品电影小说| 久久久久国产网址| 久久这里只有精品19| 亚洲五月色婷婷综合| 亚洲经典国产精华液单| 欧美另类一区| 久久精品国产鲁丝片午夜精品| 男女边摸边吃奶| 亚洲高清免费不卡视频| 日韩欧美一区视频在线观看| 亚洲国产av新网站| 欧美精品亚洲一区二区| 亚洲在久久综合| 国产麻豆69| 免费高清在线观看日韩| 日本猛色少妇xxxxx猛交久久| 成人二区视频| 中文字幕另类日韩欧美亚洲嫩草| 国产高清国产精品国产三级| 亚洲av在线观看美女高潮| 亚洲精品自拍成人| 国产在视频线精品| 免费大片黄手机在线观看| 国产精品免费大片| 韩国精品一区二区三区 | 日韩免费高清中文字幕av| 亚洲精华国产精华液的使用体验| 亚洲av成人精品一二三区| 秋霞在线观看毛片| 国产熟女欧美一区二区| 2018国产大陆天天弄谢| 熟妇人妻不卡中文字幕| 亚洲av国产av综合av卡| 国产成人免费无遮挡视频| 一区二区日韩欧美中文字幕 | 好男人视频免费观看在线| 精品一品国产午夜福利视频| 国产极品粉嫩免费观看在线| 韩国高清视频一区二区三区| 99re6热这里在线精品视频| 日韩制服骚丝袜av| 亚洲成人一二三区av| 伊人亚洲综合成人网| videossex国产| 日韩,欧美,国产一区二区三区| 在线观看三级黄色| 亚洲精品国产av蜜桃| 亚洲av中文av极速乱| 一二三四中文在线观看免费高清| 午夜视频国产福利| 国产欧美另类精品又又久久亚洲欧美| 纯流量卡能插随身wifi吗| 日韩av在线免费看完整版不卡| videossex国产| 精品少妇黑人巨大在线播放| 一边亲一边摸免费视频| 国产深夜福利视频在线观看| 亚洲av中文av极速乱| 欧美人与性动交α欧美软件 | 午夜免费鲁丝| 国产亚洲精品第一综合不卡 | 大香蕉97超碰在线| 久久久久久久久久久免费av| 综合色丁香网| 午夜视频国产福利| 777米奇影视久久| 久久久久久人人人人人| 日韩一本色道免费dvd| 亚洲av免费高清在线观看| 亚洲经典国产精华液单| 青春草亚洲视频在线观看| 亚洲精品av麻豆狂野| 在线观看免费视频网站a站| 精品久久久久久电影网| 激情视频va一区二区三区| 成人毛片60女人毛片免费| 大片免费播放器 马上看| 久久av网站| 激情五月婷婷亚洲| 又黄又爽又刺激的免费视频.| 婷婷色综合大香蕉| 性色av一级| 亚洲精品乱码久久久久久按摩| 欧美 亚洲 国产 日韩一| 老司机影院毛片| 亚洲精品av麻豆狂野| 尾随美女入室| 国产高清国产精品国产三级| 午夜免费观看性视频| 少妇的逼好多水| 高清毛片免费看| 国产精品无大码| 国产乱来视频区| 午夜日本视频在线| 一本大道久久a久久精品| 性高湖久久久久久久久免费观看| av一本久久久久| 男人添女人高潮全过程视频| 极品人妻少妇av视频| 啦啦啦啦在线视频资源| 亚洲精品日韩在线中文字幕| 日本欧美国产在线视频| 久久女婷五月综合色啪小说| 高清不卡的av网站| 色5月婷婷丁香| 最近最新中文字幕免费大全7| 亚洲精品国产av蜜桃| 日本av免费视频播放| 侵犯人妻中文字幕一二三四区| 大话2 男鬼变身卡| 最近手机中文字幕大全| 欧美最新免费一区二区三区| 国产av国产精品国产| 国产伦理片在线播放av一区| 91久久精品国产一区二区三区| 久久精品熟女亚洲av麻豆精品| 最近最新中文字幕免费大全7| 国产爽快片一区二区三区| 五月伊人婷婷丁香| 飞空精品影院首页| 欧美日韩成人在线一区二区| 亚洲国产色片| 成年美女黄网站色视频大全免费| 亚洲国产最新在线播放| 看免费成人av毛片| 欧美另类一区| 人妻 亚洲 视频| 各种免费的搞黄视频| 亚洲欧美色中文字幕在线| 自拍欧美九色日韩亚洲蝌蚪91| 欧美 日韩 精品 国产| 视频在线观看一区二区三区| 丝袜脚勾引网站| 欧美精品人与动牲交sv欧美| 2021少妇久久久久久久久久久| 欧美最新免费一区二区三区| 王馨瑶露胸无遮挡在线观看| 国产成人精品在线电影| 久久韩国三级中文字幕| 99re6热这里在线精品视频| 校园人妻丝袜中文字幕| 色5月婷婷丁香| 美女xxoo啪啪120秒动态图| kizo精华| 伦理电影免费视频| 国产av一区二区精品久久| xxx大片免费视频| av在线老鸭窝| 日日撸夜夜添| 男女午夜视频在线观看 | 丰满饥渴人妻一区二区三| 2022亚洲国产成人精品| 午夜激情av网站| 22中文网久久字幕| 日本欧美国产在线视频| 午夜福利在线观看免费完整高清在| 少妇猛男粗大的猛烈进出视频| 欧美日本中文国产一区发布| 久久综合国产亚洲精品| 18禁国产床啪视频网站| 激情五月婷婷亚洲| 亚洲三级黄色毛片| 亚洲精品色激情综合| 老女人水多毛片| 99热网站在线观看| 伊人久久国产一区二区| 我的女老师完整版在线观看| 久久久久人妻精品一区果冻| 国产淫语在线视频| 久久久精品94久久精品| 午夜免费男女啪啪视频观看| 黑人猛操日本美女一级片| 亚洲精品一二三| 久久综合国产亚洲精品| 精品少妇黑人巨大在线播放| 免费大片18禁| 欧美人与善性xxx| 亚洲精品久久午夜乱码| 看免费成人av毛片| 一区二区三区四区激情视频| 精品酒店卫生间| 婷婷色综合www| 亚洲av中文av极速乱| 午夜精品国产一区二区电影| 国产成人免费观看mmmm| 美女中出高潮动态图| 青青草视频在线视频观看| 99国产精品免费福利视频| 日韩人妻精品一区2区三区| 麻豆乱淫一区二区| 香蕉国产在线看| 乱人伦中国视频| 97在线视频观看| 伦理电影大哥的女人| 五月天丁香电影| av播播在线观看一区| 久久精品国产亚洲av天美| 国产精品欧美亚洲77777| 大香蕉久久网| 美女国产高潮福利片在线看| 日本色播在线视频| 在线观看一区二区三区激情| 中国国产av一级| 亚洲国产日韩一区二区| 18禁在线无遮挡免费观看视频| 高清黄色对白视频在线免费看| 99热这里只有是精品在线观看| 久久久久精品久久久久真实原创| 免费少妇av软件| 青青草视频在线视频观看| 精品国产一区二区三区久久久樱花| 日韩 亚洲 欧美在线| 日韩精品有码人妻一区| 午夜福利乱码中文字幕| 国产成人免费无遮挡视频| 男的添女的下面高潮视频|