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

    海量地震采集資料現(xiàn)場質(zhì)量評價方法探討

    2021-12-04 02:15:32徐雷良徐維秀
    石油地球物理勘探 2021年6期
    關鍵詞:單炮工區(qū)分類

    徐雷良 徐維秀

    (中石化石油工程地球物理有限公司勝利分公司,山東東營 257000)

    0 引言

    地震勘探技術(shù)的發(fā)展對地震資料質(zhì)量監(jiān)控技術(shù)的要求越來越高。早期,僅憑肉眼和經(jīng)驗查看監(jiān)視記錄每天即可完成二維地震數(shù)十道、上百道的單炮質(zhì)量監(jiān)控。地震勘探技術(shù)與地震采集設備以及計算機技術(shù)的發(fā)展催生了現(xiàn)場地震資料處理技術(shù)[1],并以此作為現(xiàn)場延時質(zhì)控的重要手段。隨著高精度、高密度三維地震勘探技術(shù)的出現(xiàn)以及地震采集設備的不斷升級,陸續(xù)出現(xiàn)了地震資料現(xiàn)場實時質(zhì)控技術(shù)[2-3],并研發(fā)了大量軟件成果[4-7],完全改變了依靠監(jiān)視記錄和部分資料抽檢進行人工定性監(jiān)控模式,現(xiàn)場實時質(zhì)控技術(shù)具有比較全面、定量化及相對科學的特性[7]。以小面元、大道數(shù)為典型特征的“兩寬一高”與高效地震采集技術(shù)的發(fā)展應用在持續(xù)推動技術(shù)進步的同時,也引起了數(shù)據(jù)量的急劇攀升,亦即俗稱的“海量地震數(shù)據(jù)”。近年來,三維區(qū)塊部署的單炮數(shù)據(jù)量動輒數(shù)百兆字節(jié),而相鄰兩炮激發(fā)的時差已在半分鐘內(nèi)(井炮),甚至僅有數(shù)秒(可控震源),因此傳統(tǒng)定性質(zhì)控模式很難有效監(jiān)控海量地震資料,現(xiàn)有的定量化監(jiān)控技術(shù)與評價軟件面臨巨大挑戰(zhàn)[7-8]。地震資料評價技術(shù)是質(zhì)量監(jiān)控技術(shù)的重要組成部分,二者同步發(fā)展。從最初人工定性分析發(fā)展到基于標準進行多因素評價[4-5,7],再發(fā)展為基于單個屬性的資料評價[2-3],直到如今多元地震屬性綜合統(tǒng)計分析評價[6-7],逐步形成了由激發(fā)、接收、環(huán)境噪聲及地震屬性等全方位要素參與的評價體系,特別是基于地震屬性的資料評價模式已由單炮記錄面貌轉(zhuǎn)向內(nèi)部特征,為目前油氣勘探開發(fā)所急需的高精度、高分辨率地震資料提供了更為可靠的質(zhì)量保障。地震記錄一般是在連續(xù)地表與地質(zhì)條件及相同激發(fā)與接收環(huán)境下獲得的,因此各炮之間存在天然的、隱性的聯(lián)系,需要一種自學習算法尋找其潛在的關系,以快速完成單炮記錄品質(zhì)評價。目前,人工智能[9-10]已在許多領域取得成功,顯著地改變了人們的生產(chǎn)、生活方式。它在地震資料評價方面也有一定應用成果,但其著眼點主要是針對地震成果數(shù)據(jù)[11-13],在現(xiàn)場資料評價方面的應用才剛剛起步[14]。本文首先討論了目前常用的單炮記錄評價理論方法,接著探討了基于隨機森林(Random Forest,RF)[15-16]的海量地震資料智能評價流程及其關鍵技術(shù)。

    1 多因素確定性評價模型

    國家標準[17]和一些行業(yè)標準、企業(yè)標準中對地震資料分級制定了明確的技術(shù)規(guī)范,這些規(guī)范是對地震儀器、設備、激發(fā)與接收等各種因素性能指標的限定,是多年地震勘探實踐所形成的確定性指標。例如,TB時差、不正常道數(shù)量或占比、斷排列數(shù)、震源畸變超限等。根據(jù)標準,把這些參數(shù)硬性指標的考核稱為確定性評價模型。目前,通用或商用監(jiān)控軟件[7,18-20]全部采用了該模型。這類確定性評價模型專注于施工因素監(jiān)控,但并不能充分反映地表與地質(zhì)、環(huán)境以及偶然因素對地震資料的影響,該評價模式并不對單炮記錄品質(zhì)進行評判。

    2 數(shù)理統(tǒng)計評價模型

    2.1 單屬性評價模型

    這種模型是通過新老資料對比確定地震屬性閾值,利用單一屬性的閾值是否超限進行資料分級[3],該模型的優(yōu)勢在于針對性強,對于特定油氣勘探目標資料品質(zhì)分析具有指導意義。但不同地震屬性體現(xiàn)資料的不同特征,且每種地震屬性對資料品質(zhì)的反映具有片面性和模糊性,依賴單一地震屬性判定原始單炮記錄合格與否顯然是不科學的。

    2.2 基于多元屬性判別分析的評價模型

    多元屬性單炮記錄評價的實質(zhì)是多元判別分析問題,以下先討論兩級分類。

    (1)

    (2)

    式中Σ1、Σ2分別為G1和G2屬性的協(xié)方差矩陣。

    生產(chǎn)炮合格判定的準則可描述為

    (3)

    (4)

    式(4)為單炮記錄兩級分類判別準則。如果做三級分類,可增加一個閾值。

    實際生產(chǎn)中,也可采用品質(zhì)好的生產(chǎn)炮替代試驗炮。為了確保協(xié)方差矩陣的秩存在,炮數(shù)應大于地震屬性個數(shù),一般選擇最近采集的20~30炮為宜??紤]到多線束施工模式,這些炮應兼顧到每束線、每個排列,空間分布應相對均勻,從而使其更具代表性。如果選取一條優(yōu)質(zhì)生產(chǎn)炮記錄S(稱為標準記錄)作為G1的重心,判別準則變?yōu)?/p>

    (5)

    式中:μi、si分別表示生產(chǎn)炮X和S炮的第i個屬性;θi為設定的第i個屬性的閾值,且有θi≥0。

    所有與S炮具有相似地表和地下地質(zhì)條件、相同激發(fā)與接收及環(huán)境因素的生產(chǎn)炮以式(5)為判別準則。但是,一旦與以上條件不符,如不同激發(fā)藥量,就需要建立新的炮集重心。因此,同一個工區(qū)可建立多個區(qū)域,每個區(qū)域形成一個炮集G1,每個炮集只有一條標準記錄(重心)。圖1為據(jù)此開發(fā)的監(jiān)控軟件M[8]所設計的多區(qū)域海量地震資料評價流程。

    多元屬性評價模型理念易于接受、方便計算機實現(xiàn),適于現(xiàn)場實時監(jiān)控與評價。工程技術(shù)人員利用經(jīng)驗并結(jié)合老資料可檢驗分類效果,通過標準差、方差對監(jiān)控結(jié)果進行統(tǒng)計分析,并將結(jié)果疊合在遙感影像上,便于全區(qū)資料分析,指導地震采集施工。圖2為M軟件在某工區(qū)的合格炮能量分布圖(背景為衛(wèi)星照片),可見單炮能量與沙丘高程呈正相關。

    圖2 某工區(qū)合格單炮能量分布(背景為衛(wèi)星照片)

    多元屬性評價單炮的分類模式隱含著“異類即錯”的思想,目的在于嚴格實時監(jiān)控廢炮,以便及時補炮。該模型的主要局限在于閾值θi的定義,主觀性太強,難以適應多變的地震地質(zhì)條件,生產(chǎn)中不易把握且極易引起爭議。

    3 單炮記錄智能評價模型

    人工智能是利用數(shù)字計算機或者數(shù)字計算機控制的機器開發(fā)用于模擬、延伸和擴展人類智能的理論、方法、技術(shù)及其應用的一門新興技術(shù)科學[9,21]。實現(xiàn)單炮記錄智能評價有多種途徑,但需要結(jié)合應用領域知識才能獲得正確的解決思路。

    3.1 海量地震資料特征分析

    地震屬性研究已產(chǎn)生了豐富的研究成果,定義的地震屬性多達幾百種[22],其中振幅、能量屬性就有30種之多。海量地震資料采集單炮數(shù)據(jù)量大、炮數(shù)多,但提取的地震屬性,特別是反映單炮記錄品質(zhì)的基本屬性相對較少。表征野外單炮記錄特征的屬性可分為單炮、單道、初至前、初至區(qū)及目的層等的能量、頻率、信噪比等。多個工區(qū)的實驗分析表明,同一工區(qū)的單炮記錄屬性大多呈高斯分布或偏態(tài)分布(圖3a),這一特征為地震屬性的數(shù)據(jù)預處理和統(tǒng)計學分析提供了數(shù)學基礎。

    地震屬性體現(xiàn)單炮品質(zhì)的程度不同,其間相關性差異較大。圖3b是該工區(qū)炮集能量與目的層能量對比圖,圖中藍色折線為炮集能量分布,黃色折線為目的層能量分布,顯然,這兩類能量間具有較強的相關性。一般地,機器學習應避免相關性強的屬性參與。通過對這些屬性所反映的單炮品質(zhì)深入挖掘,發(fā)現(xiàn)其潛在聯(lián)系,進而形成分類模型。

    圖3 某工區(qū)地震屬性分布特征(a)目的層能量分布;(b)炮集能量與目的層能量曲線對比

    3.2 基于RF的單炮記錄智能評價流程

    圖4為基于RF的單炮記錄智能評價流程。使用前期勘探成果及試驗炮建立初始樣本集,提取樣本屬性并構(gòu)建協(xié)方差矩陣,求解該矩陣的特征值和特征向量,利用主成分分析法(該步可選),在剔除強相關屬性后,將對地震資料品質(zhì)更具敏感的屬性挑選出來,基于這些優(yōu)化屬性并結(jié)合廢品庫,擴增樣本以滿足機器學習樣本數(shù)量;對訓練集樣本進行RF訓練,如果沒有通過驗證,則調(diào)整參數(shù)后重新訓練;輸入生產(chǎn)炮,按RF分類,在模型評估后,如果沒有達到精準度要求,則調(diào)整參數(shù)處理后重新進行分類;如果達到精準度要求,在輸出分類結(jié)果后檢測樣本集是否完備,若需要則把生產(chǎn)炮作為學習樣本補充到樣本集中,重新組成模型訓練,否則,關閉訓練模型,陸續(xù)對生產(chǎn)炮進行自動分類。

    圖4 單炮記錄智能評價流程圖

    3.3 樣本集的建立及擴充

    樣本集是人工智能學習的基礎,必須具有大量代表性及多樣性的樣本所建立的分類模型才能防止過擬合[15,23],且各類別樣本數(shù)目盡量平衡,否則,樣本數(shù)特別少的類別分類精準度低。

    利用以往勘探中地理位置鄰近或地表、地質(zhì)條件及激發(fā)與接收因素相似的地震采集工區(qū)成果,抽取具有代表性的單炮記錄,如一、二級品記錄、低頻炮、低信噪比炮;廢炮庫是由不同激發(fā)條件和接收因素以及在各種地表地質(zhì)條件下所采集的不合格單炮組成,如斷排列炮、串感炮、噪聲炮等;試驗資料具有當前工區(qū)地震激發(fā)與接收的廣泛代表性,可將其作為一、二級品樣本以減少對前期成果樣本的過度依賴。以上這些單炮組成初始樣本集。

    在智能評價建模應用初期,初始樣本集普遍存在樣本不足的情況,需要擴增樣本,在增加樣本數(shù)量時應盡量避免由此造成的過擬合問題,擴增渠道包括對部分樣本做樣本增強技術(shù)[23-24]及納入后期正確分類的生產(chǎn)炮,擴增方法一般通過數(shù)學變換或增加噪聲方式實現(xiàn)。

    由于高斯白噪的功率譜密度服從均勻分布,幅度分布服從高斯分布,利用高斯白噪聲對原始樣本集進行重構(gòu),形成新樣本集。假設地震波有效信號為W(t),時間為t,環(huán)境噪聲為H(t),高斯白噪聲為nG(t),則重構(gòu)信號Y(t)為

    Y(t)=W(t)+H(t)+nG(t)

    (6)

    在樣本擴增時,對nG(t)做如下限定

    -E(Y)≤λE(nG)≤E(Y)

    (7)

    式中:E(Y)為重構(gòu)信號能量;E(nG)為高斯噪聲能量;λ為約束因子,|λ|≤1。

    為了更好地反映特定工區(qū)中的一些特殊勘探目標要求,在試驗炮屬性提取后,按照高斯白噪聲分布規(guī)律,在限定范圍內(nèi)由試驗炮重構(gòu)新屬性。設試驗炮數(shù)為M1,新增樣本數(shù)為M2,于是,新增樣本i(i=1,2,…,M2)的屬性值w′i為

    (8)

    式中:wj為第j個試驗炮的屬性值;nGi為高斯白噪聲因子。式(8)中的λ確定了新樣本類別,由勘探目標及工區(qū)地表地質(zhì)條件等因素確定。

    圖5為LJ工區(qū)樣本擴增10倍前后單炮能量與最小優(yōu)勢頻率對比,定義|λ|≤0.20為一級品,0.25<|λ|≤0.50為二級品,0.55<|λ|<1.00為廢炮。圖中橫軸為單炮能量的自然對數(shù)值,縱軸為最小優(yōu)勢頻率(Hz)。需要說明的是,此處出現(xiàn)了0.20~0.25、0.50~0.55的數(shù)據(jù)間隙,這是因為統(tǒng)計分析認為該范圍內(nèi)生成的新樣本存在二義性(即可能為上一級,也可能為下一級)。

    圖5 LJ工區(qū)樣本擴增前后單炮能量與最小頻率交會圖

    w′p=xi+rand(0,1)×(yk-xi)

    (9)

    將原始廢炮與二級品樣本結(jié)合(式(9)),比單獨增加新廢炮樣本方式(式(8))更能改善樣本集性能,提高分類效率。

    3.4 單炮記錄RF分類建模算法

    3.4.1 分類建模算法

    基于RF單炮記錄分類建模過程如下:

    (1)從原始樣本集G1中以樣本放回方式隨機抽取M個樣本組成訓練集G,按如下過程遞歸生成決策樹:

    1)從N個地震屬性中采用放回方式隨機抽取K(K

    2)在具有K個連續(xù)屬性的M個樣本中求最優(yōu)屬性,以最優(yōu)屬性建立分支節(jié)點;

    3)分別判斷分支節(jié)點的左右支樣本集中各樣本的類別標簽(一級品、二級品和廢炮)是否屬同一類,或是葉節(jié)點,如果不是,K=K-1,返回2);

    (2)重復步驟(1)P次,生成P棵決策樹,形成RF;

    (3)對于生產(chǎn)炮,遍歷RF的每顆決策樹,統(tǒng)計每棵樹的分類結(jié)果,取最多的類別為該炮類別。

    3.4.2 連續(xù)性屬性的分支節(jié)點構(gòu)建

    基于RF的單炮記錄分類最基礎也是最重要的工作是分支節(jié)點的建立。

    通常,決策樹采用離散值作為節(jié)點,而地震屬性是連續(xù)數(shù)據(jù)。為此,對于由M個樣本組成的訓練集G,若每個樣本只取K個屬性,則G可表示為G={(xi,1,xi,2,…,xi,K,li)|i=1,2,…,M},其中,xi,j為樣本i的第j個屬性,li為樣本炮i的類別標簽。

    為了從K個地震屬性中最優(yōu)構(gòu)建決策樹的分支節(jié)點,定義信息增益作為決策樹最優(yōu)屬性的衡量指標。地震屬性j的信息增益定義如下

    (10)

    (11)

    式中pk為第k(k=1,2,3)類(分別對應一級品、二級品和廢炮)樣本在樣本集G中的占比。

    3.4.3 RF單炮分類建模參數(shù)

    決策樹數(shù)目P決定了RF規(guī)模,也體現(xiàn)RF分類性能。理論上,P越大,分類效果越好[25],但計算量會隨之提高。通常,參考樣本擴增后的樣本數(shù)及其屬性數(shù)確定P值,如果這些數(shù)目較多,決策樹的數(shù)目可相對少一些,一般以一百到數(shù)百棵為宜。

    僅從決策樹角度看,為減少異常噪聲影響,防止過擬合,需要對決策樹剪枝處理[23,25],利用以下參數(shù)進行預剪枝:最大深度、內(nèi)部節(jié)點劃分所需最小樣本數(shù)和葉節(jié)點最小樣本數(shù),這些參數(shù)的選取和調(diào)整參數(shù)的順序與具體數(shù)據(jù)分布有關,可根據(jù)局部尋優(yōu)方法依次確定[26]。文獻[16,23,25]說明:RF中兩個隨機性(隨機樣本和隨機屬性抽取)的引入使分類算法完全可避免過擬合現(xiàn)象,況且樣本集擴增已采取了多樣性增強措施。但考慮現(xiàn)場計算能力,也可對決策樹通過預剪枝以減少計算量。目前,一些開源實用開發(fā)庫[27]已提供成熟的算法較好地優(yōu)化這些參數(shù),本文不再探討。

    3.4.4 單炮記錄RF分類結(jié)果評估

    在上述算法中,構(gòu)建所有決策樹使用了P×M個樣本,但其中包含大量相同的樣本。因此,從概率上分析,樣本集中仍有36.8%的樣本未參與訓練[28],可用作驗證樣本。

    利用單炮分類正確率(C)和廢炮識別率(R)作為驗證分類標準,其中,后者必須達到對廢炮的完全識別(100%),它們分別定義為

    (12)

    (13)

    式中:M5為驗證樣本總數(shù);Uk為驗證樣本經(jīng)RF分類后k類樣本的正確分類數(shù);M6是驗證樣本中的實際廢炮總數(shù);U3是驗證樣本經(jīng)RF分類后正確分類的廢炮數(shù)。

    3.5 模型在ZH6J工區(qū)應用及效果

    人工智能技術(shù)如今發(fā)展迅速,已建立了多種開源的機器學習平臺、系統(tǒng)、框架、工具包和類庫等。利用Scikit-learn 0.21.2開源庫[25],實現(xiàn)了單炮記錄分類算法,并根據(jù)現(xiàn)場計算機配置自動調(diào)整線程并行數(shù)以提高性能,該算法已集成到M軟件中。

    利用中國西部沙漠ZH6J工區(qū)資料進行應用測試,共提取了18種地震屬性,但沒有使用圖4流程中所提的主成分分析法做參數(shù)優(yōu)化。通過樣本增強與吸收生產(chǎn)炮,建立了4500炮的樣本集。在RF模型訓練時,以分類正確率C作為分類泛化能力的檢測依據(jù)。參數(shù)按如下順序調(diào)優(yōu):首先確定決策樹的數(shù)目P使算法穩(wěn)定;再確定決策樹的最大深度和內(nèi)部節(jié)點劃分所需最小樣本數(shù)以控制算法復雜度;然后,聯(lián)合調(diào)試內(nèi)部節(jié)點劃分所需最小樣本數(shù)和葉節(jié)點最小樣本數(shù)以增強決策樹的泛化能力;最后,獲得最大隨機屬性數(shù)K。在完成訓練后對生產(chǎn)炮自動分類,并與人工分類結(jié)果對比,表1是兩次統(tǒng)計結(jié)果。對某天采集的600炮數(shù)據(jù)自動分類,與人工分類結(jié)果相比,正確率達到97.33%,且準確識別出當天全部廢炮。在參數(shù)調(diào)優(yōu)后,對工區(qū)所有56797炮自動分類(已無廢炮),正確率達98.70%。需要指出的是,人工評價與實際分類存在一定誤差。

    表1 ZH6J工區(qū)單炮記錄分類結(jié)果統(tǒng)計表

    4 模型關系及其適應性與時效性

    多因素確定性評價模型是地震采集工程現(xiàn)場質(zhì)控的重點內(nèi)容之一,是其他評價模型不可替代的;多元屬性判別分析評價模型設計思路簡明,便于發(fā)現(xiàn)廢炮,適于實時單炮監(jiān)控。以上兩種模型相結(jié)合一般能夠及時發(fā)現(xiàn)異常道、異常排列和廢炮。智能評價模型從眾多已有標簽的樣本及其各類地震屬性中學習,分類方式客觀,可用于實時監(jiān)控單炮質(zhì)量,也可用于單炮的延時分析評估,在標準記錄選擇與閾值設置困難的勘探程度相對較低地區(qū),該模型優(yōu)勢尤其明顯。

    以上三種評價模型適用于不同地震采集方式。海上或過渡帶、多波多分量地震資料等有其獨特性,主要體現(xiàn)在地震數(shù)據(jù)記錄方式,可根據(jù)各自的特點首先進行資料預處理,然后采用上述模型分類處理。例如,雙檢單炮記錄包含了陸檢和水檢分量,需把單炮記錄解編為陸檢和水檢單炮數(shù)據(jù)結(jié)構(gòu)后[28],再分別建模分類。

    影響海量地震數(shù)據(jù)采集實時質(zhì)控的因素主要包括網(wǎng)絡傳輸速率、單炮數(shù)據(jù)解編與屬性提取、分類評價等。實驗表明,目前決定實時質(zhì)控效率的關鍵是傳輸。例如,10萬道接收的單炮若采用7s長度和1ms采樣,以SEGD格式存儲單炮數(shù)據(jù)量約為2.6GB,若采用野外較通用的千兆網(wǎng)傳輸,該單炮數(shù)據(jù)傳輸與存儲約耗時31.32s;若采用先進的光纜傳輸與高效的固態(tài)硬盤存儲大約需要5.9s。因此,應專門研究這種瞬時大數(shù)據(jù)流和持續(xù)大數(shù)據(jù)量的數(shù)據(jù)傳輸與存儲模式,以適應實時質(zhì)控需求。

    不同于數(shù)據(jù)傳輸受限于網(wǎng)絡和硬盤讀寫等物理因素,單炮數(shù)據(jù)解編與屬性提取采用內(nèi)存映射、多線程并行等綜合優(yōu)化技術(shù)后實際數(shù)據(jù)處理能力顯著提升[29],耗時主要在時間域到頻率域變換過程。就單純的單炮記錄分類模型而言,多因素確定性評價模型耗時主要在一些定量化分析上[30];多元屬性判別分析模型主要耗時在區(qū)域劃分和標準記錄的選取,生產(chǎn)炮分類時僅僅是指定屬性門檻值的比對。上述兩類模型耗時幾乎都在毫秒數(shù)量級。而智能評價模型耗時主要在分類建模階段,由于需要反復建模與驗證,一般可在采集試驗后完成,但評價模型一旦建立,實際生產(chǎn)炮的分類可在秒內(nèi)完成。

    總之,在時效性方面,單炮傳輸時間在數(shù)秒到十數(shù)秒,甚至數(shù)十秒,解編和屬性提取一般2s內(nèi)可完成,而分類過程不到1s。

    5 結(jié)論及認識

    多因素確定性評價模型與各種變形的基于多元屬性判別分析的單炮評價模型,已在地震采集工程現(xiàn)場質(zhì)量監(jiān)控中發(fā)揮了重要作用。不過,多元屬性判別分析模型所基于的標準記錄與閾值定義主觀性太強,三級判別更加困難。人工智能單炮評價方法匯集以往的勘探成果,利用試驗炮和廢炮擴充樣本,既保持了各炮的獨立性,又增加了樣本的多樣性,提高了樣本集的整體性能,彌補了不平衡樣本集可能帶來的較大分類誤差;基于RF的單炮記錄分類建模利用兩個隨機性引入,避免了人工智能最易出現(xiàn)的過擬合問題,增強了算法的穩(wěn)定性,且該模型計算過程易于高度并行化處理,評價結(jié)果客觀,適用于海量地震采集現(xiàn)場質(zhì)量監(jiān)控。

    地震采集是復雜的系統(tǒng)工程,利用人工智能對地震資料品質(zhì)評價需要綜合考慮各方面的因素,原始樣本庫建設、模型分類效果評價[31]及其物理解釋等是下一步研究方向。

    特別感謝中石化石油工程地球物理有限公司勝利分公司石翠翠女士,為本文研究提供了RF單炮分類實驗結(jié)果。

    猜你喜歡
    單炮工區(qū)分類
    地震數(shù)據(jù)常用顯示方式解析
    地震數(shù)據(jù)采集現(xiàn)場實時輸出附地質(zhì)層位單炮記錄的智能方法
    分類算一算
    關于鐵路編組站減速頂工區(qū)標準化建設研究
    淺析平橋北三維工區(qū)影響單炮品質(zhì)的因素
    分類討論求坐標
    數(shù)據(jù)分析中的分類討論
    精確發(fā)現(xiàn)溢流研究及在西北工區(qū)現(xiàn)場應用
    錄井工程(2017年3期)2018-01-22 08:40:07
    鐵路隧道定額工區(qū)長度和通風管線路工程量計算規(guī)則探討
    教你一招:數(shù)的分類
    亚洲欧美清纯卡通| 直男gayav资源| 成人精品一区二区免费| 一区二区三区激情视频| 久久婷婷人人爽人人干人人爱| 俄罗斯特黄特色一大片| 成人av在线播放网站| 十八禁国产超污无遮挡网站| 三级男女做爰猛烈吃奶摸视频| 久久人人精品亚洲av| 国产精品人妻久久久影院| 日本 欧美在线| 国内精品久久久久久久电影| 亚洲成av人片在线播放无| 精品一区二区三区人妻视频| 精品乱码久久久久久99久播| 毛片一级片免费看久久久久 | 99久久精品一区二区三区| 国产老妇女一区| 在线观看一区二区三区| 久久久久久久精品吃奶| 亚洲美女搞黄在线观看 | 美女免费视频网站| 亚洲欧美日韩卡通动漫| 亚洲男人的天堂狠狠| 亚洲人与动物交配视频| 人人妻人人看人人澡| 国产精品一区二区性色av| 我的女老师完整版在线观看| 97热精品久久久久久| 亚洲久久久久久中文字幕| 国产精品不卡视频一区二区| 免费看a级黄色片| 色综合站精品国产| 两人在一起打扑克的视频| 国产在线男女| 日韩欧美国产一区二区入口| 亚洲精品粉嫩美女一区| 欧美色视频一区免费| av在线老鸭窝| 亚洲成人久久性| 99热只有精品国产| 99riav亚洲国产免费| 国内久久婷婷六月综合欲色啪| 婷婷亚洲欧美| 日日撸夜夜添| 国产精华一区二区三区| 中文字幕人妻熟人妻熟丝袜美| 国产黄片美女视频| 国产精品免费一区二区三区在线| 久9热在线精品视频| 99热这里只有是精品50| 欧美一区二区精品小视频在线| 国产精品一区二区三区四区久久| 熟女人妻精品中文字幕| 国产欧美日韩精品一区二区| 22中文网久久字幕| 听说在线观看完整版免费高清| 午夜激情欧美在线| 内射极品少妇av片p| 人人妻人人澡欧美一区二区| 日本a在线网址| 99久国产av精品| 日韩国内少妇激情av| 自拍偷自拍亚洲精品老妇| 偷拍熟女少妇极品色| 麻豆精品久久久久久蜜桃| 18+在线观看网站| 此物有八面人人有两片| 日本爱情动作片www.在线观看 | 熟妇人妻久久中文字幕3abv| 精品久久久久久,| 深爱激情五月婷婷| 亚洲avbb在线观看| 亚洲成人免费电影在线观看| 久久久久久大精品| 国产v大片淫在线免费观看| 国产免费男女视频| 亚洲av.av天堂| 18禁裸乳无遮挡免费网站照片| 天堂影院成人在线观看| 成人av一区二区三区在线看| 欧美性感艳星| 黄色丝袜av网址大全| 日本成人三级电影网站| 一进一出好大好爽视频| 男女之事视频高清在线观看| 99久久精品国产国产毛片| 久9热在线精品视频| 国产久久久一区二区三区| 亚洲天堂国产精品一区在线| 国产色爽女视频免费观看| 美女被艹到高潮喷水动态| 91狼人影院| 亚洲中文字幕一区二区三区有码在线看| 此物有八面人人有两片| 久久婷婷人人爽人人干人人爱| 少妇的逼好多水| 一个人观看的视频www高清免费观看| 午夜精品久久久久久毛片777| 精品日产1卡2卡| 国产激情偷乱视频一区二区| 黄色日韩在线| 日本与韩国留学比较| 美女黄网站色视频| 观看美女的网站| 国产人妻一区二区三区在| 亚洲国产精品久久男人天堂| 人妻制服诱惑在线中文字幕| 午夜福利欧美成人| 国产亚洲精品久久久com| 日本三级黄在线观看| 欧美日韩亚洲国产一区二区在线观看| 精品久久久久久,| 99久久久亚洲精品蜜臀av| 中文资源天堂在线| 日韩精品中文字幕看吧| 在线观看舔阴道视频| 欧美高清成人免费视频www| 亚洲色图av天堂| 国产一级毛片七仙女欲春2| 搡老岳熟女国产| 91久久精品国产一区二区三区| videossex国产| 成人鲁丝片一二三区免费| 欧美潮喷喷水| 国产视频一区二区在线看| 国产精品无大码| 99在线人妻在线中文字幕| 亚洲黑人精品在线| 乱系列少妇在线播放| 老女人水多毛片| 白带黄色成豆腐渣| 亚洲 国产 在线| avwww免费| 日日撸夜夜添| 日本爱情动作片www.在线观看 | 性欧美人与动物交配| 日本黄色片子视频| 色视频www国产| 亚洲天堂国产精品一区在线| 亚洲,欧美,日韩| 小蜜桃在线观看免费完整版高清| 国产亚洲精品久久久久久毛片| 村上凉子中文字幕在线| 日日干狠狠操夜夜爽| 性色avwww在线观看| 老司机福利观看| 色综合婷婷激情| 两个人的视频大全免费| 黄色女人牲交| 午夜亚洲福利在线播放| 亚洲人成伊人成综合网2020| 99久久精品一区二区三区| 女同久久另类99精品国产91| 91午夜精品亚洲一区二区三区 | 久久午夜福利片| 国产成人av教育| 两人在一起打扑克的视频| av女优亚洲男人天堂| 亚洲男人的天堂狠狠| 美女被艹到高潮喷水动态| 久久久久久久久久成人| 人人妻人人看人人澡| 三级国产精品欧美在线观看| 国产主播在线观看一区二区| 色综合站精品国产| 可以在线观看的亚洲视频| 精品久久久久久久久久久久久| 中国美白少妇内射xxxbb| 日韩精品青青久久久久久| 国产国拍精品亚洲av在线观看| 日本一二三区视频观看| 欧美日韩亚洲国产一区二区在线观看| 国产精品亚洲美女久久久| 亚洲专区中文字幕在线| www.色视频.com| 国产欧美日韩一区二区精品| 91久久精品国产一区二区三区| 午夜免费激情av| 亚洲avbb在线观看| 99热这里只有是精品50| 国产伦一二天堂av在线观看| 少妇熟女aⅴ在线视频| 午夜福利高清视频| 国产蜜桃级精品一区二区三区| 国产男人的电影天堂91| 一个人免费在线观看电影| 亚洲成人久久爱视频| 美女被艹到高潮喷水动态| 22中文网久久字幕| 亚洲av成人av| 亚洲欧美日韩高清专用| 久久久久久久亚洲中文字幕| 网址你懂的国产日韩在线| 成人av一区二区三区在线看| 男人舔奶头视频| 欧美又色又爽又黄视频| 国产精品久久久久久精品电影| 少妇的逼好多水| 日本欧美国产在线视频| 日本三级黄在线观看| 国产黄片美女视频| 日本 av在线| 色精品久久人妻99蜜桃| 亚洲久久久久久中文字幕| 少妇高潮的动态图| 国产不卡一卡二| 丝袜美腿在线中文| 校园春色视频在线观看| 国产黄片美女视频| 色综合亚洲欧美另类图片| 最后的刺客免费高清国语| 亚洲av不卡在线观看| 亚洲欧美日韩东京热| 黄片wwwwww| 欧美性猛交黑人性爽| 免费搜索国产男女视频| 亚洲一区高清亚洲精品| 大又大粗又爽又黄少妇毛片口| 欧美不卡视频在线免费观看| 午夜免费成人在线视频| 国产欧美日韩精品亚洲av| 色播亚洲综合网| 午夜福利成人在线免费观看| 亚洲中文字幕一区二区三区有码在线看| 日韩欧美三级三区| 18禁裸乳无遮挡免费网站照片| av在线老鸭窝| 悠悠久久av| 五月玫瑰六月丁香| 成人性生交大片免费视频hd| 国产探花在线观看一区二区| 91午夜精品亚洲一区二区三区 | 日本三级黄在线观看| 亚洲经典国产精华液单| 色噜噜av男人的天堂激情| 男女视频在线观看网站免费| 成人永久免费在线观看视频| 午夜福利18| av在线天堂中文字幕| 国产精品久久视频播放| 色视频www国产| 亚洲在线自拍视频| 国产白丝娇喘喷水9色精品| 亚洲乱码一区二区免费版| a级毛片a级免费在线| 日本免费a在线| 性插视频无遮挡在线免费观看| 18禁黄网站禁片午夜丰满| 久久99热6这里只有精品| 日本三级黄在线观看| 亚洲色图av天堂| 亚洲综合色惰| 欧美色视频一区免费| 人人妻人人看人人澡| 美女大奶头视频| 69av精品久久久久久| 内地一区二区视频在线| 九九在线视频观看精品| 亚洲18禁久久av| 国产淫片久久久久久久久| 99久久精品一区二区三区| 91在线精品国自产拍蜜月| 亚洲最大成人手机在线| 十八禁国产超污无遮挡网站| 亚洲无线观看免费| 又粗又爽又猛毛片免费看| 色吧在线观看| 色播亚洲综合网| 成熟少妇高潮喷水视频| 日本熟妇午夜| 18禁黄网站禁片免费观看直播| 亚洲18禁久久av| 麻豆成人午夜福利视频| 亚洲精华国产精华液的使用体验 | 中出人妻视频一区二区| 亚洲最大成人中文| 99久久精品热视频| www.www免费av| 可以在线观看毛片的网站| 亚洲美女黄片视频| 99久久无色码亚洲精品果冻| 男女那种视频在线观看| 国产高潮美女av| 男女之事视频高清在线观看| aaaaa片日本免费| 日韩,欧美,国产一区二区三区 | 啦啦啦观看免费观看视频高清| 久久精品国产亚洲av香蕉五月| 成人av在线播放网站| 日日摸夜夜添夜夜添av毛片 | 久久久久久久精品吃奶| 在现免费观看毛片| 亚洲在线观看片| 日日摸夜夜添夜夜添av毛片 | 美女高潮的动态| 亚洲国产日韩欧美精品在线观看| 老熟妇仑乱视频hdxx| 琪琪午夜伦伦电影理论片6080| 欧美中文日本在线观看视频| 精品福利观看| 亚洲自偷自拍三级| 免费在线观看成人毛片| 国产毛片a区久久久久| 亚洲av美国av| 国产精品美女特级片免费视频播放器| 色5月婷婷丁香| 99久久精品热视频| 国产精品一区www在线观看 | 国产精品女同一区二区软件 | 久久国产乱子免费精品| 欧美丝袜亚洲另类 | 国产69精品久久久久777片| 欧美一区二区精品小视频在线| 禁无遮挡网站| 999久久久精品免费观看国产| 精品久久久久久久末码| 三级男女做爰猛烈吃奶摸视频| 欧美激情在线99| 伦精品一区二区三区| 男女边吃奶边做爰视频| 国产欧美日韩精品一区二区| 成年女人永久免费观看视频| 看片在线看免费视频| 久久人人爽人人爽人人片va| 精品久久久噜噜| 久久中文看片网| 精品99又大又爽又粗少妇毛片 | 观看免费一级毛片| 亚洲精品在线观看二区| 亚洲欧美日韩卡通动漫| 国产日本99.免费观看| 亚洲无线观看免费| 美女免费视频网站| 久久午夜福利片| 成人鲁丝片一二三区免费| 91狼人影院| 18禁黄网站禁片免费观看直播| 国产一区二区激情短视频| 久久午夜亚洲精品久久| 真人做人爱边吃奶动态| 国产真实伦视频高清在线观看 | 国产精品98久久久久久宅男小说| 国产欧美日韩精品亚洲av| 老司机午夜福利在线观看视频| 动漫黄色视频在线观看| 一级毛片久久久久久久久女| 成年女人毛片免费观看观看9| 成人一区二区视频在线观看| 亚洲欧美日韩无卡精品| 一本一本综合久久| 欧美激情久久久久久爽电影| а√天堂www在线а√下载| 久久精品国产亚洲av涩爱 | 欧美精品国产亚洲| 成人综合一区亚洲| 亚洲国产欧洲综合997久久,| av专区在线播放| 麻豆国产97在线/欧美| 成人一区二区视频在线观看| 欧美性猛交黑人性爽| 国产在线男女| 亚洲精品国产成人久久av| 中文字幕高清在线视频| 综合色av麻豆| 欧美激情国产日韩精品一区| 嫩草影院新地址| 国产伦一二天堂av在线观看| 噜噜噜噜噜久久久久久91| 九色国产91popny在线| av在线观看视频网站免费| av在线天堂中文字幕| 精品99又大又爽又粗少妇毛片 | 精品久久国产蜜桃| 69人妻影院| 国产精品人妻久久久久久| 国产 一区精品| 亚洲精品影视一区二区三区av| 久久精品国产鲁丝片午夜精品 | 性色avwww在线观看| 国产精品野战在线观看| 亚洲乱码一区二区免费版| 久久久色成人| 国产精品1区2区在线观看.| 偷拍熟女少妇极品色| 国产在线精品亚洲第一网站| 婷婷色综合大香蕉| bbb黄色大片| 九九在线视频观看精品| 欧美日韩乱码在线| 免费在线观看日本一区| 男女做爰动态图高潮gif福利片| aaaaa片日本免费| 成人特级黄色片久久久久久久| 国内毛片毛片毛片毛片毛片| 中文字幕久久专区| 国产主播在线观看一区二区| 非洲黑人性xxxx精品又粗又长| 国产精品1区2区在线观看.| 亚洲精品色激情综合| 1024手机看黄色片| 最近视频中文字幕2019在线8| 午夜激情欧美在线| 看片在线看免费视频| 欧美日本亚洲视频在线播放| 一个人观看的视频www高清免费观看| 亚洲美女搞黄在线观看 | av在线蜜桃| 最近最新免费中文字幕在线| 看黄色毛片网站| 美女高潮喷水抽搐中文字幕| 校园春色视频在线观看| 亚洲精品在线观看二区| 免费观看的影片在线观看| 老司机深夜福利视频在线观看| 国产男人的电影天堂91| 亚洲图色成人| 国产午夜精品久久久久久一区二区三区 | 内射极品少妇av片p| 色综合站精品国产| 国产免费男女视频| 午夜福利在线观看吧| 日本-黄色视频高清免费观看| 久久久色成人| 熟女电影av网| 日韩欧美免费精品| 亚洲久久久久久中文字幕| 亚洲成人久久性| 免费大片18禁| 国产 一区 欧美 日韩| 亚洲最大成人中文| 国产精品一区二区三区四区久久| 男人狂女人下面高潮的视频| 麻豆成人av在线观看| 国产精品一区二区三区四区久久| 中文字幕免费在线视频6| 欧美日本亚洲视频在线播放| 人妻少妇偷人精品九色| 亚洲人成网站在线播放欧美日韩| 日本a在线网址| 免费av观看视频| 免费电影在线观看免费观看| 欧美另类亚洲清纯唯美| 色在线成人网| 久9热在线精品视频| 99视频精品全部免费 在线| 女人十人毛片免费观看3o分钟| 欧美日韩国产亚洲二区| 久久久久久大精品| 联通29元200g的流量卡| 国产精品免费一区二区三区在线| 国产久久久一区二区三区| 美女大奶头视频| 午夜激情欧美在线| 精品不卡国产一区二区三区| 哪里可以看免费的av片| 国产高清三级在线| 露出奶头的视频| 欧美绝顶高潮抽搐喷水| 在线看三级毛片| netflix在线观看网站| 少妇的逼好多水| 99热这里只有是精品50| 一进一出好大好爽视频| 国产精品久久电影中文字幕| 国产精品无大码| 久久99热6这里只有精品| 亚洲精品亚洲一区二区| 在线免费十八禁| 日本黄大片高清| 中亚洲国语对白在线视频| 伊人久久精品亚洲午夜| 欧美一区二区亚洲| 91在线观看av| 丰满乱子伦码专区| 麻豆一二三区av精品| 亚洲国产欧洲综合997久久,| 又黄又爽又刺激的免费视频.| 色吧在线观看| 日韩欧美在线乱码| 亚洲精品影视一区二区三区av| 色哟哟哟哟哟哟| 国内精品久久久久久久电影| 国产不卡一卡二| 久久久精品大字幕| 在线观看美女被高潮喷水网站| av天堂中文字幕网| 国产麻豆成人av免费视频| 亚洲成人精品中文字幕电影| 国产综合懂色| 听说在线观看完整版免费高清| 伦理电影大哥的女人| videossex国产| 亚洲人成网站在线播放欧美日韩| 色精品久久人妻99蜜桃| 黄色配什么色好看| 我要看日韩黄色一级片| 成人欧美大片| 久久中文看片网| 九九爱精品视频在线观看| 最好的美女福利视频网| 三级毛片av免费| 国产高潮美女av| 亚洲欧美日韩高清专用| 悠悠久久av| 精品一区二区三区视频在线| 成人特级黄色片久久久久久久| 中亚洲国语对白在线视频| 久久久久国产精品人妻aⅴ院| 岛国在线免费视频观看| 伦精品一区二区三区| 日本免费a在线| 免费看光身美女| av视频在线观看入口| 日韩精品中文字幕看吧| 亚洲国产精品成人综合色| 欧美日韩瑟瑟在线播放| 欧美激情在线99| 亚洲av二区三区四区| 午夜福利在线观看吧| 三级男女做爰猛烈吃奶摸视频| 久久午夜福利片| 在线天堂最新版资源| 毛片女人毛片| 国产高清视频在线播放一区| 在线免费十八禁| 男女啪啪激烈高潮av片| 热99re8久久精品国产| 日日啪夜夜撸| 最新在线观看一区二区三区| 精华霜和精华液先用哪个| 亚洲av成人精品一区久久| 少妇裸体淫交视频免费看高清| 全区人妻精品视频| 又爽又黄无遮挡网站| 日本精品一区二区三区蜜桃| 免费看光身美女| 精品一区二区三区视频在线| 少妇丰满av| 99精品在免费线老司机午夜| 亚洲黑人精品在线| 日本黄大片高清| 女的被弄到高潮叫床怎么办 | 真人做人爱边吃奶动态| 午夜福利在线观看吧| 日本五十路高清| 观看免费一级毛片| 观看美女的网站| 精品一区二区三区人妻视频| 色av中文字幕| 男女边吃奶边做爰视频| av天堂中文字幕网| 夜夜夜夜夜久久久久| 亚洲精品色激情综合| 天天一区二区日本电影三级| 欧美色视频一区免费| 五月伊人婷婷丁香| 窝窝影院91人妻| 嫩草影院入口| 波多野结衣高清无吗| 国产综合懂色| 嫩草影院新地址| 日本黄大片高清| 亚洲人成网站在线播| 国产一区二区激情短视频| 伦理电影大哥的女人| 国产中年淑女户外野战色| 色吧在线观看| 天天一区二区日本电影三级| 可以在线观看的亚洲视频| 日韩欧美一区二区三区在线观看| 91狼人影院| 国产欧美日韩一区二区精品| 国产主播在线观看一区二区| 99久久无色码亚洲精品果冻| 亚洲av二区三区四区| 两人在一起打扑克的视频| 亚洲精品色激情综合| 国产精品一区二区三区四区免费观看 | 亚洲国产色片| 男人舔女人下体高潮全视频| 午夜免费男女啪啪视频观看 | 国产 一区 欧美 日韩| 日本爱情动作片www.在线观看 | 国产精品三级大全| 麻豆成人av在线观看| 成人二区视频| 男女边吃奶边做爰视频| 国产视频一区二区在线看| 免费看日本二区| 色综合站精品国产| 夜夜看夜夜爽夜夜摸| 色av中文字幕| 中文字幕熟女人妻在线| 欧美日韩黄片免| 九九爱精品视频在线观看| 黄色视频,在线免费观看| 麻豆久久精品国产亚洲av| 日本三级黄在线观看| 日韩欧美免费精品| 22中文网久久字幕| 亚洲第一电影网av| 一个人免费在线观看电影| 在线观看av片永久免费下载| 老司机福利观看| 亚洲最大成人av| 国产精品一区二区性色av| 国产精品久久久久久久久免| 亚洲av日韩精品久久久久久密| 国产亚洲精品久久久久久毛片| 欧美成人免费av一区二区三区| 国国产精品蜜臀av免费| 亚洲成人免费电影在线观看| 欧美另类亚洲清纯唯美|