• <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ù)的分類
    国产精品女同一区二区软件| 婷婷成人精品国产| 大又大粗又爽又黄少妇毛片口| 国产成人freesex在线| 欧美三级亚洲精品| 亚洲美女视频黄频| 人人妻人人添人人爽欧美一区卜| 少妇人妻精品综合一区二区| 黄色欧美视频在线观看| 免费观看av网站的网址| xxx大片免费视频| 一二三四中文在线观看免费高清| 晚上一个人看的免费电影| 成人午夜精彩视频在线观看| av.在线天堂| 18在线观看网站| 菩萨蛮人人尽说江南好唐韦庄| 一边亲一边摸免费视频| 国产精品熟女久久久久浪| 自线自在国产av| 日韩在线高清观看一区二区三区| 精品人妻偷拍中文字幕| 亚洲精品色激情综合| 日日摸夜夜添夜夜爱| 麻豆乱淫一区二区| 人妻 亚洲 视频| 最近中文字幕2019免费版| 成人黄色视频免费在线看| av免费在线看不卡| 好男人视频免费观看在线| 日本黄色片子视频| 熟女av电影| 精品熟女少妇av免费看| 日本wwww免费看| 日本爱情动作片www.在线观看| a 毛片基地| 亚洲人成网站在线播| 一级爰片在线观看| 亚洲国产毛片av蜜桃av| av在线app专区| 91成人精品电影| 十八禁高潮呻吟视频| 日韩制服骚丝袜av| 少妇被粗大猛烈的视频| 精品一区二区免费观看| 国产男人的电影天堂91| 久久久精品94久久精品| 亚洲,一卡二卡三卡| 精品国产国语对白av| 日本wwww免费看| av免费在线看不卡| 男女边吃奶边做爰视频| 18禁裸乳无遮挡动漫免费视频| 亚洲av成人精品一二三区| 五月伊人婷婷丁香| 中文欧美无线码| 欧美亚洲 丝袜 人妻 在线| 亚洲精品自拍成人| 汤姆久久久久久久影院中文字幕| 成人毛片60女人毛片免费| 欧美精品一区二区大全| 男女啪啪激烈高潮av片| 午夜福利网站1000一区二区三区| 欧美三级亚洲精品| 国产免费现黄频在线看| 性高湖久久久久久久久免费观看| 欧美激情国产日韩精品一区| 国产精品99久久久久久久久| 十八禁网站网址无遮挡| 黄色配什么色好看| 精品国产露脸久久av麻豆| 亚洲天堂av无毛| 菩萨蛮人人尽说江南好唐韦庄| 少妇人妻久久综合中文| 精品国产一区二区久久| 2018国产大陆天天弄谢| av卡一久久| 高清不卡的av网站| 久久精品久久精品一区二区三区| 欧美日韩国产mv在线观看视频| 亚洲成人手机| 简卡轻食公司| 亚洲av在线观看美女高潮| 欧美成人精品欧美一级黄| 最近中文字幕高清免费大全6| 欧美性感艳星| 午夜老司机福利剧场| 99久久精品国产国产毛片| 国产视频内射| 热re99久久精品国产66热6| 成人二区视频| 中文字幕制服av| 伊人久久国产一区二区| 乱码一卡2卡4卡精品| 日本黄色片子视频| 特大巨黑吊av在线直播| 成年av动漫网址| 在线观看一区二区三区激情| 麻豆成人av视频| www.av在线官网国产| 两个人免费观看高清视频| 国产伦理片在线播放av一区| 国产在线免费精品| 插逼视频在线观看| a级毛片在线看网站| 国产精品国产三级专区第一集| 日韩中字成人| 亚洲av福利一区| 日韩欧美精品免费久久| 久久精品久久久久久久性| 涩涩av久久男人的天堂| 欧美变态另类bdsm刘玥| 大香蕉久久网| 国产免费福利视频在线观看| 亚洲精品日韩av片在线观看| 99久久中文字幕三级久久日本| 97在线视频观看| 婷婷色综合www| 婷婷色av中文字幕| 狠狠精品人妻久久久久久综合| 久久久久久久大尺度免费视频| 亚洲国产成人一精品久久久| 国产在线一区二区三区精| 中文字幕制服av| 成年av动漫网址| 免费黄色在线免费观看| 国产精品一二三区在线看| 国产无遮挡羞羞视频在线观看| 国产精品久久久久久久久免| 大又大粗又爽又黄少妇毛片口| 日韩熟女老妇一区二区性免费视频| 欧美精品人与动牲交sv欧美| 69精品国产乱码久久久| 七月丁香在线播放| 女人精品久久久久毛片| 日韩制服骚丝袜av| 日本爱情动作片www.在线观看| 久久久精品区二区三区| 久久久精品免费免费高清| 国产免费一级a男人的天堂| 日本与韩国留学比较| 日韩av免费高清视频| 最近中文字幕高清免费大全6| 亚洲欧美日韩另类电影网站| 蜜桃久久精品国产亚洲av| 精品久久久噜噜| 国产伦理片在线播放av一区| 毛片一级片免费看久久久久| 国产精品久久久久成人av| 婷婷色av中文字幕| 男人添女人高潮全过程视频| 亚洲无线观看免费| 亚洲国产精品一区二区三区在线| 少妇人妻 视频| 国产精品欧美亚洲77777| xxxhd国产人妻xxx| 99九九线精品视频在线观看视频| 国产精品熟女久久久久浪| 涩涩av久久男人的天堂| 国产精品99久久久久久久久| 天天操日日干夜夜撸| av播播在线观看一区| 亚洲精品色激情综合| 国产精品女同一区二区软件| 一本久久精品| 亚洲精品456在线播放app| 少妇精品久久久久久久| 日韩成人伦理影院| 亚洲第一区二区三区不卡| 看非洲黑人一级黄片| 亚洲欧美成人综合另类久久久| 丰满乱子伦码专区| 亚洲欧美清纯卡通| 欧美日韩视频高清一区二区三区二| 欧美精品国产亚洲| 国产欧美日韩综合在线一区二区| 超色免费av| 欧美精品一区二区大全| 免费av不卡在线播放| 欧美97在线视频| 久久久久久久大尺度免费视频| 久久久久国产精品人妻一区二区| av有码第一页| 麻豆精品久久久久久蜜桃| 人人澡人人妻人| 欧美精品国产亚洲| 有码 亚洲区| 久久久久视频综合| av在线观看视频网站免费| 国产高清不卡午夜福利| 少妇精品久久久久久久| 久久人人爽人人片av| 能在线免费看毛片的网站| 国产男人的电影天堂91| 亚洲精品色激情综合| 国产高清有码在线观看视频| 国产成人免费无遮挡视频| 亚洲av免费高清在线观看| 女的被弄到高潮叫床怎么办| 夜夜爽夜夜爽视频| 免费高清在线观看日韩| 国产欧美日韩综合在线一区二区| 少妇丰满av| √禁漫天堂资源中文www| 黄色毛片三级朝国网站| 嘟嘟电影网在线观看| 男人添女人高潮全过程视频| www.av在线官网国产| 欧美精品一区二区大全| 国产极品粉嫩免费观看在线 | 极品人妻少妇av视频| 亚洲精品乱久久久久久| 国产精品人妻久久久影院| 草草在线视频免费看| 国产日韩一区二区三区精品不卡 | 秋霞在线观看毛片| 乱码一卡2卡4卡精品| 国产在线视频一区二区| 免费高清在线观看视频在线观看| 97超视频在线观看视频| 亚洲成人手机| 如何舔出高潮| kizo精华| 三级国产精品欧美在线观看| 中国三级夫妇交换| 男男h啪啪无遮挡| 亚洲精品456在线播放app| 97超碰精品成人国产| 亚洲av日韩在线播放| xxx大片免费视频| 我要看黄色一级片免费的| 在线观看人妻少妇| 高清午夜精品一区二区三区| 国产成人av激情在线播放 | 欧美精品一区二区免费开放| 中文字幕精品免费在线观看视频 | 日产精品乱码卡一卡2卡三| 中文字幕av电影在线播放| 久久久久久久亚洲中文字幕| 3wmmmm亚洲av在线观看| 最近2019中文字幕mv第一页| 亚洲成人av在线免费| 日韩强制内射视频| 极品少妇高潮喷水抽搐| 哪个播放器可以免费观看大片| 国产精品久久久久成人av| 热99久久久久精品小说推荐| 女人精品久久久久毛片| 飞空精品影院首页| 2022亚洲国产成人精品| 99久久精品国产国产毛片| 高清视频免费观看一区二区| 一本—道久久a久久精品蜜桃钙片| 秋霞在线观看毛片| 国产免费福利视频在线观看| 18禁在线无遮挡免费观看视频| 久热这里只有精品99| 国产精品一区www在线观看| 99re6热这里在线精品视频| 五月玫瑰六月丁香| 嘟嘟电影网在线观看| 国产精品久久久久久久电影| 综合色丁香网| 久热这里只有精品99| 最新中文字幕久久久久| 女性被躁到高潮视频| 午夜影院在线不卡| 看非洲黑人一级黄片| 国产精品久久久久久精品电影小说| 国产探花极品一区二区| 国国产精品蜜臀av免费| 十八禁高潮呻吟视频| 人成视频在线观看免费观看| 一区二区三区免费毛片| 91精品三级在线观看| 亚洲一级一片aⅴ在线观看| 国产精品一二三区在线看| 天堂中文最新版在线下载| 成人18禁高潮啪啪吃奶动态图 | 国产精品一国产av| a级毛片黄视频| 日本黄色日本黄色录像| 97在线视频观看| 欧美成人午夜免费资源| 只有这里有精品99| 人妻少妇偷人精品九色| 亚洲欧美一区二区三区黑人 | 久久久精品94久久精品| 日本91视频免费播放| 久久99精品国语久久久| 国产精品久久久久久精品电影小说| videosex国产| 中国国产av一级| 老女人水多毛片| 亚洲熟女精品中文字幕| 97精品久久久久久久久久精品| 女人久久www免费人成看片| 欧美另类一区| 免费黄频网站在线观看国产| 最黄视频免费看| 黄色视频在线播放观看不卡| 狠狠婷婷综合久久久久久88av| 校园人妻丝袜中文字幕| 国产精品国产av在线观看| 精品一区二区免费观看| 中文字幕亚洲精品专区| 亚洲精品久久久久久婷婷小说| 亚洲婷婷狠狠爱综合网| www.av在线官网国产| av一本久久久久| 丁香六月天网| 亚洲国产精品专区欧美| 亚洲不卡免费看| 一区二区三区精品91| 国语对白做爰xxxⅹ性视频网站| 亚洲,欧美,日韩| 777米奇影视久久| 另类精品久久| 视频在线观看一区二区三区| 一本一本综合久久| 99国产综合亚洲精品| 综合色丁香网| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 亚洲欧美成人精品一区二区| 免费av中文字幕在线| 国产亚洲av片在线观看秒播厂| 中文欧美无线码| 水蜜桃什么品种好| 大话2 男鬼变身卡| 午夜福利影视在线免费观看| 精品国产一区二区三区久久久樱花| 午夜免费鲁丝| 中国美白少妇内射xxxbb| 午夜激情久久久久久久| 91成人精品电影| 免费播放大片免费观看视频在线观看| 欧美亚洲日本最大视频资源| 一区二区三区精品91| a级毛片黄视频| 久久国产精品大桥未久av| 一区二区三区四区激情视频| 精品人妻在线不人妻| 亚洲国产欧美在线一区| 亚洲不卡免费看| 欧美日韩在线观看h| 亚洲国产精品成人久久小说| 国产免费一级a男人的天堂| 免费观看的影片在线观看| 亚洲av二区三区四区| 免费观看性生交大片5| 国产午夜精品一二区理论片| 午夜激情久久久久久久| 午夜影院在线不卡| 国产成人精品婷婷| 丁香六月天网| 国产精品国产三级国产av玫瑰| 欧美 日韩 精品 国产| 最新的欧美精品一区二区| 王馨瑶露胸无遮挡在线观看| videosex国产| 卡戴珊不雅视频在线播放| 欧美+日韩+精品| 欧美少妇被猛烈插入视频| 亚洲av成人精品一区久久| 亚洲一级一片aⅴ在线观看| av黄色大香蕉| 一级毛片 在线播放| 久久精品国产亚洲av天美| 国产精品成人在线| 国产视频内射| 欧美日韩成人在线一区二区| 曰老女人黄片| 在线看a的网站| 精品卡一卡二卡四卡免费| 久久鲁丝午夜福利片| 亚洲一级一片aⅴ在线观看| 国产黄色视频一区二区在线观看| 国产欧美亚洲国产| 欧美少妇被猛烈插入视频| 亚洲av综合色区一区| 亚洲人成77777在线视频| 男女边吃奶边做爰视频| 久久久久久久久久久丰满| 午夜日本视频在线| av国产久精品久网站免费入址| 人妻 亚洲 视频| 亚洲欧洲国产日韩| 啦啦啦啦在线视频资源| 男的添女的下面高潮视频| 在线观看免费高清a一片| 精品视频人人做人人爽| 亚洲av.av天堂| 国产日韩欧美视频二区| 人妻夜夜爽99麻豆av| 久久国产亚洲av麻豆专区| h视频一区二区三区| 99热国产这里只有精品6| 97精品久久久久久久久久精品| 亚洲成人av在线免费| 丰满迷人的少妇在线观看| 精品亚洲成国产av| 免费播放大片免费观看视频在线观看| 久久精品久久久久久久性| 亚洲欧美精品自产自拍| 精品国产国语对白av| 亚洲,欧美,日韩| 日本猛色少妇xxxxx猛交久久| 久久精品国产a三级三级三级| 国产伦精品一区二区三区视频9| 新久久久久国产一级毛片| 日韩三级伦理在线观看| 一级毛片aaaaaa免费看小| 国产精品免费大片| 色5月婷婷丁香| av国产精品久久久久影院| 日韩熟女老妇一区二区性免费视频| 精品国产国语对白av| 免费观看av网站的网址| 国产黄色视频一区二区在线观看| av播播在线观看一区| 久久久久久久大尺度免费视频| 精品国产一区二区久久| 免费少妇av软件| 国产亚洲午夜精品一区二区久久| 国产成人一区二区在线| 久久韩国三级中文字幕| 熟女电影av网| av女优亚洲男人天堂| 交换朋友夫妻互换小说| 性色av一级| av天堂久久9| 亚洲无线观看免费| 毛片一级片免费看久久久久| 大香蕉久久网| 久久午夜福利片| 一级毛片黄色毛片免费观看视频| 91在线精品国自产拍蜜月| 在线精品无人区一区二区三| 国产精品人妻久久久影院| 久久精品久久精品一区二区三区| 两个人免费观看高清视频| 久久国产亚洲av麻豆专区| 成人国产麻豆网| 国产成人午夜福利电影在线观看| 麻豆乱淫一区二区| 十八禁高潮呻吟视频| 久久精品久久久久久久性| 999精品在线视频| 国产一区有黄有色的免费视频| 老女人水多毛片| 中文字幕制服av| av不卡在线播放| 久久久久国产精品人妻一区二区| 夜夜爽夜夜爽视频| 成年人免费黄色播放视频| 中文字幕亚洲精品专区| 久久精品国产亚洲网站| 人妻 亚洲 视频| 国产高清不卡午夜福利| 天天躁夜夜躁狠狠久久av| 777米奇影视久久| 国产高清有码在线观看视频| 丝袜脚勾引网站| 久久国产精品男人的天堂亚洲 | 欧美老熟妇乱子伦牲交| www.av在线官网国产| 永久免费av网站大全| 亚洲激情五月婷婷啪啪| 国产一区二区在线观看av| 又大又黄又爽视频免费| 2022亚洲国产成人精品| 男人操女人黄网站| 国产欧美另类精品又又久久亚洲欧美| 丰满迷人的少妇在线观看| 一级a做视频免费观看| 国产极品粉嫩免费观看在线 | 国精品久久久久久国模美| 免费观看性生交大片5| 男女边吃奶边做爰视频| 一级片'在线观看视频| 亚洲精品日韩av片在线观看| 久久99蜜桃精品久久| 制服诱惑二区| 国产 一区精品| 久久精品国产自在天天线| 亚州av有码| 成人国语在线视频| 亚洲精品乱久久久久久| 一本—道久久a久久精品蜜桃钙片| 亚洲av男天堂| 中文字幕最新亚洲高清| √禁漫天堂资源中文www| 91久久精品国产一区二区三区| 中国国产av一级| 精品一区二区免费观看| 最后的刺客免费高清国语| 日本av免费视频播放| 亚洲国产av影院在线观看| 一本大道久久a久久精品| 欧美变态另类bdsm刘玥| 只有这里有精品99| videossex国产| 狠狠精品人妻久久久久久综合| 内地一区二区视频在线| 99国产综合亚洲精品| 少妇丰满av| 91精品国产国语对白视频| 丝瓜视频免费看黄片| 国产男女超爽视频在线观看| 午夜激情久久久久久久| 人妻少妇偷人精品九色| 日韩成人伦理影院| freevideosex欧美| av国产久精品久网站免费入址| 91精品国产九色| 看十八女毛片水多多多| 日韩大片免费观看网站| 欧美精品亚洲一区二区| 久久久久久久久久人人人人人人| 色哟哟·www| 免费黄网站久久成人精品| 午夜日本视频在线| 2021少妇久久久久久久久久久| 国产一区二区三区综合在线观看 | 亚洲一区二区三区欧美精品| 最后的刺客免费高清国语| 80岁老熟妇乱子伦牲交| 成人国语在线视频| 欧美 日韩 精品 国产| 69精品国产乱码久久久| 少妇猛男粗大的猛烈进出视频| 国产欧美日韩一区二区三区在线 | 母亲3免费完整高清在线观看 | 国产午夜精品一二区理论片| 日韩亚洲欧美综合| 国产乱来视频区| av电影中文网址| 爱豆传媒免费全集在线观看| 国产男女内射视频| 亚洲婷婷狠狠爱综合网| 熟女人妻精品中文字幕| 日韩中文字幕视频在线看片| 久久综合国产亚洲精品| 日韩免费高清中文字幕av| 热re99久久精品国产66热6| 国产一区有黄有色的免费视频| 亚洲精品美女久久av网站| 欧美日韩一区二区视频在线观看视频在线| 日韩中字成人| 丝袜美足系列| 大又大粗又爽又黄少妇毛片口| 久久久精品国产亚洲av高清涩受| 欧美另类亚洲清纯唯美| 熟女少妇亚洲综合色aaa.| 精品国产一区二区三区久久久樱花| 久久中文字幕人妻熟女| 无人区码免费观看不卡 | 久久婷婷成人综合色麻豆| 免费在线观看影片大全网站| 视频区欧美日本亚洲| 99久久精品国产亚洲精品| 色播在线永久视频| 午夜视频精品福利| 午夜福利乱码中文字幕| 在线观看免费日韩欧美大片| a级毛片在线看网站| 欧美 亚洲 国产 日韩一| 午夜福利欧美成人| aaaaa片日本免费| 亚洲熟女毛片儿| 一区二区三区精品91| 极品少妇高潮喷水抽搐| 两人在一起打扑克的视频| 91字幕亚洲| 超色免费av| 日韩一区二区三区影片| 桃花免费在线播放| 日韩成人在线观看一区二区三区| 我的亚洲天堂| 亚洲精品美女久久久久99蜜臀| 欧美精品av麻豆av| 午夜福利免费观看在线| 国产av国产精品国产| 日本黄色日本黄色录像| 五月天丁香电影| 国产伦理片在线播放av一区| 黄色毛片三级朝国网站| 大香蕉久久网| 免费一级毛片在线播放高清视频 | 亚洲成a人片在线一区二区| 精品国产亚洲在线| 亚洲五月婷婷丁香| 一边摸一边做爽爽视频免费| 欧美精品人与动牲交sv欧美| 飞空精品影院首页| tocl精华| 欧美成狂野欧美在线观看| 激情视频va一区二区三区| 肉色欧美久久久久久久蜜桃| 国产一区二区三区视频了| 波多野结衣一区麻豆| 老熟妇仑乱视频hdxx| 精品国产一区二区久久| 欧美大码av| 中文字幕高清在线视频| 日韩免费高清中文字幕av| 成年动漫av网址| 亚洲精品中文字幕一二三四区 | 国产有黄有色有爽视频| 国产区一区二久久| 精品久久久精品久久久| 在线永久观看黄色视频| 欧美日韩一级在线毛片| 久久久久国产一级毛片高清牌| 黄网站色视频无遮挡免费观看|