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

    管單元用于縮減巖體滲透張量計(jì)算規(guī)模

    2022-01-24 02:16:54王俊奇劉博偉
    關(guān)鍵詞:張量滲流裂隙

    王俊奇,劉博偉,岳 瀟

    (1.華北電力大學(xué) 水利與水電工程學(xué)院,北京 102206;2.山東電力工程咨詢院有限公司,濟(jì)南 250014)

    裂隙巖體滲流是導(dǎo)致水利、土木工程事故的主要原因之一[1-2],還關(guān)系到采油、核廢料污染的防治、煤田瓦斯危害的防治和填埋垃圾的污水下滲等方面,工程中需要進(jìn)行細(xì)致的考察和分析。

    常用的裂隙巖體滲流分析數(shù)學(xué)模型有離散裂隙網(wǎng)絡(luò)模型、等效連續(xù)介質(zhì)模型、雙重介質(zhì)模型以及混合滲流模型[3]。因滲流主要在裂隙網(wǎng)絡(luò)內(nèi)發(fā)生,用離散裂隙網(wǎng)絡(luò)模型來模擬巖體滲流規(guī)律時精度較高,并且由于離散型裂隙網(wǎng)絡(luò)模型忽略了巖塊本身的孔隙系統(tǒng)及裂隙巖體之間的水交替過程,分析時更加簡便。二維離散裂隙網(wǎng)絡(luò)模型由于構(gòu)造簡單發(fā)展較為成熟,但計(jì)算結(jié)果和實(shí)際情況相比誤差較大。三維離散型裂隙網(wǎng)絡(luò)模型方面,李新強(qiáng)等[4]使用邊界元法計(jì)算裂隙網(wǎng)絡(luò)的滲透張量,并與有限元法流量分析結(jié)果比較驗(yàn)證方法的有效性。于青春等[5-7]基于管單元模型通過對裂隙網(wǎng)絡(luò)的滲透張量進(jìn)行求解來分析巖體的滲流情況,Ren等[8]提出用統(tǒng)一管網(wǎng)法對裂隙巖體內(nèi)的滲流情況進(jìn)行研究,Sun等[9]用3D統(tǒng)一管網(wǎng)法模擬含裂隙的多孔巖石中的裂隙情況。Wang等[6-10]從滲透張量的角度對裂隙網(wǎng)絡(luò)的REV尺寸進(jìn)行研究。王振偉等[11]通過巖體的滲透張量來確定滲流的REV尺寸,并研究了裂隙的密度和跡長對裂隙巖體的滲透張量和REV大小的影響。也有研究發(fā)現(xiàn)巖體裂隙存在優(yōu)勢水力路徑,Zhao等[12-13]發(fā)現(xiàn)在裂隙巖體中只有10%~20%的裂隙可以發(fā)生滲流,更多的裂隙僅起到存水或持水的作用。Zhang等[14-15]利用數(shù)值分析方法以及現(xiàn)場試驗(yàn)測試進(jìn)行滲流分析,對特定巖體中優(yōu)勢水流的形成機(jī)理及其路徑進(jìn)行了相關(guān)研究,林宏奕[16]利用質(zhì)點(diǎn)追蹤技術(shù)來對巖體裂隙中存在的優(yōu)勢水流路徑進(jìn)行分析,劉華梅等[17]通過對裂隙網(wǎng)絡(luò)滲流路徑進(jìn)行優(yōu)化,發(fā)現(xiàn)剔除對滲流無貢獻(xiàn)的裂隙,可以減少計(jì)算量和縮短計(jì)算時長,倪紹虎等[18-19]通過對裂隙巖體的滲透特性進(jìn)行研究發(fā)現(xiàn)在裂隙巖體中的滲流的確存在優(yōu)勢水力路徑。另外,也有一些學(xué)者采用不同的方法對砂巖油藏中優(yōu)勢滲流通道進(jìn)行識別和分析,如雷霆等[20]提出改進(jìn)的CRM模型結(jié)合生產(chǎn)動態(tài)數(shù)據(jù)來對優(yōu)勢滲透通道進(jìn)行識別描述。現(xiàn)有研究工作表明,確定滲透張量方面,一維管單元模型可以顯著縮減計(jì)算工作量,而且由于巖體中存在較大裂隙,溝通巖體范圍較大,一定程度也起到優(yōu)勢水力路徑作用,因此,在確定滲透張量時,一定程度上可以忽略掉連通范圍小的裂隙,但縮減的幅度沒有相關(guān)研究。

    本研究采用三維裂隙網(wǎng)絡(luò)模型對巖體裂隙滲流進(jìn)行分析,根據(jù)巖體露頭裂隙統(tǒng)計(jì)參數(shù),用蒙特卡羅法生成巖體的圓盤裂隙網(wǎng)絡(luò)樣本,將滲流通道簡化為管單元模型。一定程度上,巖體主干裂隙控制滲透性能,通過對裂隙規(guī)模的研究,探索不同縮減規(guī)模對滲透張量和巖體REV大小的影響,以及計(jì)算規(guī)模可縮減的幅度,并通過實(shí)際工程案例對此方法進(jìn)行了驗(yàn)證。

    1 三維裂隙網(wǎng)絡(luò)生成與滲透張量求解

    1.1 模型的建立與參數(shù)生成

    采用Baecher模型[21],建立離散網(wǎng)絡(luò)模型時采用三維整體坐標(biāo)系,將生成的所有裂隙通過該坐標(biāo)系確定。將坐標(biāo)系的主軸方向與地理位置一一對應(yīng),即X軸正半軸指向東,Y軸正半軸指向北,Z軸正半軸的方向指向上(與高程一致),按照右手法則來進(jìn)行方向的角度旋轉(zhuǎn)。生成裂隙所需要確定的幾何參數(shù)包括裂隙的方向(傾向α,傾角β)、大小(半徑r)、位置(圓盤中心點(diǎn)坐標(biāo))、隙寬b,即(α,β,r,x,y,z,b),每個參數(shù)變量的生成都遵循其自身的概率密度函數(shù),利用積分法產(chǎn)生隨機(jī)變量。

    裂隙的位置用裂隙面上有代表性的點(diǎn)來描述,在統(tǒng)計(jì)過程中用裂隙圓盤圓心來確定裂隙在空間中的位置,在三維空間中即圓心點(diǎn)(x,y,z)。假設(shè)圓心點(diǎn)在研究域各處出現(xiàn)的概率是相等的,即裂隙圓盤中心點(diǎn)的生成服從均勻分布。在生成裂隙網(wǎng)絡(luò)時,應(yīng)先確定生成域的大小,記錄圓心在該生成域內(nèi)的所有裂隙,然后在生成域內(nèi)再創(chuàng)建一定范圍內(nèi)的研究域,記錄圓心在該研究域內(nèi)的所有裂隙和與研究域邊界相交的裂隙。

    生成裂隙三維裂隙網(wǎng)絡(luò)后,求出各圓盤的交線,傳統(tǒng)方法是在此基礎(chǔ)上對裂隙圓盤剖分成平面單元進(jìn)行裂隙網(wǎng)絡(luò)水力學(xué)分析,但是一個研究域的裂隙多達(dá)上萬條,用平面單元分析會產(chǎn)生十幾萬的自由度,使求解效率低下。為簡化分析計(jì)算過程,Tsang[22]和Cacas等[23]提出水流在裂隙面內(nèi)以溝槽流形態(tài)運(yùn)動,因此,可以對裂隙滲流的分析計(jì)算進(jìn)行簡化,如圖1所示。張有天[1]認(rèn)為可以用一維管單元進(jìn)行三維裂隙網(wǎng)絡(luò)水力學(xué)分析。具體簡化模型為將相交裂隙節(jié)理圓盤的圓心和交線中心相互連接起來,然后生成管單元,只在管單元內(nèi)產(chǎn)生水流,圓心與交線的中點(diǎn)即為節(jié)點(diǎn)。對于滲流區(qū)域內(nèi)所有節(jié)點(diǎn)建立有限元方程組,用矩陣表示為

    圖1 三維裂隙網(wǎng)絡(luò)簡化分析概念簡圖Fig.1 Conceptual simplified analysis diagram of 3D fracture network

    KH=Q

    (1)

    式中:K為總滲透矩陣,在三維空間中有9個分量;H為節(jié)點(diǎn)水頭矩陣;Q為節(jié)點(diǎn)流量矩陣。由于是恒定流,方程可解。

    1.2 滲透張量的數(shù)值求法

    數(shù)值求法的基本原理為:在模型的某一相對面施加等水頭邊界,兩個面之間存在一定的水頭差,另外兩對相對面施加變水頭邊界[24](如圖2所示),滲流邊界的節(jié)點(diǎn)水頭分布情況如圖3所示。利用該模型共模擬了3組實(shí)驗(yàn),第一組實(shí)驗(yàn)中令其中一對相對面施加等水頭,在另外兩對相對面施加變水頭,進(jìn)行計(jì)算并保存相應(yīng)的計(jì)算結(jié)果。另外兩組實(shí)驗(yàn)只改變等水頭邊界施加的方向,其他不變。由于對稱性,每組實(shí)驗(yàn)的滲流量只計(jì)算6個即可,3組實(shí)驗(yàn)共產(chǎn)生的流量數(shù)據(jù)為18個。根據(jù)滲透張量的定義得[21]

    圖2 模型的邊界條件Fig.2 Boundary condition of the model

    圖3 裂隙網(wǎng)絡(luò)滲流水頭分布示意Fig.3 Schematic of seepage head distribution in fracture network

    (2)

    式中:kij為滲透張量在j方向滲透坡降下產(chǎn)生的i方向分量,qij為滲流量在j方向滲透坡降下產(chǎn)生的i方向分量,ΔPi為在i方向的水力坡降。

    對于均質(zhì)體,在各個方向上的入滲量和出滲量都相等。然而,對于節(jié)理巖體,由于巖體中裂隙的分布具有很強(qiáng)的隨機(jī)性,巖體的入滲流量和出滲流量并不總是相等。

    2 計(jì)算規(guī)??s減探索

    于青春等[25]研究發(fā)現(xiàn),滲流主要發(fā)生在較大的裂隙中,而較小的裂隙對滲流的影響十分微小,因此,在進(jìn)行裂隙滲流計(jì)算時可以在精度允許的情況下將較小的裂隙刪除掉,從而減少計(jì)算量,縮短計(jì)算時間。本文主要研究了在進(jìn)行滲流計(jì)算時可縮減的規(guī)模及其相應(yīng)的誤差,并用實(shí)際工程進(jìn)行了驗(yàn)證,當(dāng)然在實(shí)際工程設(shè)計(jì)中,應(yīng)根據(jù)工程的精度要求選擇合適的縮減量。在此約定:縮減規(guī)模為0.1即先將滲流裂隙網(wǎng)絡(luò)中的裂隙直徑按照由小到大的順序進(jìn)行排列,然后去掉排在前1/10直徑的裂隙,只用剩下的裂隙來進(jìn)行計(jì)算。

    規(guī)模的縮減一般為全排列縮減,即對裂隙網(wǎng)絡(luò)空間內(nèi)所有組數(shù)的裂隙從小到大進(jìn)行整體的規(guī)??s減。例如,若縮減規(guī)模為0.1,則將3組裂隙從小到大排列,刪除排在前10%的裂隙,保留剩下的90%進(jìn)行計(jì)算。全排列縮減方式的弊端在于容易較多地刪除直徑最小的那組裂隙,可能會降低滲透張量的準(zhǔn)確性??紤]到力學(xué)成因,此組裂隙可能開度較大,滲透能力并不弱,因此,提出另一種縮減方法,即對各組裂隙按指標(biāo)分別從小到大規(guī)??s減,例如,若縮減規(guī)模為0.1,分別對每組裂隙從小到大排列,再分別刪除每組前10%的裂隙,然后將各組的剩余裂隙重新統(tǒng)計(jì)在一起,進(jìn)行滲透張量計(jì)算。該方法在計(jì)算滲透張量前保證了各組裂隙縮減的平衡性,使各組裂隙數(shù)量均勻縮減。

    三維裂隙網(wǎng)絡(luò)的參數(shù)信息采用李新強(qiáng)等[4]研究的三維裂隙網(wǎng)絡(luò)滲流模型,如表1所示,在該區(qū)域共產(chǎn)生的裂隙組數(shù)為3。

    表1 三維裂隙網(wǎng)絡(luò)模型參數(shù)Tab.1 Model parameters of 3D fracture network

    由邊界元法計(jì)算確定初值,再用現(xiàn)場壓水試驗(yàn)校核得到滲透主系數(shù)(m/s)和滲透主方向?yàn)閇4]

    (3)

    應(yīng)用自編的C++有限元滲流程序,控制管徑與隙寬的比值為11.506 9[26],應(yīng)用上述數(shù)據(jù)進(jìn)行滲流模擬,選取100 m×100 m×100 m的立方體區(qū)域作為生成域,在生成域的中心選擇20 m×20 m×20 m的立方體作為研究域。對生成的裂隙進(jìn)行如上所述兩種方式的規(guī)??s減,表2為每組裂隙分別按相同比例的規(guī)??s減,表3為3組裂隙全排列情況下的規(guī)模縮減。表2和表3中滲透主值的誤差為

    (4)

    表2和表3中滲透主方向的誤差不可以直接求出,將計(jì)算所得滲透主值方向與現(xiàn)場壓水試驗(yàn)主值方向的夾角作為主方向的絕對誤差。兩滲透主值方向的夾角余弦為

    表2 每組裂隙分別按相同比例的規(guī)??s減及相應(yīng)誤差Tab.2 Scale reduction of each fracture group according to the same proportion and corresponding errors

    表3 3組裂隙全排列規(guī)模縮減及相應(yīng)誤差Tab.3 Scale reduction of three fracture groups with full arrangement and corresponding errors

    (5)

    將未縮減時所得到的滲透主系數(shù)和滲透主方向與式(3)進(jìn)行比較,相應(yīng)的誤差如下:

    (6)

    式中:εi為主值相對誤差,θi為主方向絕對誤差,i=1,2,3。

    由式(6)模擬可見,模擬計(jì)算的結(jié)果與現(xiàn)場壓水試驗(yàn)校正結(jié)果的誤差并不太大,在工程實(shí)踐可接受的范圍內(nèi)。

    為了便于比較兩種排列方式對滲透張量的影響,將表2和表3中數(shù)據(jù)分別繪制兩種方式的滲透主值和滲透主方向的誤差圖像,結(jié)果如圖4和圖5所示。

    圖4 兩種縮減方式滲透主值相對誤差對比Fig.4 Comparison of relative errors of principal permeability for two scale reduction methods

    圖5 兩種縮減方式滲透主方向絕對誤差對比Fig.5 Comparison of absolute errors of permeability directions for two scale reduction methods

    由圖4可知,兩種縮減方式滲透主值誤差在縮減規(guī)模0.1~0.7基本一致。當(dāng)縮減規(guī)模達(dá)到0.8時,分組排列的誤差絕對值為11%,全排列的誤差絕對值為9.7%。當(dāng)縮減規(guī)模達(dá)到0.8以后,誤差絕對值開始驟然增大,分組排列的誤差絕對值為32.6%,全排列的誤差絕對值為41.3%,此時縮減效果已成病態(tài),不適宜再縮減裂隙規(guī)模進(jìn)行計(jì)算。同樣根據(jù)圖5可知,兩種縮減方式滲透主方向均值在縮減規(guī)模0.1~0.7波動較小。當(dāng)縮減規(guī)模達(dá)到0.7以后,分組排列的角度誤差絕對值為2.03°,全排列的誤差絕對值為2.60°。當(dāng)縮減規(guī)模達(dá)到0.8以后,誤差絕對值增幅增大,此時分組排列的誤差絕對值為9°,全排列的誤差絕對值為4.9°,此時方向誤差已經(jīng)不夠準(zhǔn)確,不適宜再縮減裂隙規(guī)模進(jìn)行計(jì)算。

    綜上,對于兩種縮減方式效果基本一致,實(shí)際工程中可以采用任意一種縮減方式。當(dāng)最大縮減規(guī)模為0.7時,最為適宜本例,當(dāng)縮減規(guī)模大于0.7時,裂隙滲流網(wǎng)絡(luò)中重要的裂隙將被刪除,導(dǎo)致滲流分析的結(jié)果與裂隙巖體的實(shí)際滲流情況差別較大。

    3 巖體表征單元體(REV)的確定

    由于在天然巖體中存在的裂隙是隨機(jī)分布的,國內(nèi)外學(xué)者在進(jìn)行巖體力學(xué)參數(shù)選取時經(jīng)常會對巖體的表征單元體(REV)進(jìn)行分析。巖體的各項(xiàng)等效參數(shù)具有十分明顯的尺寸效應(yīng),當(dāng)巖體的試樣體積增加到一定值V0時,各項(xiàng)滲透參數(shù)趨于穩(wěn)定,這一V0稱為巖體滲透性代表單元體積(REV),確定巖體REV的尺寸是進(jìn)行裂隙滲流研究的一個重要方面。目前,國內(nèi)外學(xué)者對于REV尺寸有從等效力學(xué)參數(shù)進(jìn)行分析確定的,也有從水力學(xué)參數(shù)方面對REV進(jìn)行研究的,還有從裂隙塊體幾何參數(shù)對巖體的REV進(jìn)行研究的?;诓煌瑤r體參數(shù)的尺寸效應(yīng)而得到的表征單元體的尺寸都不太相同,周創(chuàng)兵等[27]提出了3種確定巖體REV的方法,即數(shù)值模擬法、地質(zhì)統(tǒng)計(jì)法和能量疊加法。

    本研究運(yùn)用數(shù)值模擬的方法來計(jì)算巖體的滲透張量,根據(jù)滲透張量3個滲透主系數(shù)的變化情況來確定巖體的REV尺寸,通過對不同縮減規(guī)模情況下裂隙網(wǎng)絡(luò)REV的研究,進(jìn)一步探索規(guī)??s減對裂隙巖體REV的影響。在縮減規(guī)模確定的情況下,對研究域從9 m×9 m×9 m開始逐漸增加到23 m×23 m×23 m,每次研究域的邊長增加1 m。以計(jì)算規(guī)模未進(jìn)行縮減時的情況為例,管徑與隙寬的比值為11.506 9,控制生成域的大小為100 m×100 m×100 m,通過不斷改變研究域的大小,得到不同研究域的滲透張量,將所得的裂隙網(wǎng)絡(luò)滲透主系數(shù)繪制成折線圖,如圖6所示??梢钥闯?,當(dāng)所生成的裂隙網(wǎng)絡(luò)研究域的邊長增大到15 m后,隨著研究域邊長的增大,巖體滲透張量的滲透主系數(shù)變化情況基本趨于穩(wěn)定,故當(dāng)縮減規(guī)模為0時的REV尺寸可取15 m×15 m×15 m。同理,對其他縮減規(guī)模下的REV尺寸進(jìn)行分析,并將得到的REV的邊長繪制在折線圖上,如圖7所示??梢钥闯觯S著裂隙網(wǎng)絡(luò)縮減規(guī)模數(shù)值的增大,巖體的REV尺寸逐漸變大,說明裂隙密度對巖體的表征單元體大小有一定的影響,REV尺寸隨裂隙密度減小而增大,該認(rèn)識與李錦輝等[28]對裂隙數(shù)量與REV尺寸關(guān)系[11]的研究結(jié)論相一致。

    圖6 計(jì)算規(guī)模不縮減時的滲透主系數(shù)Fig.6 Principal permeability curves without calculation scale reduction

    圖7 不同縮減規(guī)模下的REV邊長Fig.7 REV side lengths for different scale reductions

    4 工程驗(yàn)證

    中國水利水電科學(xué)研究院的科研人員在新疆對某水利樞紐的左岸邊坡巖體進(jìn)行現(xiàn)場調(diào)研后,根據(jù)巖體中結(jié)構(gòu)面的分布和組合情況得到如表4所示的統(tǒng)計(jì)結(jié)果[29]。

    表4 巖體結(jié)構(gòu)面幾何參數(shù)統(tǒng)計(jì)結(jié)果Tab.4 Statistical results of geometric parameters of rock mass structural surface

    對此工程實(shí)例,用自編的程序進(jìn)行數(shù)值模擬,將生成域的大小設(shè)為40 m×40 m×40 m,研究域的大小為15 m×15 m×15 m。不縮減時在生成的研究域內(nèi)裂隙共有8 129條,在裂隙網(wǎng)絡(luò)中節(jié)點(diǎn)數(shù)為53 650,單元數(shù)為99 659,經(jīng)過數(shù)值仿真模擬得到相應(yīng)的滲透主系數(shù)及滲透主方向?yàn)?/p>

    (7)

    本工程實(shí)例采用縮減規(guī)模為0.7,經(jīng)縮減后共剩2 394條裂隙,在裂隙網(wǎng)絡(luò)中共有15 986個節(jié)點(diǎn)、29 716個單元,經(jīng)數(shù)值仿真模擬得到相對應(yīng)的滲透主系數(shù)及滲透主方向?yàn)?/p>

    (8)

    為形象化滲透特性,根據(jù)滲透橢球方程

    (9)

    將式(7)、(8)用MATLAB軟件繪制縮減前后的滲透橢球圖,如圖8所示,由于橢球具有對稱性,為更好地顯示,取橢球的一半進(jìn)行觀察。

    圖8 縮減前后滲透橢球?qū)Ρ菷ig.8 Comparison of permeability ellipsoid before and after reduction

    由式(7)和(8)對比及圖8可以看出,縮減前與縮減規(guī)模為0.7的滲透張量及滲透橢球相差不大,說明采用縮減規(guī)模為0.7時,計(jì)算出的滲透主系數(shù)和滲透主方向與未縮減前計(jì)算的結(jié)果比較接近,這一方法對于工程實(shí)際具有一定的參考價(jià)值。

    對縮減前后裂隙巖體模型的REV尺寸進(jìn)行計(jì)算得:

    1)未縮減時裂隙巖體的REV邊長為4 m,相應(yīng)的滲透主系數(shù)和滲透主方向?yàn)?/p>

    (10)

    2)縮減規(guī)模為0.7時的REV邊長為6 m,相應(yīng)的滲透主系數(shù)及滲透主方向?yàn)?/p>

    (11)

    將式(11)與(10)進(jìn)行比較,滲透主系數(shù)及滲透主方向的誤差分別為

    (12)

    式中:εi為主值相對誤差,θi為主方向絕對誤差,i=1,2,3。

    由式(12)可知,未縮減時邊長為4 m的REV與縮減為0.7時邊長為6 m的REV滲透主系數(shù)相對誤差的最大值為0.257,滲透主方向絕對誤差的最大值為4.62°,縮減前后REV的滲透主系數(shù)和滲透主方向誤差均較小。

    5 結(jié) 論

    1)通過刪減次要裂隙、保留主干裂隙進(jìn)行滲透張量計(jì)算,結(jié)果表明,裂隙網(wǎng)絡(luò)確實(shí)存在優(yōu)勢路徑,保留主干裂隙、縮減計(jì)算規(guī)模0.7時,計(jì)算得到的滲透張量仍然具有較好的精度,可以滿足工程需要。

    2)對全排列和分組排列兩種方法的規(guī)??s減程度與結(jié)果精度進(jìn)行分析,裂隙的滲透主值與主方向與現(xiàn)場壓水試驗(yàn)校核后的結(jié)果誤差在可以接受的范圍內(nèi),而且兩種方法的縮減效果基本一致。

    3)裂隙巖體REV尺寸的大小會隨著縮減規(guī)模的增大而變大,通過對巖體表征單元體的討論,發(fā)現(xiàn)縮減規(guī)模從0變化到0.7時,裂隙巖體的REV體積由15 m×15 m×15 m逐漸增大到21 m×21 m×21 m,說明REV的尺寸受縮減規(guī)模的影響。工程實(shí)例REV的邊長隨著縮減規(guī)模的增加從4 m增大到了6 m。

    4)通過實(shí)例驗(yàn)證縮減計(jì)算的實(shí)用性。巖體裂隙網(wǎng)絡(luò)滲流滿足逾滲條件,存在REV情況下,對于那些對精度要求不高、地質(zhì)條件不復(fù)雜的情況,滲透張量可以采用縮減規(guī)模數(shù)為0.7進(jìn)行計(jì)算。對于工程精度要求較高、地質(zhì)條件復(fù)雜的重要工程項(xiàng)目,應(yīng)結(jié)合現(xiàn)場的實(shí)測資料與試驗(yàn),選擇合適的縮減規(guī)模數(shù)進(jìn)行分析。在誤差允許的條件下相對地減少最大縮減量,可以達(dá)到減少計(jì)算量的目的,提高三維裂隙離散網(wǎng)絡(luò)模型在工程中的實(shí)用性。

    猜你喜歡
    張量滲流裂隙
    裂隙腦室綜合征的診斷治療新進(jìn)展
    偶數(shù)階張量core逆的性質(zhì)和應(yīng)用
    四元數(shù)張量方程A*NX=B 的通解
    裂隙燈檢查的個性化應(yīng)用(下)
    擴(kuò)散張量成像MRI 在CO中毒后遲發(fā)腦病中的應(yīng)用
    《老炮兒》:在時代裂隙中揚(yáng)棄焦慮
    電影新作(2016年1期)2016-02-27 09:16:48
    工程中張量概念的思考
    河南科技(2014年19期)2014-02-27 14:15:33
    簡述滲流作用引起的土體破壞及防治措施
    河南科技(2014年12期)2014-02-27 14:10:26
    關(guān)于渠道滲流計(jì)算方法的選用
    河南科技(2014年11期)2014-02-27 14:09:48
    非完全充填裂隙的應(yīng)力應(yīng)變特性研究
    河南科技(2014年11期)2014-02-27 14:09:42
    精品久久久久久久末码| 国产亚洲欧美98| 亚洲av成人av| 亚洲 欧美 日韩 在线 免费| 欧美精品亚洲一区二区| 欧美色欧美亚洲另类二区| 热99re8久久精品国产| 亚洲五月天丁香| 亚洲成人国产一区在线观看| 少妇粗大呻吟视频| 国产免费男女视频| 久久人妻av系列| 久久午夜亚洲精品久久| 99久久无色码亚洲精品果冻| 日本a在线网址| 一个人免费在线观看的高清视频| 在线观看一区二区三区| 欧美在线黄色| 午夜免费成人在线视频| 成人一区二区视频在线观看| 一级毛片女人18水好多| 国产aⅴ精品一区二区三区波| 色尼玛亚洲综合影院| 亚洲第一电影网av| 国产精品永久免费网站| 法律面前人人平等表现在哪些方面| 在线看三级毛片| bbb黄色大片| 精品一区二区三区四区五区乱码| e午夜精品久久久久久久| 一级作爱视频免费观看| 欧美人与性动交α欧美精品济南到| tocl精华| 最近最新免费中文字幕在线| 精品国产乱码久久久久久男人| 我要搜黄色片| 欧美性猛交黑人性爽| 一本精品99久久精品77| 亚洲国产精品999在线| 国产av又大| 国产在线观看jvid| 中出人妻视频一区二区| 后天国语完整版免费观看| 日韩欧美免费精品| 国产黄a三级三级三级人| 巨乳人妻的诱惑在线观看| 日韩中文字幕欧美一区二区| 看片在线看免费视频| 日韩有码中文字幕| 亚洲精品粉嫩美女一区| 国产主播在线观看一区二区| 成人av一区二区三区在线看| 亚洲精品国产一区二区精华液| 久久这里只有精品19| 日韩 欧美 亚洲 中文字幕| 亚洲一区二区三区不卡视频| 久久 成人 亚洲| 欧美乱妇无乱码| 人成视频在线观看免费观看| 国产精品,欧美在线| 亚洲无线在线观看| 久久精品国产清高在天天线| 美女午夜性视频免费| 婷婷精品国产亚洲av在线| √禁漫天堂资源中文www| 久久婷婷成人综合色麻豆| 宅男免费午夜| 大型黄色视频在线免费观看| 国产亚洲av高清不卡| 久久久久免费精品人妻一区二区| 制服诱惑二区| 国产精品久久视频播放| 亚洲精品色激情综合| 2021天堂中文幕一二区在线观| 亚洲电影在线观看av| 色老头精品视频在线观看| 一级毛片女人18水好多| 人妻丰满熟妇av一区二区三区| 国产精品亚洲av一区麻豆| 国产一区二区三区在线臀色熟女| 12—13女人毛片做爰片一| 又紧又爽又黄一区二区| 亚洲最大成人中文| 老司机深夜福利视频在线观看| 亚洲在线自拍视频| 久久久久精品国产欧美久久久| 欧美色欧美亚洲另类二区| 99精品久久久久人妻精品| 给我免费播放毛片高清在线观看| 国产一区二区在线av高清观看| videosex国产| 嫩草影院精品99| 叶爱在线成人免费视频播放| 久久久国产精品麻豆| АⅤ资源中文在线天堂| 亚洲国产欧洲综合997久久,| 久久久久久久午夜电影| 久久香蕉国产精品| 亚洲熟妇熟女久久| av在线天堂中文字幕| 黄色视频,在线免费观看| 长腿黑丝高跟| 亚洲人成77777在线视频| 国产欧美日韩精品亚洲av| 深夜精品福利| 中文字幕人妻丝袜一区二区| 欧美日韩中文字幕国产精品一区二区三区| 亚洲九九香蕉| 母亲3免费完整高清在线观看| 91国产中文字幕| 无人区码免费观看不卡| 久久久久性生活片| 国产v大片淫在线免费观看| 成人三级做爰电影| 最近最新中文字幕大全电影3| 窝窝影院91人妻| 在线观看免费视频日本深夜| 免费看a级黄色片| 99精品欧美一区二区三区四区| 成在线人永久免费视频| 啪啪无遮挡十八禁网站| 欧美日韩一级在线毛片| 亚洲精品国产一区二区精华液| 国产精品一区二区免费欧美| 亚洲中文日韩欧美视频| 欧美乱码精品一区二区三区| 成人国产综合亚洲| 又紧又爽又黄一区二区| 亚洲专区中文字幕在线| 18禁黄网站禁片免费观看直播| 欧美黑人精品巨大| 亚洲精品国产精品久久久不卡| 欧美日韩乱码在线| 日韩高清综合在线| 国内精品久久久久久久电影| 欧美中文综合在线视频| 亚洲 欧美一区二区三区| 欧美日本视频| 国产精品免费视频内射| 99国产极品粉嫩在线观看| 国产精品av视频在线免费观看| 女生性感内裤真人,穿戴方法视频| 国产亚洲精品久久久久5区| 国产成人精品久久二区二区91| 午夜福利视频1000在线观看| 无遮挡黄片免费观看| 久99久视频精品免费| 可以免费在线观看a视频的电影网站| 淫妇啪啪啪对白视频| 久久精品亚洲精品国产色婷小说| 婷婷亚洲欧美| 国产高清有码在线观看视频 | 国内揄拍国产精品人妻在线| 嫩草影视91久久| 麻豆久久精品国产亚洲av| 99国产极品粉嫩在线观看| 国内揄拍国产精品人妻在线| 国产午夜精品久久久久久| 欧美性猛交╳xxx乱大交人| 狠狠狠狠99中文字幕| 久久午夜亚洲精品久久| 亚洲第一欧美日韩一区二区三区| 高潮久久久久久久久久久不卡| 国产亚洲精品一区二区www| 国产亚洲精品久久久久久毛片| 日本免费一区二区三区高清不卡| 亚洲美女黄片视频| 国产私拍福利视频在线观看| 国产精品精品国产色婷婷| 久久精品成人免费网站| 国产成人av教育| 国产精品永久免费网站| 啦啦啦观看免费观看视频高清| 国产精品亚洲美女久久久| 亚洲成人免费电影在线观看| 午夜免费成人在线视频| 免费电影在线观看免费观看| 露出奶头的视频| 午夜福利视频1000在线观看| 女人被狂操c到高潮| 日本一本二区三区精品| 99热这里只有是精品50| 国产乱人伦免费视频| 嫩草影视91久久| 日韩精品免费视频一区二区三区| 97人妻精品一区二区三区麻豆| 国产av一区二区精品久久| a级毛片a级免费在线| 99riav亚洲国产免费| 女人被狂操c到高潮| 免费观看精品视频网站| 黄色成人免费大全| 亚洲乱码一区二区免费版| 久久精品国产亚洲av香蕉五月| 麻豆国产97在线/欧美 | 99在线视频只有这里精品首页| 久久草成人影院| 亚洲av中文字字幕乱码综合| 亚洲国产高清在线一区二区三| 亚洲美女黄片视频| 亚洲成人国产一区在线观看| 亚洲成人精品中文字幕电影| 国产真实乱freesex| 精品日产1卡2卡| 男男h啪啪无遮挡| 国产精品日韩av在线免费观看| 欧美黑人精品巨大| 免费看美女性在线毛片视频| 亚洲aⅴ乱码一区二区在线播放 | 日日夜夜操网爽| 亚洲色图av天堂| 日韩精品青青久久久久久| 国产精品免费一区二区三区在线| 亚洲人成77777在线视频| 成年人黄色毛片网站| 最新在线观看一区二区三区| 国内精品久久久久精免费| 欧美不卡视频在线免费观看 | 欧美日韩亚洲综合一区二区三区_| 十八禁人妻一区二区| 91麻豆av在线| 免费一级毛片在线播放高清视频| 欧美激情久久久久久爽电影| 亚洲av成人不卡在线观看播放网| 美女高潮喷水抽搐中文字幕| 51午夜福利影视在线观看| 久99久视频精品免费| 日本一本二区三区精品| 欧美不卡视频在线免费观看 | 麻豆国产av国片精品| 国产精品永久免费网站| 亚洲av成人一区二区三| 999久久久精品免费观看国产| 美女午夜性视频免费| 日韩欧美免费精品| 一本一本综合久久| aaaaa片日本免费| 高潮久久久久久久久久久不卡| 色综合欧美亚洲国产小说| 久久 成人 亚洲| 妹子高潮喷水视频| 在线永久观看黄色视频| а√天堂www在线а√下载| aaaaa片日本免费| 啦啦啦韩国在线观看视频| 国产精品98久久久久久宅男小说| 久久国产乱子伦精品免费另类| 91成年电影在线观看| 在线观看免费日韩欧美大片| 中出人妻视频一区二区| 亚洲av片天天在线观看| 麻豆久久精品国产亚洲av| 美女黄网站色视频| 久久亚洲真实| 成年女人毛片免费观看观看9| 在线观看舔阴道视频| 国产三级中文精品| 亚洲av美国av| 桃红色精品国产亚洲av| 法律面前人人平等表现在哪些方面| 两性午夜刺激爽爽歪歪视频在线观看 | 啦啦啦观看免费观看视频高清| 亚洲va日本ⅴa欧美va伊人久久| 99久久久亚洲精品蜜臀av| 欧美黑人欧美精品刺激| 男人的好看免费观看在线视频 | 男女下面进入的视频免费午夜| 99国产综合亚洲精品| 国产熟女午夜一区二区三区| 国产高清videossex| 亚洲欧美精品综合一区二区三区| 国产一区在线观看成人免费| 岛国在线免费视频观看| 亚洲一区高清亚洲精品| 人人妻人人看人人澡| 免费在线观看影片大全网站| 国产精品一及| 男女之事视频高清在线观看| 又爽又黄无遮挡网站| 亚洲va日本ⅴa欧美va伊人久久| 日本五十路高清| 久久精品国产清高在天天线| 69av精品久久久久久| 国产单亲对白刺激| 国产熟女xx| 全区人妻精品视频| 色噜噜av男人的天堂激情| 国产成人精品久久二区二区91| 色综合欧美亚洲国产小说| 婷婷亚洲欧美| 午夜a级毛片| 一个人观看的视频www高清免费观看 | 一进一出抽搐gif免费好疼| 香蕉av资源在线| 高清毛片免费观看视频网站| 欧美久久黑人一区二区| 精品不卡国产一区二区三区| 久久精品aⅴ一区二区三区四区| 亚洲精华国产精华精| 成人18禁在线播放| 国产精品国产高清国产av| 一a级毛片在线观看| 岛国在线免费视频观看| 校园春色视频在线观看| 欧美精品亚洲一区二区| 他把我摸到了高潮在线观看| 中国美女看黄片| 亚洲在线自拍视频| 12—13女人毛片做爰片一| 国产精品一区二区三区四区久久| 午夜精品久久久久久毛片777| 黄色毛片三级朝国网站| 男人舔奶头视频| 99国产综合亚洲精品| 成人av在线播放网站| 色综合婷婷激情| 国产av在哪里看| 欧美日韩亚洲综合一区二区三区_| 欧美色欧美亚洲另类二区| 亚洲国产欧美网| 国产三级中文精品| 叶爱在线成人免费视频播放| 男人舔奶头视频| 男女下面进入的视频免费午夜| 69av精品久久久久久| 桃色一区二区三区在线观看| 欧美绝顶高潮抽搐喷水| 欧美又色又爽又黄视频| 久久久精品欧美日韩精品| 久久久国产成人精品二区| 日韩精品中文字幕看吧| 国内揄拍国产精品人妻在线| 精品久久久久久久末码| 两个人看的免费小视频| 每晚都被弄得嗷嗷叫到高潮| 久久香蕉国产精品| 久久午夜亚洲精品久久| 看片在线看免费视频| 午夜免费成人在线视频| 久久精品国产亚洲av高清一级| 丰满人妻熟妇乱又伦精品不卡| 在线a可以看的网站| 亚洲av美国av| 午夜精品久久久久久毛片777| 男男h啪啪无遮挡| 国内少妇人妻偷人精品xxx网站 | 丰满的人妻完整版| 久久九九热精品免费| 99久久99久久久精品蜜桃| 精品久久久久久,| а√天堂www在线а√下载| 在线播放国产精品三级| 精品国产超薄肉色丝袜足j| 免费高清视频大片| 制服诱惑二区| 亚洲 欧美一区二区三区| 午夜免费成人在线视频| 99国产综合亚洲精品| e午夜精品久久久久久久| 欧美一级毛片孕妇| 亚洲va日本ⅴa欧美va伊人久久| 国产不卡一卡二| 成人手机av| 日日摸夜夜添夜夜添小说| 村上凉子中文字幕在线| 51午夜福利影视在线观看| www.熟女人妻精品国产| 妹子高潮喷水视频| 中文字幕熟女人妻在线| 免费看日本二区| 欧美精品亚洲一区二区| 2021天堂中文幕一二区在线观| av有码第一页| 最新在线观看一区二区三区| 欧美不卡视频在线免费观看 | 在线观看美女被高潮喷水网站 | 午夜福利成人在线免费观看| 黄色毛片三级朝国网站| 亚洲人成网站在线播放欧美日韩| 久久久水蜜桃国产精品网| 女警被强在线播放| 亚洲一区二区三区不卡视频| 日韩精品免费视频一区二区三区| 母亲3免费完整高清在线观看| 91大片在线观看| 一区二区三区高清视频在线| 天天添夜夜摸| 三级毛片av免费| 男男h啪啪无遮挡| 亚洲成人中文字幕在线播放| 高清毛片免费观看视频网站| 亚洲国产日韩欧美精品在线观看 | 国产区一区二久久| 婷婷亚洲欧美| 在线视频色国产色| 99久久无色码亚洲精品果冻| 伦理电影免费视频| 久久久久久九九精品二区国产 | 人妻丰满熟妇av一区二区三区| 国产三级中文精品| 欧美在线一区亚洲| 啦啦啦观看免费观看视频高清| 亚洲欧美日韩无卡精品| 欧美成人一区二区免费高清观看 | 国产精品久久电影中文字幕| 老熟妇乱子伦视频在线观看| 亚洲男人天堂网一区| 国模一区二区三区四区视频 | 国产精品永久免费网站| 99久久久亚洲精品蜜臀av| 黄色视频,在线免费观看| 法律面前人人平等表现在哪些方面| 国产乱人伦免费视频| 日韩三级视频一区二区三区| 在线看三级毛片| 精品久久久久久成人av| 亚洲一区二区三区不卡视频| 99re在线观看精品视频| 色精品久久人妻99蜜桃| 91老司机精品| 日韩欧美精品v在线| 亚洲av熟女| 99热这里只有是精品50| 在线观看www视频免费| 欧美日韩乱码在线| 观看免费一级毛片| 午夜免费成人在线视频| 久久久久九九精品影院| 国产成+人综合+亚洲专区| 亚洲av美国av| 精品无人区乱码1区二区| a在线观看视频网站| 亚洲国产精品sss在线观看| 日本 欧美在线| 亚洲午夜理论影院| 一a级毛片在线观看| 无人区码免费观看不卡| 麻豆成人av在线观看| 国产精品爽爽va在线观看网站| 中文在线观看免费www的网站 | 别揉我奶头~嗯~啊~动态视频| 激情在线观看视频在线高清| 欧美日韩中文字幕国产精品一区二区三区| 国产精品永久免费网站| 国内少妇人妻偷人精品xxx网站 | 国产精品久久视频播放| 中文字幕高清在线视频| 免费av毛片视频| 特大巨黑吊av在线直播| 国产精品综合久久久久久久免费| 成人午夜高清在线视频| 黑人操中国人逼视频| www.999成人在线观看| 99热这里只有精品一区 | 欧美绝顶高潮抽搐喷水| 亚洲熟妇熟女久久| 狂野欧美激情性xxxx| 老熟妇仑乱视频hdxx| 最近最新免费中文字幕在线| 亚洲专区字幕在线| 久久久久久久午夜电影| www.熟女人妻精品国产| 巨乳人妻的诱惑在线观看| 97超级碰碰碰精品色视频在线观看| 三级国产精品欧美在线观看 | 亚洲五月天丁香| 老司机深夜福利视频在线观看| 久久久久久久午夜电影| xxx96com| 国产不卡一卡二| 成人国语在线视频| 人妻夜夜爽99麻豆av| 这个男人来自地球电影免费观看| 久久久久国产精品人妻aⅴ院| 中文字幕熟女人妻在线| 精品国产乱子伦一区二区三区| 国产免费av片在线观看野外av| 久久草成人影院| 婷婷精品国产亚洲av| 亚洲av第一区精品v没综合| 一区二区三区国产精品乱码| 亚洲男人天堂网一区| 欧美中文综合在线视频| 国产av在哪里看| 狂野欧美白嫩少妇大欣赏| 亚洲 欧美 日韩 在线 免费| 日本 欧美在线| 天天一区二区日本电影三级| 久久这里只有精品中国| 亚洲第一欧美日韩一区二区三区| 女人高潮潮喷娇喘18禁视频| 操出白浆在线播放| 日本黄色视频三级网站网址| 亚洲精华国产精华精| 国产精品国产高清国产av| 亚洲成人国产一区在线观看| 国产精品久久电影中文字幕| 日韩欧美精品v在线| 久久久精品国产亚洲av高清涩受| 麻豆久久精品国产亚洲av| 国产成人影院久久av| av福利片在线观看| 久久久久久久精品吃奶| 国产精品亚洲美女久久久| 琪琪午夜伦伦电影理论片6080| 久久久久精品国产欧美久久久| avwww免费| 午夜精品久久久久久毛片777| 午夜视频精品福利| 1024手机看黄色片| 妹子高潮喷水视频| 成人手机av| 久久久久久免费高清国产稀缺| 欧美在线一区亚洲| 好看av亚洲va欧美ⅴa在| 久久久久久久午夜电影| 日日夜夜操网爽| 国产精品一区二区三区四区久久| 最近最新中文字幕大全电影3| 男女下面进入的视频免费午夜| 免费看日本二区| 久久人人精品亚洲av| 久久这里只有精品19| 成人av一区二区三区在线看| 黄片大片在线免费观看| 在线a可以看的网站| 1024视频免费在线观看| 又紧又爽又黄一区二区| 亚洲精品久久国产高清桃花| 又紧又爽又黄一区二区| 丝袜人妻中文字幕| 三级男女做爰猛烈吃奶摸视频| 午夜福利免费观看在线| 性欧美人与动物交配| 国产精品99久久99久久久不卡| 天堂√8在线中文| 久久精品国产99精品国产亚洲性色| 最新美女视频免费是黄的| 成年版毛片免费区| 亚洲精品色激情综合| www.www免费av| 亚洲avbb在线观看| 久久精品亚洲精品国产色婷小说| 亚洲成a人片在线一区二区| 亚洲va日本ⅴa欧美va伊人久久| 黄色视频,在线免费观看| 亚洲午夜理论影院| 床上黄色一级片| 免费在线观看日本一区| 美女黄网站色视频| 两个人看的免费小视频| 黄片小视频在线播放| 99热这里只有是精品50| 国产高清视频在线观看网站| 久久国产精品人妻蜜桃| 黄片小视频在线播放| 午夜精品一区二区三区免费看| 亚洲精品色激情综合| 91老司机精品| 日韩精品免费视频一区二区三区| www.精华液| 黄色女人牲交| 91大片在线观看| 岛国视频午夜一区免费看| 日本一本二区三区精品| 国产欧美日韩一区二区三| 久久久久国产一级毛片高清牌| 婷婷精品国产亚洲av在线| 最新在线观看一区二区三区| 亚洲国产精品999在线| 在线看三级毛片| 亚洲成人精品中文字幕电影| 欧美中文日本在线观看视频| 成人三级黄色视频| 欧美中文日本在线观看视频| 99re在线观看精品视频| 欧美国产日韩亚洲一区| 一进一出好大好爽视频| 91大片在线观看| 波多野结衣高清作品| 国产激情欧美一区二区| 九九热线精品视视频播放| 国产精品亚洲美女久久久| 国产高清有码在线观看视频 | 欧美最黄视频在线播放免费| 国产精品香港三级国产av潘金莲| 国产成人av教育| 久久香蕉精品热| 我要搜黄色片| 啦啦啦观看免费观看视频高清| 国产主播在线观看一区二区| 国内精品久久久久久久电影| 黑人巨大精品欧美一区二区mp4| 国产精品九九99| 熟女少妇亚洲综合色aaa.| 男人舔女人的私密视频| 亚洲精品久久国产高清桃花| 黑人操中国人逼视频| 欧美黄色片欧美黄色片| 又爽又黄无遮挡网站| 又粗又爽又猛毛片免费看| 亚洲五月婷婷丁香| 亚洲av五月六月丁香网| 精品久久久久久久毛片微露脸| 夜夜躁狠狠躁天天躁| 国产高清激情床上av| 99久久精品国产亚洲精品| 久久婷婷成人综合色麻豆| 久久久久国内视频| 天天一区二区日本电影三级| 九九热线精品视视频播放| 丰满人妻熟妇乱又伦精品不卡| 成人三级黄色视频| 中文字幕精品亚洲无线码一区| 日本成人三级电影网站|