• <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
    狂野欧美白嫩少妇大欣赏| 久久久久久国产a免费观看| 国产亚洲精品久久久久久毛片| 国产成人啪精品午夜网站| 一夜夜www| 亚洲精品成人久久久久久| 亚洲国产精品合色在线| 欧美成人免费av一区二区三区| 亚洲欧美一区二区三区黑人| 欧美在线黄色| 国产精品98久久久久久宅男小说| 国产久久久一区二区三区| 国产成人影院久久av| 深夜精品福利| 老熟妇仑乱视频hdxx| 国产精品野战在线观看| 国产麻豆成人av免费视频| 国产伦人伦偷精品视频| 亚洲精品色激情综合| 我的老师免费观看完整版| 99久久精品一区二区三区| 免费在线观看亚洲国产| 日韩欧美免费精品| 欧美bdsm另类| 欧美区成人在线视频| 高潮久久久久久久久久久不卡| 三级毛片av免费| 日韩欧美免费精品| 久久婷婷人人爽人人干人人爱| 亚洲人成网站在线播放欧美日韩| 免费观看精品视频网站| 99热6这里只有精品| 一夜夜www| 久久中文看片网| 老熟妇仑乱视频hdxx| 搞女人的毛片| 亚洲第一电影网av| 午夜精品在线福利| 淫秽高清视频在线观看| 女人高潮潮喷娇喘18禁视频| www日本在线高清视频| 国产高清videossex| 国产淫片久久久久久久久 | 天堂av国产一区二区熟女人妻| 99久久精品国产亚洲精品| 欧美成人性av电影在线观看| 国产精品一及| 大型黄色视频在线免费观看| 搞女人的毛片| 九色国产91popny在线| 亚洲成人久久爱视频| 天美传媒精品一区二区| 久久中文看片网| 人人妻,人人澡人人爽秒播| 欧美日韩瑟瑟在线播放| 国产高清三级在线| 国产日本99.免费观看| 久久久久久久精品吃奶| 中文字幕人妻丝袜一区二区| 女人高潮潮喷娇喘18禁视频| 国语自产精品视频在线第100页| 久久久久久久亚洲中文字幕 | 欧美色视频一区免费| 叶爱在线成人免费视频播放| 久久久久久久久大av| 亚洲 国产 在线| 99久久九九国产精品国产免费| 免费观看精品视频网站| 婷婷丁香在线五月| 亚洲成av人片免费观看| 久久久精品欧美日韩精品| 日本 av在线| 欧美一区二区精品小视频在线| 两性午夜刺激爽爽歪歪视频在线观看| av欧美777| 日本黄色视频三级网站网址| 在线国产一区二区在线| 亚洲不卡免费看| 国产成人aa在线观看| 两个人视频免费观看高清| 亚洲18禁久久av| 国产午夜精品久久久久久一区二区三区 | 少妇人妻一区二区三区视频| 在线天堂最新版资源| 国产亚洲精品一区二区www| 午夜久久久久精精品| 亚洲人成伊人成综合网2020| 嫩草影院入口| 琪琪午夜伦伦电影理论片6080| 国产不卡一卡二| 一区二区三区高清视频在线| 网址你懂的国产日韩在线| 成人一区二区视频在线观看| 精品日产1卡2卡| 好男人在线观看高清免费视频| 国产伦人伦偷精品视频| 亚洲av免费在线观看| 国产av在哪里看| 亚洲国产高清在线一区二区三| 国产亚洲精品一区二区www| 亚洲人成网站高清观看| 首页视频小说图片口味搜索| 精品欧美国产一区二区三| 性色av乱码一区二区三区2| 亚洲国产欧洲综合997久久,| 午夜福利欧美成人| 国产乱人视频| 日本黄色视频三级网站网址| www日本在线高清视频| 色综合婷婷激情| 日本一本二区三区精品| 欧美另类亚洲清纯唯美| x7x7x7水蜜桃| 久久性视频一级片| 91在线精品国自产拍蜜月 | 制服人妻中文乱码| 亚洲一区高清亚洲精品| av女优亚洲男人天堂| 亚洲五月天丁香| 亚洲七黄色美女视频| 成年女人永久免费观看视频| 午夜免费观看网址| 国产成人系列免费观看| 欧美+亚洲+日韩+国产| 久久香蕉精品热| 成人高潮视频无遮挡免费网站| 亚洲中文日韩欧美视频| 国产高清三级在线| 国产精品久久电影中文字幕| 亚洲欧美激情综合另类| 观看免费一级毛片| 国产精品一区二区免费欧美| 蜜桃久久精品国产亚洲av| 无人区码免费观看不卡| 午夜福利视频1000在线观看| 国产精品久久久久久人妻精品电影| 露出奶头的视频| 日韩欧美三级三区| 日本 av在线| 一级a爱片免费观看的视频| 村上凉子中文字幕在线| 亚洲人成网站在线播放欧美日韩| 亚洲欧美日韩东京热| 国产熟女xx| 免费一级毛片在线播放高清视频| svipshipincom国产片| 久久久国产成人精品二区| 99精品久久久久人妻精品| 久久亚洲真实| 精品免费久久久久久久清纯| or卡值多少钱| 变态另类丝袜制服| 亚洲熟妇中文字幕五十中出| 国产精品久久久久久精品电影| 在线视频色国产色| 嫩草影视91久久| 成年女人看的毛片在线观看| 国产精品1区2区在线观看.| 伊人久久精品亚洲午夜| 日韩欧美在线乱码| 成人欧美大片| 成人欧美大片| 窝窝影院91人妻| 90打野战视频偷拍视频| 国产高清视频在线观看网站| 免费av毛片视频| 国产精品香港三级国产av潘金莲| 人人妻人人澡欧美一区二区| 久久精品人妻少妇| 每晚都被弄得嗷嗷叫到高潮| 三级男女做爰猛烈吃奶摸视频| 婷婷亚洲欧美| 欧美日韩瑟瑟在线播放| 在线十欧美十亚洲十日本专区| 欧美日韩国产亚洲二区| 久久久久久久久中文| 波多野结衣巨乳人妻| 成人欧美大片| 成人精品一区二区免费| 一级毛片女人18水好多| 夜夜躁狠狠躁天天躁| xxx96com| 91九色精品人成在线观看| 在线国产一区二区在线| netflix在线观看网站| 亚洲中文日韩欧美视频| 免费看日本二区| 久久九九热精品免费| 国产蜜桃级精品一区二区三区| 精品国产三级普通话版| 久久亚洲精品不卡| 他把我摸到了高潮在线观看| www.色视频.com| 丰满人妻一区二区三区视频av | 国产精品精品国产色婷婷| 欧美一区二区亚洲| 国产黄色小视频在线观看| 亚洲欧美一区二区三区黑人| 岛国在线免费视频观看| 日日干狠狠操夜夜爽| 两个人视频免费观看高清| 免费搜索国产男女视频| 亚洲黑人精品在线| 在线观看午夜福利视频| 一边摸一边抽搐一进一小说| 午夜福利视频1000在线观看| 欧美色视频一区免费| 精品无人区乱码1区二区| xxxwww97欧美| 首页视频小说图片口味搜索| 国产v大片淫在线免费观看| 99热这里只有精品一区| av黄色大香蕉| 精品一区二区三区av网在线观看| 国产精品电影一区二区三区| 日日摸夜夜添夜夜添小说| 欧美高清成人免费视频www| 亚洲国产日韩欧美精品在线观看 | 少妇的逼好多水| 久久精品综合一区二区三区| 舔av片在线| 18禁在线播放成人免费| 久久6这里有精品| 亚洲精品久久国产高清桃花| 美女cb高潮喷水在线观看| 亚洲成av人片在线播放无| 国产午夜精品论理片| 伊人久久精品亚洲午夜| 少妇人妻一区二区三区视频| 久久香蕉国产精品| 国产黄色小视频在线观看| av女优亚洲男人天堂| 国产久久久一区二区三区| 国产免费男女视频| 美女免费视频网站| 欧美一区二区国产精品久久精品| 51午夜福利影视在线观看| 天天躁日日操中文字幕| 久久性视频一级片| 午夜激情欧美在线| 欧美另类亚洲清纯唯美| 亚洲无线观看免费| 999久久久精品免费观看国产| 夜夜看夜夜爽夜夜摸| 精品人妻偷拍中文字幕| 少妇人妻一区二区三区视频| 国产伦精品一区二区三区视频9 | 18禁在线播放成人免费| 老鸭窝网址在线观看| 日韩亚洲欧美综合| 美女高潮的动态| 精品人妻一区二区三区麻豆 | 久久久久九九精品影院| 成人av在线播放网站| 99久久精品国产亚洲精品| 欧美日韩亚洲国产一区二区在线观看| 九色国产91popny在线| 午夜日韩欧美国产| 国产乱人伦免费视频| 1000部很黄的大片| 国产真实伦视频高清在线观看 | 国产av不卡久久| 极品教师在线免费播放| 日韩av在线大香蕉| 午夜精品在线福利| 国产精品亚洲av一区麻豆| 国产在视频线在精品| 成人一区二区视频在线观看| 久久精品国产亚洲av涩爱 | 18禁裸乳无遮挡免费网站照片| 99在线视频只有这里精品首页| 女人被狂操c到高潮| 久久久久精品国产欧美久久久| 亚洲av中文字字幕乱码综合| 成人高潮视频无遮挡免费网站| 午夜福利在线观看免费完整高清在 | 午夜免费观看网址| 淫秽高清视频在线观看| 色精品久久人妻99蜜桃| 国产成人av教育| 亚洲不卡免费看| 天堂av国产一区二区熟女人妻| 免费无遮挡裸体视频| 午夜影院日韩av| 精品乱码久久久久久99久播| 偷拍熟女少妇极品色| a在线观看视频网站| 桃色一区二区三区在线观看| 午夜久久久久精精品| 国产午夜精品论理片| 精品国产超薄肉色丝袜足j| 日韩成人在线观看一区二区三区| 97超视频在线观看视频| 99国产综合亚洲精品| 欧美乱妇无乱码| 成熟少妇高潮喷水视频| 国产欧美日韩精品一区二区| 欧美激情久久久久久爽电影| 91久久精品电影网| 久久精品国产亚洲av香蕉五月| 久久久久久久午夜电影| 老鸭窝网址在线观看| 非洲黑人性xxxx精品又粗又长| 法律面前人人平等表现在哪些方面| 日日夜夜操网爽| 亚洲色图av天堂| 又粗又爽又猛毛片免费看| 色综合亚洲欧美另类图片| 极品教师在线免费播放| 欧美乱妇无乱码| 亚洲avbb在线观看| 9191精品国产免费久久| 亚洲精品在线美女| 国产精品香港三级国产av潘金莲| 丰满乱子伦码专区| 亚洲精品影视一区二区三区av| 欧美色欧美亚洲另类二区| 亚洲在线自拍视频| 久久久久国产精品人妻aⅴ院| ponron亚洲| 99久久精品热视频| 午夜亚洲福利在线播放| 麻豆国产av国片精品| 国产精品 欧美亚洲| 在线观看舔阴道视频| 三级国产精品欧美在线观看| 久久这里只有精品中国| 亚洲av电影在线进入| 亚洲欧美日韩高清专用| 亚洲不卡免费看| 熟妇人妻久久中文字幕3abv| 老汉色∧v一级毛片| www.色视频.com| 69人妻影院| 热99re8久久精品国产| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 国产成人影院久久av| 午夜福利成人在线免费观看| 天堂av国产一区二区熟女人妻| 精品不卡国产一区二区三区| 欧美在线黄色| 国产精品自产拍在线观看55亚洲| 精品欧美国产一区二区三| 国内久久婷婷六月综合欲色啪| 欧美+日韩+精品| 欧美一区二区精品小视频在线| 免费电影在线观看免费观看| 欧美成人性av电影在线观看| 在线看三级毛片| 真人做人爱边吃奶动态| 国产一区二区三区视频了| 他把我摸到了高潮在线观看| 日日干狠狠操夜夜爽| 桃色一区二区三区在线观看| a级一级毛片免费在线观看| 动漫黄色视频在线观看| 久久99热这里只有精品18| 在线观看美女被高潮喷水网站 | www.999成人在线观看| 99久久无色码亚洲精品果冻| 久久久久性生活片| 国产欧美日韩精品一区二区| 精品国产亚洲在线| 香蕉av资源在线| 亚洲精华国产精华精| 成人亚洲精品av一区二区| 特级一级黄色大片| 2021天堂中文幕一二区在线观| 18禁美女被吸乳视频| e午夜精品久久久久久久| 看黄色毛片网站| 日韩人妻高清精品专区| 亚洲熟妇熟女久久| 亚洲av免费高清在线观看| 国产 一区 欧美 日韩| 欧美日韩综合久久久久久 | 午夜亚洲福利在线播放| e午夜精品久久久久久久| 久久精品国产综合久久久| 国产三级中文精品| 哪里可以看免费的av片| 亚洲av免费高清在线观看| 黑人欧美特级aaaaaa片| 亚洲精品在线美女| 欧美黄色片欧美黄色片| 亚洲精品久久国产高清桃花| 高清在线国产一区| 国语自产精品视频在线第100页| 欧美bdsm另类| 欧美在线一区亚洲| 99热这里只有精品一区| 少妇高潮的动态图| av欧美777| 嫩草影院入口| 99久久精品一区二区三区| 国产精品综合久久久久久久免费| 精品乱码久久久久久99久播| 桃红色精品国产亚洲av| av黄色大香蕉| 亚洲中文字幕日韩| 99精品欧美一区二区三区四区| 久久九九热精品免费| 欧美一级a爱片免费观看看| 我的老师免费观看完整版| 国产av不卡久久| 国产成人aa在线观看| 国产亚洲精品综合一区在线观看| 国产av在哪里看| 欧美日韩乱码在线| 色av中文字幕| av福利片在线观看| 欧美午夜高清在线| 欧美成狂野欧美在线观看| 免费在线观看日本一区| 亚洲 国产 在线| av在线天堂中文字幕| 蜜桃亚洲精品一区二区三区| 成人国产综合亚洲| 男女做爰动态图高潮gif福利片| 在线观看午夜福利视频| 99热这里只有是精品50| 国产精品综合久久久久久久免费| 色综合站精品国产| 日本撒尿小便嘘嘘汇集6| 国产亚洲精品一区二区www| 18禁黄网站禁片免费观看直播| 日日干狠狠操夜夜爽| 国产av一区在线观看免费| 999久久久精品免费观看国产| 少妇的逼好多水| 老司机午夜十八禁免费视频| 欧美乱色亚洲激情| 国产精品av视频在线免费观看| 亚洲精品影视一区二区三区av| 国产综合懂色| 国内精品美女久久久久久| 90打野战视频偷拍视频| 久久久国产成人精品二区| 国内精品一区二区在线观看| 国产免费av片在线观看野外av| 成年人黄色毛片网站| 中文字幕av在线有码专区| 欧美av亚洲av综合av国产av| 亚洲黑人精品在线| 99热这里只有是精品50| 波多野结衣高清无吗| 91在线精品国自产拍蜜月 | 97人妻精品一区二区三区麻豆| 丝袜美腿在线中文| 99热这里只有是精品50| 波多野结衣高清无吗| 岛国在线免费视频观看| 欧美成人一区二区免费高清观看| 亚洲久久久久久中文字幕| 精品一区二区三区视频在线观看免费| 床上黄色一级片| 少妇的丰满在线观看| 变态另类成人亚洲欧美熟女| 国产国拍精品亚洲av在线观看 | 免费观看精品视频网站| 午夜免费成人在线视频| 亚洲欧美日韩卡通动漫| 日本一本二区三区精品| 99国产极品粉嫩在线观看| 亚洲aⅴ乱码一区二区在线播放| 最近最新中文字幕大全电影3| 成人国产一区最新在线观看| 日日干狠狠操夜夜爽| 99热这里只有是精品50| 99国产极品粉嫩在线观看| 夜夜爽天天搞| 可以在线观看毛片的网站| 男人舔奶头视频| 国产精品久久久久久人妻精品电影| 亚洲欧美日韩高清在线视频| 日本与韩国留学比较| 免费看十八禁软件| 少妇的逼水好多| 精品国产三级普通话版| av女优亚洲男人天堂| 亚洲va日本ⅴa欧美va伊人久久| 黄片小视频在线播放| 两人在一起打扑克的视频| 亚洲av二区三区四区| 国产一区二区激情短视频| 成人亚洲精品av一区二区| 国内精品美女久久久久久| 男女之事视频高清在线观看| 亚洲人成网站高清观看| 亚洲自拍偷在线| 国产亚洲精品久久久com| 在线播放无遮挡| 国产精品久久久久久精品电影| 美女 人体艺术 gogo| 白带黄色成豆腐渣| 色综合站精品国产| 精品日产1卡2卡| 免费在线观看影片大全网站| 男女之事视频高清在线观看| 国产成人a区在线观看| 91av网一区二区| av女优亚洲男人天堂| 在线天堂最新版资源| 国产精品女同一区二区软件 | 亚洲美女视频黄频| 女警被强在线播放| 免费av不卡在线播放| 久久久精品欧美日韩精品| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 成人18禁在线播放| 人人妻,人人澡人人爽秒播| 十八禁人妻一区二区| 欧美+日韩+精品| 国产av不卡久久| 日本在线视频免费播放| 午夜免费成人在线视频| 青草久久国产| 亚洲欧美日韩东京热| 亚洲国产欧美人成| 精品久久久久久久久久久久久| 最好的美女福利视频网| 老司机午夜十八禁免费视频| 不卡一级毛片| 欧美zozozo另类| 亚洲国产欧美网| 男女床上黄色一级片免费看| 亚洲av不卡在线观看| 亚洲 国产 在线| 男人的好看免费观看在线视频| 伊人久久精品亚洲午夜| 女生性感内裤真人,穿戴方法视频| 国产久久久一区二区三区| 日本五十路高清| 国产欧美日韩一区二区精品| 国产黄片美女视频| 香蕉丝袜av| 免费一级毛片在线播放高清视频| 国产精品久久久久久精品电影| 蜜桃亚洲精品一区二区三区| 国产免费一级a男人的天堂| 午夜福利在线观看吧| 欧美绝顶高潮抽搐喷水| 婷婷六月久久综合丁香| 男人的好看免费观看在线视频| 熟妇人妻久久中文字幕3abv| 成年女人永久免费观看视频| 日本一本二区三区精品| 久久这里只有精品中国| 首页视频小说图片口味搜索| 精品人妻偷拍中文字幕| 国产亚洲精品久久久com| 婷婷精品国产亚洲av| 免费看美女性在线毛片视频| a级一级毛片免费在线观看| www.色视频.com| 亚洲人成网站在线播| 亚洲国产精品999在线| x7x7x7水蜜桃| 日韩免费av在线播放| 欧美成狂野欧美在线观看| 在线观看66精品国产| 日日干狠狠操夜夜爽| 亚洲一区二区三区不卡视频| 丰满乱子伦码专区| 成人欧美大片| 国产国拍精品亚洲av在线观看 | 婷婷精品国产亚洲av| 最近最新中文字幕大全免费视频| 久久精品国产综合久久久| 国产一区二区亚洲精品在线观看| 天堂网av新在线| 在线免费观看的www视频| 尤物成人国产欧美一区二区三区| 在线观看一区二区三区| 亚洲第一欧美日韩一区二区三区| 亚洲精品一区av在线观看| 三级毛片av免费| 国内少妇人妻偷人精品xxx网站| 久久久色成人| 日本黄大片高清| 国产亚洲av嫩草精品影院| 中文亚洲av片在线观看爽| 一边摸一边抽搐一进一小说| 黄色日韩在线| 久久亚洲真实| 欧美日韩精品网址| 亚洲黑人精品在线| 99久久九九国产精品国产免费| 日日干狠狠操夜夜爽| 日本熟妇午夜| 两个人视频免费观看高清| 国产精品乱码一区二三区的特点| 少妇裸体淫交视频免费看高清| 国产熟女xx| 亚洲av一区综合| 久久精品国产亚洲av涩爱 | 香蕉久久夜色| 免费观看人在逋| 欧美绝顶高潮抽搐喷水| 一个人免费在线观看电影| 欧美+日韩+精品| 一夜夜www| 天天一区二区日本电影三级| 岛国在线免费视频观看| 999久久久精品免费观看国产| 亚洲成a人片在线一区二区| 黄色女人牲交| 真实男女啪啪啪动态图| 在线免费观看的www视频| 男女那种视频在线观看| www.999成人在线观看| 欧美最黄视频在线播放免费| 久久精品综合一区二区三区| 久久久久久久久中文| 亚洲性夜色夜夜综合|