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

    基于空間近似概率約束的混合密度網(wǎng)絡(luò)砂體厚度預(yù)測

    2020-07-25 02:55:48張運(yùn)龍
    石油物探 2020年4期
    關(guān)鍵詞:概率分布砂體儲(chǔ)層

    劉 可,尹 成,趙 虎,張運(yùn)龍

    (1.西南石油大學(xué)地球科學(xué)與技術(shù)學(xué)院,四川成都610500;2.中國石油新疆油田分公司勘探開發(fā)研究院,新疆克拉瑪依834000)

    神經(jīng)網(wǎng)絡(luò)作為一種由數(shù)據(jù)驅(qū)動(dòng)的非線性算法,因其“端對(duì)端”的特點(diǎn)以及強(qiáng)大的學(xué)習(xí)能力,廣泛應(yīng)用于圖像處理、語音識(shí)別、自然語言處理等領(lǐng)域[1]。根據(jù)神經(jīng)網(wǎng)絡(luò)所承擔(dān)任務(wù)的數(shù)學(xué)特點(diǎn)可分為:聚類任務(wù)、分類任務(wù)、回歸任務(wù)、重構(gòu)任務(wù)等。

    在地球物理領(lǐng)域,聚類任務(wù)主要包括地震相劃分[2],常使用自組織映射(SOM)神經(jīng)網(wǎng)絡(luò);分類任務(wù)包括儲(chǔ)層識(shí)別[3]、初至波識(shí)別[4]、斷層識(shí)別[5]等,常使用概率神經(jīng)網(wǎng)絡(luò)(PNN)和具有二維空間適應(yīng)能力的卷積神經(jīng)網(wǎng)絡(luò)(CNN)等;回歸任務(wù)包括巖石物理參數(shù)預(yù)測[6]、儲(chǔ)層參數(shù)預(yù)測[7]、測井曲線預(yù)測[8]等,往往使用前饋神經(jīng)網(wǎng)絡(luò)(BP)、徑向基神經(jīng)網(wǎng)絡(luò)(RBF)、循環(huán)神經(jīng)網(wǎng)絡(luò)(RNN)等;重構(gòu)任務(wù)包括疊前、疊后地震數(shù)據(jù)插值重構(gòu)[9]、地震數(shù)據(jù)去噪[10]等,主要使用卷積神經(jīng)網(wǎng)絡(luò)(CNN)、深度殘差網(wǎng)絡(luò)(RESnet)以及U型神經(jīng)網(wǎng)絡(luò)(U-net)等。

    在回歸任務(wù)中的儲(chǔ)層參數(shù)預(yù)測研究方面,除神經(jīng)網(wǎng)絡(luò)方法外,樹模型類方法[11]、支持向量機(jī)類方法[12]、多元回歸類方法[13]等方法也被廣泛應(yīng)用。在傳統(tǒng)方法基礎(chǔ)上,近年提出了以深度神經(jīng)決策森林(DNDF)[14]、自適應(yīng)神經(jīng)決策樹(ANT)[15]為代表的多方法融合深度學(xué)習(xí)模型。雖然這些模型多適用于分類任務(wù),但只須將模型進(jìn)行少許變形即可適用于回歸任務(wù),如由深度神經(jīng)決策森林(DNDF)發(fā)展而來的深度回歸森林(DRF)[16]。本質(zhì)而言,這些方法或模型均可視為一個(gè)復(fù)雜的、非線性的函數(shù),即每個(gè)輸入都有確定唯一的輸出與之對(duì)應(yīng),因此這些方法或模型僅適用于復(fù)雜唯一性問題。

    地震數(shù)據(jù)反演通常具有非唯一解,并且受噪聲等因素影響較大。由地震數(shù)據(jù)估計(jì)儲(chǔ)層參數(shù)具有明顯的不確定性,因此直接應(yīng)用上述方法或模型估計(jì)儲(chǔ)層參數(shù)會(huì)存在一定程度的誤差。利用概率分布解決不確定性問題效果較好。結(jié)合神經(jīng)網(wǎng)絡(luò)和概率分布兩者優(yōu)勢的混合密度網(wǎng)絡(luò)是目前解決不確定性問題的較好選擇。

    混合密度網(wǎng)絡(luò)(MDN)最早由BISHOP等[17]和CORNFORD等[18]提出,該網(wǎng)絡(luò)不僅擁有神經(jīng)網(wǎng)絡(luò)的學(xué)習(xí)能力,還擁有概率分布描述不確定性問題能力,已經(jīng)成功應(yīng)用于具有不確定性的海面風(fēng)場反演。陳坤堂等[19]和馮倩等[20]聯(lián)合BP神經(jīng)網(wǎng)絡(luò)和混合密度網(wǎng)絡(luò)成功完成海面風(fēng)場反演,證明了該網(wǎng)絡(luò)的有效性。SHAHRAEENI等[21]將該網(wǎng)絡(luò)引入巖石物理反演。殷習(xí)容等[22]通過該網(wǎng)絡(luò),利用橫縱波速度反演含水飽和度、泥質(zhì)含量和孔隙度等巖石物理參數(shù),同時(shí)測試該網(wǎng)絡(luò)抗噪聲能力,證明該網(wǎng)絡(luò)對(duì)巖石物理參數(shù)反演的有效性。但混合密度網(wǎng)絡(luò)在儲(chǔ)層參數(shù)預(yù)測方面的研究文獻(xiàn)較少。

    由于地震數(shù)據(jù)與儲(chǔ)層參數(shù)的不確定性關(guān)系,通過混合密度網(wǎng)絡(luò)得到的儲(chǔ)層參數(shù)“觀測”概率分布往往呈現(xiàn)“多峰”,為提升儲(chǔ)層參數(shù)預(yù)測的精確度,需要對(duì)該概率分布進(jìn)行約束。由于地下介質(zhì)具有連續(xù)分布的、空間橫向漸變性的特征,因而可以假設(shè)相鄰地震道儲(chǔ)層參數(shù)之間具有一定的相似性[23]。以該假設(shè)為基礎(chǔ),可以從已知井點(diǎn)的準(zhǔn)確值出發(fā),建立相鄰道的“估計(jì)”概率分布來約束“觀測”概率分布。相鄰道“估計(jì)”概率分布的約束,使得“觀測”概率分布的“多峰”性得到一定程度壓制。以往研究表明受“過擬合”等影響,單獨(dú)使用神經(jīng)網(wǎng)絡(luò)得到的結(jié)果往往分布雜亂,空間橫向連續(xù)性差[24]。因此反演中需考慮相鄰道儲(chǔ)層參數(shù),使得相鄰道之間儲(chǔ)層參數(shù)過渡更為平滑。

    本文在混合密度網(wǎng)絡(luò)的基礎(chǔ)上,以相鄰道橫向漸變?yōu)榧s束條件,提出了基于空間近似概率約束的混合密度網(wǎng)絡(luò)儲(chǔ)層參數(shù)預(yù)測方法。以砂體厚度預(yù)測為例,分別利用本文方法對(duì)理論模型和實(shí)際資料進(jìn)行了試驗(yàn),同時(shí)與BP神經(jīng)網(wǎng)絡(luò)和混合密度網(wǎng)絡(luò)的預(yù)測結(jié)果進(jìn)行了對(duì)比,并結(jié)合預(yù)測結(jié)果分析了3種方法在地震儲(chǔ)層參數(shù)預(yù)測方面的效果。

    1 算法原理

    1.1 混合密度網(wǎng)絡(luò)

    1.1.1 混合模型

    將多個(gè)不同頻率的正弦或余弦信號(hào)進(jìn)行組合可以得到一個(gè)復(fù)雜信號(hào)。同理,將多個(gè)簡單概率分布進(jìn)行組合也可以得到一個(gè)復(fù)雜概率分布,即混合概率分布模型。假設(shè)簡單概率分布為高斯概率分布(N),在地震數(shù)據(jù)反演中地震屬性(s)已知,則儲(chǔ)層參數(shù)(r)的概率分布可表述如下:

    (1)

    1.1.2 混合密度網(wǎng)絡(luò)結(jié)構(gòu)

    混合密度網(wǎng)絡(luò)數(shù)學(xué)模型來自于混合概率分布模型。與BP神經(jīng)網(wǎng)絡(luò)相似,混合密度網(wǎng)絡(luò)一般由4部分組成:輸入層,隱藏層,預(yù)處理層和輸出層(圖1),但其輸出層輸出的值并不是儲(chǔ)層參數(shù),而是一系列與混合概率分布模型相關(guān)的參數(shù)[25]。

    圖1 混合密度網(wǎng)絡(luò)的一般網(wǎng)絡(luò)結(jié)構(gòu)

    對(duì)于上述混合密度網(wǎng)絡(luò)結(jié)構(gòu),其前向過程可分為:輸入層-隱藏層、隱藏層-預(yù)處理層以及預(yù)處理層-輸出層3個(gè)部分。

    輸入層-隱藏層與BP神經(jīng)網(wǎng)絡(luò)一致,采取全連接方式計(jì)算:

    (2)

    隱藏層-預(yù)處理層也采取全連接方式計(jì)算,但不使用激活函數(shù):

    (3)

    根據(jù)高斯混合復(fù)雜概率分布的限制條件,對(duì)預(yù)處理層-輸出層不同類別節(jié)點(diǎn)采取不同處理方式:各高斯分布權(quán)值為正數(shù),權(quán)值之和為1以及方差為正數(shù)。則由混合密度網(wǎng)絡(luò)得到的“觀測”概率分布(Pobs)可表述如下:

    (4)

    根據(jù)最大似然估計(jì),混合密度網(wǎng)絡(luò)的誤差函數(shù)L一般為:

    (5)

    式中:Xd和Yd分別為訓(xùn)練集中第d個(gè)訓(xùn)練樣本的輸入和輸出;D為參與訓(xùn)練樣本數(shù)目。

    根據(jù)誤差函數(shù),使用BP算法進(jìn)行訓(xùn)練?;旌厦芏染W(wǎng)絡(luò)可視為BP神經(jīng)網(wǎng)絡(luò)的一種概率上的推廣[25],在單獨(dú)使用混合密度網(wǎng)絡(luò)時(shí),往往將預(yù)測區(qū)間內(nèi)概率最大處作為預(yù)測輸出。

    1.2 相鄰道估計(jì)模型

    以往神經(jīng)網(wǎng)絡(luò)僅僅單獨(dú)應(yīng)用于每一地震道,未考慮地層介質(zhì)在橫向上具有的漸變性,由此可假設(shè)相鄰地震道儲(chǔ)層參數(shù)之間具有一定相似性[23],在此基礎(chǔ)上建立相鄰道“估計(jì)”數(shù)學(xué)模型。對(duì)于二維地震數(shù)據(jù),每一地震道均具有唯一的一維坐標(biāo)n,其中n為該地震道的inline號(hào)。若第n道儲(chǔ)層參數(shù)rn未知,而第n-1道儲(chǔ)層參數(shù)rn-1已知,則rn的“估計(jì)”概率分布(Pest)可表述如下:

    (6)

    (7)

    式中:sn和sn-1為第n道和第n-1道的地震屬性;C為整個(gè)工區(qū)地震屬性組成矩陣的協(xié)方差矩陣;const為常數(shù),表示整個(gè)工區(qū)的基礎(chǔ)差異。

    當(dāng)?shù)趎-1道儲(chǔ)層參數(shù)rn-1不確定,僅知其概率分布為bel(rn-1)時(shí),則rn的“估計(jì)”概率分布(Pest)可表述如下:

    (8)

    實(shí)際地震數(shù)據(jù)多為三維數(shù)據(jù),對(duì)此每一地震道均具有唯一的二維坐標(biāo)(n,m),其中n和m分別代表其inline號(hào)和crossline號(hào)。仿照卡爾曼濾波相鄰估計(jì)方法,將其相鄰估計(jì)方法應(yīng)用于二維空間(圖2),其中藍(lán)色表示第一輪相鄰估計(jì),紅色表示第二輪相鄰估計(jì)??柭鼮V波相鄰估計(jì)多用于一維時(shí)間域,是單向、單點(diǎn)且不交叉的(圖2a),而二維空間域相鄰估計(jì)是多向、集體且存在交叉的(圖2b)。

    1.3 儲(chǔ)層參數(shù)預(yù)測

    對(duì)于儲(chǔ)層參數(shù)未知的(n,m)道,使用混合密度網(wǎng)絡(luò)可以得到該道儲(chǔ)層參數(shù)rn,m的“觀測”概率分布(Pobs),同時(shí)根據(jù)相鄰道“估計(jì)”關(guān)系((6)式和(8)式)得到(n,m)道儲(chǔ)層參數(shù)rn,m的“估計(jì)”概率分布(Pest)。由于二維空間域相鄰估計(jì)具有交叉的特點(diǎn)(圖2b),為簡便表達(dá),設(shè)(n,m)道的相鄰道隸屬于集合w,將兩種概率分布進(jìn)行融合,得:

    圖2 兩種相鄰估計(jì)a 卡爾曼濾波相鄰估計(jì); b 本文二維相鄰估計(jì)

    (9)

    式中;η為歸一化參數(shù);K′為融合后概率分布的高斯分布數(shù)目。

    顯然,Pmerge擁有的高斯分布數(shù)目遠(yuǎn)大于Pest和Pobs擁有的高斯分布數(shù)目。為運(yùn)算簡便,采取近似方法,即在Pmerge中選取具有最大的K″個(gè)權(quán)值的高斯分布,重新歸一化后組合,得到(n,m)道儲(chǔ)層參數(shù)rn,m的概率分布bel(rn,m):

    (10)

    以井位處作為起始位置,根據(jù)二維空間外推方案(圖2b),以及公式(6)、公式(8)、公式(9)及公式(10)求取相鄰道儲(chǔ)層參數(shù)概率分布,再以這些相鄰道為起始位置重復(fù)以上過程,直到完成全三維工區(qū)地震道的計(jì)算,即完成了全工區(qū)儲(chǔ)層參數(shù)概率預(yù)測。單井儲(chǔ)層參數(shù)概率預(yù)測流程如圖3所示。

    圖3 單井儲(chǔ)層參數(shù)概率預(yù)測流程

    實(shí)際工區(qū)有多口井,往往會(huì)產(chǎn)生多個(gè)不同初始點(diǎn)的全工區(qū)儲(chǔ)層參數(shù)概率預(yù)測結(jié)果。概率理論下的相鄰道估計(jì)過程中會(huì)不斷積累誤差,為了在減少誤差的同時(shí)融合不同初始點(diǎn)的反演結(jié)果,根據(jù)預(yù)測道與各初始位置的距離進(jìn)行加權(quán)平均,得到加權(quán)后的預(yù)測道儲(chǔ)層參數(shù)概率分布belavg(rn,m):

    (11)

    (12)

    2 理論模型試驗(yàn)

    設(shè)計(jì)的理論模型如圖4a所示,包括砂泥巖薄互層以及上、下兩套圍巖,其中目的層厚度50m。砂體分為類菱形、類平行四邊形以及凹凸形3類,地震波速度2350m/s,密度2.1g/cm3;泥巖地震波速度2420m/s,密度2.2g/cm3;圍巖地震波速度2500m/s,密度2.0g/cm3。圖4b為采用主頻為50Hz的零相位雷克子波進(jìn)行正演模擬的地震數(shù)據(jù)。

    對(duì)于該地震數(shù)據(jù)(圖4b),首先根據(jù)頂、底界面提取了振幅類、瞬時(shí)類、頻譜類等八大類60種地震屬性[26],同時(shí)采用粗糙集(RS)[27]、三參數(shù)敏感屬性分析(SDC)[28]等方法,聯(lián)合挑選出7種敏感地震屬性。然后隨機(jī)選取其中30%,20%,10%以及5%道數(shù)據(jù)作為訓(xùn)練數(shù)據(jù),最終得到如圖5至圖8所示的不同訓(xùn)練樣本數(shù)下砂體厚度(d)的預(yù)測結(jié)果及絕對(duì)誤差曲線。

    圖4 理論模型(a)及正演地震數(shù)據(jù)(b)

    圖5 3種方法的砂體厚度預(yù)測(a)及絕對(duì)誤差曲線(b)(30%訓(xùn)練數(shù)據(jù))

    圖6 3種方法的砂體厚度預(yù)測(a)及絕對(duì)誤差曲線(b)(20%訓(xùn)練數(shù)據(jù))

    對(duì)比本文方法(紅色曲線)、BP神經(jīng)網(wǎng)絡(luò)(藍(lán)色曲線)以及混合密度網(wǎng)絡(luò)(綠色曲線)這3種方法的預(yù)測結(jié)果可以看出,隨著訓(xùn)練樣本的減少,3種方法預(yù)測曲線越來越偏離實(shí)際厚度曲線。但是在相同條件下,本文方法與實(shí)際厚度曲線更為接近,橫向過渡更為連續(xù)。

    考慮到單次隨機(jī)選取訓(xùn)練數(shù)據(jù)不具有代表性,進(jìn)行50次訓(xùn)練數(shù)據(jù)的Monte Carlo隨機(jī)挑選,進(jìn)行3種方法對(duì)比預(yù)測試驗(yàn)。選擇預(yù)測砂體厚度曲線與實(shí)際砂體厚度曲線的平均誤差和相關(guān)系數(shù)作為預(yù)測效果的評(píng)價(jià)指標(biāo)。采用3種方法分別計(jì)算50次預(yù)測試驗(yàn)的平均誤差和平均相關(guān)系數(shù),結(jié)果如圖9所示。

    由圖9a可見,隨著訓(xùn)練樣本的減少,BP神經(jīng)網(wǎng)絡(luò)預(yù)測的平均誤差由0.6m上升到3.1m;混合密度網(wǎng)絡(luò)平均誤差由1.3m上升到3.9m;本文方法平均誤差由0.2m上升到2.1m。由圖9b可見,隨著訓(xùn)練樣本的減少,BP神經(jīng)網(wǎng)絡(luò)平均相關(guān)系數(shù)由0.98下降至0.83;混合密度網(wǎng)絡(luò)平均相關(guān)系數(shù)由0.89下降至0.67;本文方法平均相關(guān)系數(shù)由0.99下降至0.92。綜合圖5至圖9可見,在相同訓(xùn)練樣本下,單純混合密度網(wǎng)絡(luò)預(yù)測效果最差,BP神經(jīng)網(wǎng)絡(luò)較好,本文方法最好。

    圖7 3種方法的砂體厚度預(yù)測(a)及絕對(duì)誤差曲線(b)(15%訓(xùn)練數(shù)據(jù))

    圖8 3種方法的砂體厚度預(yù)測(a)及絕對(duì)誤差曲線(b)(5%訓(xùn)練數(shù)據(jù))

    圖9 50次Monte Carlo試驗(yàn)結(jié)果統(tǒng)計(jì)a 50次試驗(yàn)平均誤差; b 50次試驗(yàn)平均相關(guān)系數(shù)

    3 實(shí)際數(shù)據(jù)應(yīng)用

    基于上述實(shí)驗(yàn)結(jié)果,將本文方法應(yīng)用于X工區(qū)某實(shí)際地震數(shù)據(jù)。從工區(qū)提取的多種屬性中,同樣通過粗糙集(RS)、三參數(shù)敏感屬性分析(SDC)等方法聯(lián)合優(yōu)選出8種敏感屬性(圖10至圖13)。將本工區(qū)的56口井分為41口訓(xùn)練井(圖10至圖13中含Tr前綴井)和15口測試井(圖10至圖13中含Te前綴井)。

    圖10 基于敏感性分析得到的包絡(luò)差(a)和累計(jì)能量比(b)

    圖11 基于敏感性分析得到的工區(qū)能量半時(shí)間(a)和變異系數(shù)(b)

    圖12 基于敏感性分析得到的波形對(duì)稱度(a)和有效帶寬(b)

    圖13 基于敏感性分析得到的平均頻率(a)和平均振幅(b)

    在訓(xùn)練井位置處抽取地震道敏感屬性與訓(xùn)練井組合成訓(xùn)練集,訓(xùn)練BP神經(jīng)網(wǎng)絡(luò)和混合密度網(wǎng)絡(luò),將訓(xùn)練完成的網(wǎng)絡(luò)應(yīng)用于整個(gè)工區(qū),得到工區(qū)的砂體厚度分布預(yù)測結(jié)果(圖14a和圖14b)。再采用本文方法對(duì)整個(gè)工區(qū)進(jìn)行預(yù)測,得到工區(qū)的砂體厚度分布(圖14c)。從圖14可以看出,BP神經(jīng)網(wǎng)絡(luò)和混合密度網(wǎng)絡(luò)預(yù)測結(jié)果砂體厚度分布較為雜亂,橫向不連續(xù),而本文方法預(yù)測的砂體厚度橫向更為連續(xù),分布集中。

    圖14 3種方法預(yù)測的砂體厚度平面分布a BP神經(jīng)網(wǎng)絡(luò)預(yù)測結(jié)果; b 混合密度網(wǎng)絡(luò)預(yù)測結(jié)果; c 本文方法預(yù)測結(jié)果

    選取該工區(qū)測試井位處砂體預(yù)測厚度與已知的測試井位處的砂體厚度值進(jìn)行對(duì)比,其結(jié)果如圖15 所示。統(tǒng)計(jì)各個(gè)測試井?dāng)?shù)據(jù),BP神經(jīng)網(wǎng)絡(luò)預(yù)測的平均誤差為1.08m,誤差大于1m的井有7口,最大誤差為3.6m,位于井Te2處;混合密度網(wǎng)絡(luò)平均誤差為1.59m,誤差大于1m的井有7口,最大誤差為5.8m,位于井Te13處;本文方法平均誤差為0.7m,誤差大于1m的井有3口,最大誤差為2.5m,位于Te3處。由此可見本文方法預(yù)測的誤差最低。

    圖15 測試井點(diǎn)三種方法的預(yù)測結(jié)果a 3種方法預(yù)測的砂體厚度與井上實(shí)際厚度; b 3種方法的絕對(duì)誤差

    綜上所述,綜合判斷本文方法預(yù)測結(jié)果能夠反映工區(qū)砂體厚度的分布情況。

    4 結(jié)論

    傳統(tǒng)神經(jīng)網(wǎng)絡(luò)可視為一個(gè)復(fù)雜的、非線性的函數(shù),即每個(gè)輸入都有確定的輸出與之對(duì)應(yīng)。但是利用地震數(shù)據(jù)估計(jì)儲(chǔ)層參數(shù)具有明顯的不確定性,因而傳統(tǒng)神經(jīng)網(wǎng)絡(luò)不能很好適應(yīng)。對(duì)地震預(yù)測的不確定性問題,最好使用概率分布來間接解決。結(jié)合了神經(jīng)網(wǎng)絡(luò)和概率分布優(yōu)勢的混合密度網(wǎng)絡(luò),雖然從理論上能很好地解決多解性,但實(shí)際上的表現(xiàn)卻不如BP神經(jīng)網(wǎng)絡(luò),這是由于混合密度網(wǎng)絡(luò)需要訓(xùn)練的參數(shù)過多,以及概率分布本身所具有的不確定性造成的。

    本文提出的基于空間近似概率約束的混合密度網(wǎng)絡(luò)砂體厚度預(yù)測方法,結(jié)合了混合密度網(wǎng)絡(luò)和橫向漸變假設(shè),借鑒地質(zhì)統(tǒng)計(jì)學(xué)從井出發(fā)的觀點(diǎn),實(shí)現(xiàn)了概率分布下的儲(chǔ)層參數(shù)最優(yōu)估計(jì)。通過理論數(shù)據(jù)和實(shí)際資料的應(yīng)用測試,對(duì)比了BP神經(jīng)網(wǎng)絡(luò)、混合密度網(wǎng)絡(luò)以及本文方法的預(yù)測效果。結(jié)果表明,本文方法具有更高的精確性和較大的應(yīng)用潛力。

    本文方法中,地震道儲(chǔ)層參數(shù)概率分布需要由各已知井外推得到,而各已知井外推是相互獨(dú)立的運(yùn)算,所以需要較大次數(shù)的迭代運(yùn)算才能實(shí)現(xiàn)所有已知井的外推,如何實(shí)現(xiàn)多井同步外推是下一步研究的重點(diǎn)。

    猜你喜歡
    概率分布砂體儲(chǔ)層
    輸導(dǎo)層
    ——北美又一種非常規(guī)儲(chǔ)層類型
    河流相復(fù)合砂體不連續(xù)界限等效表征方法探討
    離散型概率分布的ORB圖像特征點(diǎn)誤匹配剔除算法
    CSAMT法在柴北緣砂巖型鈾礦勘查砂體探測中的應(yīng)用
    基于儲(chǔ)層構(gòu)型研究的儲(chǔ)層平面非均質(zhì)性表征
    關(guān)于概率分布函數(shù)定義的辨析
    科技視界(2016年19期)2017-05-18 10:18:46
    基于概率分布的PPP項(xiàng)目風(fēng)險(xiǎn)承擔(dān)支出測算
    低滲透儲(chǔ)層核磁共振可動(dòng)流體研究——以姬塬地區(qū)長6儲(chǔ)層為例
    砂體構(gòu)型對(duì)剩余油分布控制研究—以文中油田文25東油藏為例
    河南科技(2014年4期)2014-02-27 14:07:25
    一種相依極小P值統(tǒng)計(jì)量概率分布的近似計(jì)算方法
    国产精品一区二区精品视频观看| av有码第一页| 这个男人来自地球电影免费观看| www.精华液| 黑人操中国人逼视频| 飞空精品影院首页| 亚洲熟妇熟女久久| 两性夫妻黄色片| 国内久久婷婷六月综合欲色啪| 在线观看免费高清a一片| 麻豆成人av在线观看| av在线播放免费不卡| 亚洲九九香蕉| 中文字幕色久视频| 欧美性长视频在线观看| 国产xxxxx性猛交| 国内毛片毛片毛片毛片毛片| 欧美黄色淫秽网站| 丝袜美足系列| 纯流量卡能插随身wifi吗| 在线免费观看的www视频| 亚洲国产欧美一区二区综合| 免费av中文字幕在线| 色老头精品视频在线观看| 国产野战对白在线观看| 精品国产一区二区三区久久久樱花| 国产欧美日韩综合在线一区二区| 精品少妇久久久久久888优播| 国产亚洲av高清不卡| 久久久精品免费免费高清| 国产成人精品无人区| 精品久久久久久电影网| av电影中文网址| 国产在线一区二区三区精| 黄色视频不卡| bbb黄色大片| 国产99久久九九免费精品| 侵犯人妻中文字幕一二三四区| 国产真人三级小视频在线观看| 久久国产精品大桥未久av| 国产欧美日韩综合在线一区二区| 国产成人精品久久二区二区免费| 国产亚洲精品第一综合不卡| а√天堂www在线а√下载 | 亚洲中文字幕日韩| 美女高潮喷水抽搐中文字幕| 成年人黄色毛片网站| 国产在线精品亚洲第一网站| 在线观看66精品国产| 美女视频免费永久观看网站| 免费在线观看完整版高清| 桃红色精品国产亚洲av| 黑人欧美特级aaaaaa片| 日韩免费av在线播放| 亚洲精品美女久久av网站| 两性夫妻黄色片| 黄色a级毛片大全视频| 亚洲精品一卡2卡三卡4卡5卡| 亚洲第一av免费看| 亚洲av欧美aⅴ国产| 午夜成年电影在线免费观看| 一级a爱视频在线免费观看| 亚洲国产看品久久| 在线观看午夜福利视频| 精品久久久久久久久久免费视频 | 精品少妇一区二区三区视频日本电影| 他把我摸到了高潮在线观看| 后天国语完整版免费观看| 免费观看精品视频网站| av线在线观看网站| 亚洲欧美一区二区三区黑人| 国产精品 欧美亚洲| 午夜激情av网站| 国产主播在线观看一区二区| 两性夫妻黄色片| 免费观看精品视频网站| 国产亚洲精品第一综合不卡| 亚洲黑人精品在线| xxx96com| 亚洲av成人不卡在线观看播放网| 99精国产麻豆久久婷婷| 国产激情久久老熟女| 国产精品98久久久久久宅男小说| 国产精品二区激情视频| 欧美日韩一级在线毛片| 国精品久久久久久国模美| 久久久久久久国产电影| 中文字幕人妻丝袜制服| 欧美国产精品一级二级三级| 999精品在线视频| 国产一区二区激情短视频| 在线观看日韩欧美| 亚洲专区字幕在线| 少妇被粗大的猛进出69影院| 亚洲人成电影免费在线| 精品卡一卡二卡四卡免费| 日本五十路高清| 亚洲av第一区精品v没综合| 无人区码免费观看不卡| 91麻豆av在线| 久久国产精品人妻蜜桃| 王馨瑶露胸无遮挡在线观看| 在线观看舔阴道视频| 亚洲精品成人av观看孕妇| 乱人伦中国视频| 国产精品国产av在线观看| 亚洲av电影在线进入| 国产亚洲精品久久久久5区| 18禁美女被吸乳视频| 久久午夜亚洲精品久久| 久久精品国产99精品国产亚洲性色 | 国产成人啪精品午夜网站| 精品福利观看| 天堂动漫精品| 成人国语在线视频| 老司机在亚洲福利影院| 啦啦啦 在线观看视频| x7x7x7水蜜桃| 成年人午夜在线观看视频| 精品人妻在线不人妻| 久久精品亚洲av国产电影网| 美国免费a级毛片| 成年动漫av网址| 99riav亚洲国产免费| 两个人看的免费小视频| 久久精品国产a三级三级三级| 免费女性裸体啪啪无遮挡网站| 人妻久久中文字幕网| 成年版毛片免费区| 亚洲熟妇熟女久久| 国产亚洲欧美精品永久| 激情视频va一区二区三区| 国产欧美亚洲国产| 欧美成人免费av一区二区三区 | 岛国毛片在线播放| 亚洲精品国产一区二区精华液| 变态另类成人亚洲欧美熟女 | 精品一品国产午夜福利视频| 青草久久国产| 国产免费男女视频| 一本一本久久a久久精品综合妖精| 18禁黄网站禁片午夜丰满| 无限看片的www在线观看| 亚洲片人在线观看| 波多野结衣一区麻豆| 黑人操中国人逼视频| 国产精品香港三级国产av潘金莲| 在线观看日韩欧美| 国产乱人伦免费视频| 亚洲第一青青草原| 性少妇av在线| 亚洲三区欧美一区| 欧美激情高清一区二区三区| 91老司机精品| 国产一区在线观看成人免费| 黄片大片在线免费观看| 国产在视频线精品| 男女床上黄色一级片免费看| 女人被狂操c到高潮| 国产精品乱码一区二三区的特点 | 纯流量卡能插随身wifi吗| 欧美激情高清一区二区三区| 精品久久久久久久久久免费视频 | av欧美777| 18在线观看网站| 美女福利国产在线| 一进一出好大好爽视频| 色婷婷av一区二区三区视频| 黄网站色视频无遮挡免费观看| videosex国产| 久久热在线av| 国产野战对白在线观看| 欧美色视频一区免费| 99久久99久久久精品蜜桃| 精品一区二区三卡| 精品国产一区二区久久| 变态另类成人亚洲欧美熟女 | 美女高潮喷水抽搐中文字幕| 久久久久久亚洲精品国产蜜桃av| 精品一区二区三区视频在线观看免费 | avwww免费| 久久九九热精品免费| 欧美日韩国产mv在线观看视频| 国产高清激情床上av| 免费观看a级毛片全部| 超碰成人久久| 美女高潮到喷水免费观看| 成年女人毛片免费观看观看9 | 亚洲黑人精品在线| 91在线观看av| 人人妻人人爽人人添夜夜欢视频| 国产亚洲欧美精品永久| 人人妻,人人澡人人爽秒播| 日韩人妻精品一区2区三区| 国精品久久久久久国模美| 午夜福利欧美成人| 亚洲第一青青草原| 两人在一起打扑克的视频| 97人妻天天添夜夜摸| 国产一区二区三区综合在线观看| 国产精品永久免费网站| 91字幕亚洲| 一级片免费观看大全| 国精品久久久久久国模美| 午夜福利欧美成人| av片东京热男人的天堂| 又紧又爽又黄一区二区| 老汉色av国产亚洲站长工具| 国产精品久久电影中文字幕 | 十八禁人妻一区二区| 王馨瑶露胸无遮挡在线观看| 老汉色av国产亚洲站长工具| 777米奇影视久久| 婷婷成人精品国产| 精品免费久久久久久久清纯 | 欧美日韩乱码在线| 中文字幕av电影在线播放| 国产又色又爽无遮挡免费看| 午夜老司机福利片| 麻豆乱淫一区二区| 啪啪无遮挡十八禁网站| 香蕉久久夜色| 亚洲精品国产精品久久久不卡| 熟女少妇亚洲综合色aaa.| 老熟妇仑乱视频hdxx| 大型av网站在线播放| 黄色毛片三级朝国网站| 人人妻人人爽人人添夜夜欢视频| 多毛熟女@视频| 好男人电影高清在线观看| 电影成人av| 在线国产一区二区在线| 十八禁网站免费在线| 校园春色视频在线观看| 亚洲 国产 在线| 久久中文字幕人妻熟女| 欧美丝袜亚洲另类 | 人人澡人人妻人| 欧美成人午夜精品| 久久天堂一区二区三区四区| 最近最新免费中文字幕在线| 人妻久久中文字幕网| 欧美日韩瑟瑟在线播放| 91成人精品电影| 91av网站免费观看| 操出白浆在线播放| 看免费av毛片| 飞空精品影院首页| 久久精品aⅴ一区二区三区四区| 乱人伦中国视频| 久久中文字幕一级| 精品久久蜜臀av无| 高清欧美精品videossex| 天天影视国产精品| 天天操日日干夜夜撸| 久久草成人影院| 精品少妇一区二区三区视频日本电影| 亚洲aⅴ乱码一区二区在线播放 | 亚洲三区欧美一区| 1024视频免费在线观看| 亚洲专区字幕在线| 天天添夜夜摸| 热re99久久国产66热| 宅男免费午夜| 国产成人啪精品午夜网站| 极品少妇高潮喷水抽搐| 嫁个100分男人电影在线观看| 最近最新中文字幕大全电影3 | 露出奶头的视频| 人人澡人人妻人| 午夜福利乱码中文字幕| 91麻豆av在线| 极品少妇高潮喷水抽搐| 欧美日韩福利视频一区二区| 一级a爱视频在线免费观看| 91av网站免费观看| 精品久久久精品久久久| videosex国产| 美女福利国产在线| 亚洲成av片中文字幕在线观看| 亚洲精品国产一区二区精华液| 十八禁高潮呻吟视频| 欧美日韩视频精品一区| 久久久精品区二区三区| 久久久久国内视频| 老司机亚洲免费影院| 午夜福利在线免费观看网站| 91精品国产国语对白视频| 精品乱码久久久久久99久播| 亚洲精品国产一区二区精华液| 国产日韩欧美亚洲二区| av片东京热男人的天堂| 精品人妻1区二区| 欧美色视频一区免费| 精品少妇一区二区三区视频日本电影| 色精品久久人妻99蜜桃| 成年版毛片免费区| av网站免费在线观看视频| 午夜91福利影院| 99re在线观看精品视频| 午夜精品在线福利| 人成视频在线观看免费观看| 亚洲精品国产一区二区精华液| 麻豆乱淫一区二区| 天天添夜夜摸| 两个人看的免费小视频| 亚洲欧美色中文字幕在线| 国产av又大| 亚洲国产精品合色在线| 十八禁网站免费在线| 在线av久久热| 久久久精品区二区三区| 中文字幕另类日韩欧美亚洲嫩草| 精品一区二区三区av网在线观看| 欧美在线黄色| 18禁黄网站禁片午夜丰满| 一级毛片高清免费大全| 精品一区二区三卡| 国产欧美日韩综合在线一区二区| 久久精品亚洲熟妇少妇任你| 美女高潮喷水抽搐中文字幕| 十分钟在线观看高清视频www| 下体分泌物呈黄色| 宅男免费午夜| 国产成人av教育| 亚洲va日本ⅴa欧美va伊人久久| 天天躁夜夜躁狠狠躁躁| 国产有黄有色有爽视频| 午夜成年电影在线免费观看| 日本黄色视频三级网站网址 | 精品欧美一区二区三区在线| 下体分泌物呈黄色| 日韩欧美一区视频在线观看| 两性夫妻黄色片| 在线观看舔阴道视频| 国产主播在线观看一区二区| 欧美日韩视频精品一区| 好看av亚洲va欧美ⅴa在| av欧美777| e午夜精品久久久久久久| aaaaa片日本免费| 一进一出抽搐gif免费好疼 | 窝窝影院91人妻| 80岁老熟妇乱子伦牲交| 男女高潮啪啪啪动态图| 亚洲精品一二三| 亚洲欧洲精品一区二区精品久久久| 欧美黑人精品巨大| 亚洲欧洲精品一区二区精品久久久| 久久久精品免费免费高清| 免费久久久久久久精品成人欧美视频| 久热这里只有精品99| 国产激情久久老熟女| 午夜福利欧美成人| 欧美日韩亚洲国产一区二区在线观看 | 久久精品人人爽人人爽视色| 一区二区三区国产精品乱码| 大陆偷拍与自拍| 在线国产一区二区在线| 久久中文字幕人妻熟女| 国产精品乱码一区二三区的特点 | 亚洲av电影在线进入| 精品高清国产在线一区| 成在线人永久免费视频| 无限看片的www在线观看| 欧美激情高清一区二区三区| 免费日韩欧美在线观看| 亚洲黑人精品在线| 亚洲精品成人av观看孕妇| 久久精品成人免费网站| 亚洲中文字幕日韩| 欧美色视频一区免费| 亚洲精品国产区一区二| 波多野结衣一区麻豆| 欧美亚洲日本最大视频资源| 久99久视频精品免费| 好男人电影高清在线观看| 99久久国产精品久久久| 村上凉子中文字幕在线| 婷婷丁香在线五月| 国产高清激情床上av| 国产精品永久免费网站| 亚洲男人天堂网一区| 人人妻人人添人人爽欧美一区卜| 精品久久久久久,| 别揉我奶头~嗯~啊~动态视频| 欧美另类亚洲清纯唯美| 91精品国产国语对白视频| 99re在线观看精品视频| 久久久国产成人免费| 亚洲专区国产一区二区| 成人手机av| 999精品在线视频| 欧美日韩一级在线毛片| 十八禁高潮呻吟视频| 久久久久久亚洲精品国产蜜桃av| 亚洲国产精品sss在线观看 | 亚洲avbb在线观看| 亚洲视频免费观看视频| 91老司机精品| 久久久国产精品麻豆| 91在线观看av| 久久久久久久精品吃奶| 午夜两性在线视频| 在线观看免费午夜福利视频| 涩涩av久久男人的天堂| 婷婷丁香在线五月| 国产不卡av网站在线观看| 国产精品一区二区免费欧美| 成年人免费黄色播放视频| 看黄色毛片网站| 欧美日韩亚洲综合一区二区三区_| 欧美 亚洲 国产 日韩一| 久久99一区二区三区| 国产精品免费视频内射| 亚洲熟女精品中文字幕| 日韩欧美国产一区二区入口| 精品欧美一区二区三区在线| 国产欧美日韩一区二区三| 黄片播放在线免费| 欧美日韩中文字幕国产精品一区二区三区 | 午夜亚洲福利在线播放| 欧美精品亚洲一区二区| 久久九九热精品免费| 中文字幕制服av| 亚洲欧美激情综合另类| 亚洲国产看品久久| av天堂久久9| 中文字幕人妻丝袜一区二区| 王馨瑶露胸无遮挡在线观看| 满18在线观看网站| 建设人人有责人人尽责人人享有的| 欧美人与性动交α欧美精品济南到| 纯流量卡能插随身wifi吗| 操美女的视频在线观看| 国产主播在线观看一区二区| 人人妻人人添人人爽欧美一区卜| 老司机影院毛片| 老司机亚洲免费影院| 亚洲一区二区三区欧美精品| 欧美日韩亚洲高清精品| 国产单亲对白刺激| 另类亚洲欧美激情| 国产激情欧美一区二区| 久久人妻熟女aⅴ| 亚洲精品成人av观看孕妇| av网站免费在线观看视频| 亚洲成人免费电影在线观看| 一级作爱视频免费观看| 国产亚洲精品第一综合不卡| 国产激情久久老熟女| 亚洲第一青青草原| 在线国产一区二区在线| 免费人成视频x8x8入口观看| 黑人巨大精品欧美一区二区mp4| 欧美成狂野欧美在线观看| 嫩草影视91久久| 女人被狂操c到高潮| av欧美777| 人人澡人人妻人| 露出奶头的视频| 色婷婷久久久亚洲欧美| 久久性视频一级片| 亚洲精品在线观看二区| av网站在线播放免费| 免费在线观看完整版高清| 一级毛片高清免费大全| 99国产精品一区二区蜜桃av | 极品人妻少妇av视频| 美女午夜性视频免费| 最近最新中文字幕大全电影3 | 免费观看精品视频网站| 国产有黄有色有爽视频| 国产精品亚洲一级av第二区| 热99re8久久精品国产| 久久国产精品男人的天堂亚洲| 丰满人妻熟妇乱又伦精品不卡| 老汉色av国产亚洲站长工具| 精品久久久精品久久久| 在线观看午夜福利视频| 成人国产一区最新在线观看| 美女扒开内裤让男人捅视频| 久久久久精品人妻al黑| 欧美 亚洲 国产 日韩一| 满18在线观看网站| 999久久久精品免费观看国产| 亚洲人成伊人成综合网2020| 国产又色又爽无遮挡免费看| 一级毛片女人18水好多| 日韩欧美免费精品| 日韩中文字幕欧美一区二区| 操美女的视频在线观看| 一区二区日韩欧美中文字幕| 欧美成人免费av一区二区三区 | 亚洲精品国产区一区二| 大陆偷拍与自拍| 精品午夜福利视频在线观看一区| videos熟女内射| 两人在一起打扑克的视频| 精品电影一区二区在线| 精品熟女少妇八av免费久了| 久久人人97超碰香蕉20202| 亚洲欧美激情综合另类| 久久国产乱子伦精品免费另类| 成年女人毛片免费观看观看9 | 老司机影院毛片| 欧美 亚洲 国产 日韩一| 免费高清在线观看日韩| 一级毛片精品| 国产一区二区激情短视频| 窝窝影院91人妻| 在线天堂中文资源库| 欧美乱色亚洲激情| 一区二区三区精品91| 人人妻,人人澡人人爽秒播| 色尼玛亚洲综合影院| 国产亚洲欧美在线一区二区| 美女国产高潮福利片在线看| 亚洲人成电影免费在线| av欧美777| 丁香六月欧美| 一级片免费观看大全| 亚洲全国av大片| 亚洲精品国产一区二区精华液| 成年人免费黄色播放视频| 国精品久久久久久国模美| 亚洲av日韩在线播放| 国产无遮挡羞羞视频在线观看| 国产精品亚洲av一区麻豆| 日韩欧美三级三区| 亚洲自偷自拍图片 自拍| 久久午夜亚洲精品久久| 巨乳人妻的诱惑在线观看| 精品一区二区三区av网在线观看| 男女之事视频高清在线观看| 超碰97精品在线观看| 精品欧美一区二区三区在线| 亚洲精品美女久久av网站| 欧美激情久久久久久爽电影 | 国产精品 国内视频| 丰满饥渴人妻一区二区三| 国产成人免费观看mmmm| 美女高潮喷水抽搐中文字幕| 超碰成人久久| 中文字幕人妻熟女乱码| 日韩欧美一区视频在线观看| 色播在线永久视频| 国产成人影院久久av| 精品人妻1区二区| 最近最新中文字幕大全电影3 | 精品国产国语对白av| 这个男人来自地球电影免费观看| 精品亚洲成a人片在线观看| 久久影院123| 欧美激情高清一区二区三区| 身体一侧抽搐| 人妻久久中文字幕网| 美女高潮到喷水免费观看| 欧美日韩av久久| 精品无人区乱码1区二区| 亚洲成人国产一区在线观看| 高清欧美精品videossex| 久9热在线精品视频| av欧美777| 亚洲av片天天在线观看| 男人的好看免费观看在线视频 | 久久亚洲真实| 久久精品国产a三级三级三级| 久久影院123| 51午夜福利影视在线观看| 中文字幕制服av| av一本久久久久| 精品国产一区二区三区久久久樱花| 精品国产超薄肉色丝袜足j| e午夜精品久久久久久久| 99re6热这里在线精品视频| 激情在线观看视频在线高清 | 久久精品aⅴ一区二区三区四区| 亚洲av片天天在线观看| 免费看十八禁软件| 老司机靠b影院| 热99re8久久精品国产| 国产亚洲av高清不卡| 久热这里只有精品99| www.熟女人妻精品国产| 国产精品二区激情视频| 亚洲午夜理论影院| 日韩欧美免费精品| 欧美精品亚洲一区二区| 男人的好看免费观看在线视频 | 国产高清videossex| 亚洲精品美女久久av网站| 韩国av一区二区三区四区| 69av精品久久久久久| 亚洲欧美激情在线| 18禁黄网站禁片午夜丰满| 搡老岳熟女国产| 最新美女视频免费是黄的| 久久香蕉国产精品| 不卡av一区二区三区| 久久久久国内视频| 午夜精品国产一区二区电影| 狂野欧美激情性xxxx| 午夜免费观看网址| 中文字幕高清在线视频| 人成视频在线观看免费观看| 妹子高潮喷水视频| 日韩一卡2卡3卡4卡2021年| 久久精品熟女亚洲av麻豆精品| 久久国产精品人妻蜜桃| 一个人免费在线观看的高清视频| 又黄又爽又免费观看的视频| 又大又爽又粗| 精品一区二区三卡| 咕卡用的链子|