一種壓性斷裂帶結(jié)構(gòu)測(cè)井識(shí)別方法
【專利摘要】本發(fā)明涉及一種壓性斷裂帶結(jié)構(gòu)測(cè)井識(shí)別方法，對(duì)壓性應(yīng)力環(huán)境進(jìn)行分析，選取目標(biāo)鉆井和測(cè)井曲線，計(jì)算裂縫敏感性參數(shù)，對(duì)各個(gè)參數(shù)進(jìn)行歸一化計(jì)算，構(gòu)建斷裂帶內(nèi)部結(jié)構(gòu)綜合判別參數(shù)FZI。本發(fā)明可以有效的利用測(cè)井資料對(duì)斷裂帶內(nèi)部結(jié)構(gòu)進(jìn)行識(shí)別，并且選取的曲線類型為常規(guī)的標(biāo)準(zhǔn)測(cè)井，具有很強(qiáng)的可操作性。
【專利說(shuō)明】
一種壓性斷裂帶結(jié)構(gòu)測(cè)井識(shí)別方法
技術(shù)領(lǐng)域
[0001] 本發(fā)明屬于石油天然氣勘探與開發(fā)技術(shù)領(lǐng)域，具體涉及到一種壓性斷裂帶結(jié)構(gòu)測(cè) 井識(shí)別方法。
【背景技術(shù)】
[0002] 含油氣區(qū)斷裂發(fā)育，且斷裂控制了油氣的運(yùn)移與聚集。斷裂是巖層在應(yīng)力作用下 發(fā)生破裂并沿破裂面產(chǎn)生明顯相對(duì)移動(dòng)的一種構(gòu)造變形現(xiàn)象。破裂面兩側(cè)巖層在錯(cuò)動(dòng)過(guò)程 中發(fā)生擠壓研磨，巖石破碎強(qiáng)烈，同時(shí)在其兩側(cè)產(chǎn)生大量與斷層伴生的裂縫（宋到福等， 2010)。張慶蓮等(2010)通過(guò)野外實(shí)測(cè)發(fā)現(xiàn)，斷裂控制的構(gòu)造裂縫的密度與距斷裂的距離呈 指數(shù)遞減關(guān)系，越靠近破裂面，裂縫越發(fā)育。Billi A(2003)等把這種斷裂結(jié)構(gòu)劃分為斷層 核和破碎帶兩部分，斷層核是一狹窄帶，由主滑動(dòng)面和斷層巖構(gòu)成；破碎帶位于斷層核周 緣，主要由裂縫發(fā)育帶和小斷層組成（Agosta F et alJOOehTveranger J(2005)、 Braathen A(2009)等把這種受斷裂作用影響的巖石三維空間定義為斷層包絡(luò)體，包括破碎 帶和誘導(dǎo)裂縫帶2個(gè)部分，其中破碎帶以發(fā)育斷層巖和伴生裂縫為主要特征，誘導(dǎo)裂縫帶內(nèi) 主要發(fā)育各種類型的誘導(dǎo)裂縫(付曉飛等，2005)?？傊?，斷層不是一個(gè)簡(jiǎn)單的"面"，而是一 個(gè)"帶"（付廣等，2008)，具有復(fù)雜的內(nèi)部結(jié)構(gòu)(Caine，1996;Gudmundsson A，2001;付曉飛 等，2005)。吳智平等(2010)依據(jù)變形程度，將斷裂帶劃分為滑動(dòng)破碎帶和誘導(dǎo)裂縫帶兩類 結(jié)構(gòu)單元。
[0003] 壓性斷裂是在擠壓應(yīng)力場(chǎng)作用下形成的逆斷層。其受力作用強(qiáng)，沿滑動(dòng)面兩側(cè)巖 石破碎嚴(yán)重，斷裂帶結(jié)構(gòu)完整，對(duì)油氣成藏控制明顯，是我國(guó)西部盆地油氣運(yùn)聚的關(guān)鍵因 素。通常，壓性斷裂帶的滑動(dòng)破碎帶內(nèi)巖石破碎嚴(yán)重，呈巖粉狀，多發(fā)育斷層泥，部分糜棱巖 化，流體通過(guò)能力差，對(duì)油氣具有較強(qiáng)的封堵性(吳孔友等，2012)。而誘導(dǎo)裂縫帶發(fā)育大量 裂縫，流體通過(guò)能力強(qiáng)，為油氣運(yùn)移提供了通道(付曉飛等，2005)?？傊?，斷裂帶對(duì)油氣具有 輸導(dǎo)和封堵雙重能力。即斷裂帶內(nèi)部結(jié)構(gòu)的時(shí)空差異將導(dǎo)致其在油氣運(yùn)聚過(guò)程中扮演不同 的角色。因此，預(yù)測(cè)含油氣區(qū)地下深部發(fā)育的斷裂帶結(jié)構(gòu)特征，對(duì)油氣藏勘探具有極為重要 的意義。
[0004] 對(duì)于覆蓋區(qū)壓性斷裂而言，由于不能直接觀察，因此，如何準(zhǔn)確的識(shí)別斷裂結(jié)構(gòu)一 直是石油地質(zhì)勘探中的難題之一。劉偉等(2013)曾提出利用電阻率差比、孔隙度比值、曲線 變化率及井徑增大率四個(gè)參數(shù)，并結(jié)合主成分分析法來(lái)定量的識(shí)別斷裂帶內(nèi)部結(jié)構(gòu)。然而， 該方法由于參數(shù)選擇少、未考慮曲線變化幅度等原因?qū)е略趹?yīng)用過(guò)程中不能有效的對(duì)斷裂 帶內(nèi)部結(jié)構(gòu)進(jìn)行識(shí)別。趙永剛等（2013)提出三孔隙度比值、等效彈性模量差比、次生孔隙 度、雙感應(yīng)幅度差指標(biāo)、龜裂系數(shù)、井徑相對(duì)異常、電阻率侵入校正及膠結(jié)指數(shù)八種參數(shù)為 判斷地層中裂縫發(fā)育的有效參數(shù)。由斷裂帶內(nèi)部誘導(dǎo)裂縫帶和滑動(dòng)破碎帶的特征可知，裂 縫是否發(fā)育是這兩部分結(jié)構(gòu)重要的區(qū)別之一，因此，挑選裂縫敏感性指標(biāo)，構(gòu)建斷裂結(jié)構(gòu)綜 合判別參數(shù)，可以從測(cè)井曲線上有效的判斷斷裂帶內(nèi)部結(jié)構(gòu)。

【發(fā)明內(nèi)容】

[0005] 本發(fā)明的目的在于提供一種壓性斷裂帶結(jié)構(gòu)測(cè)井識(shí)別方法，通過(guò)多參數(shù)測(cè)井曲線 擬合，構(gòu)建斷裂帶結(jié)構(gòu)綜合判別參數(shù)，結(jié)合斷裂帶內(nèi)部結(jié)構(gòu)的特點(diǎn)，定量的劃分覆蓋區(qū)壓性 斷裂內(nèi)部結(jié)構(gòu)，為準(zhǔn)確評(píng)價(jià)斷裂的啟閉性，提高油氣勘探與預(yù)測(cè)的有效性提供技術(shù)支撐。
[0006] 為解決上述技術(shù)問(wèn)題，本發(fā)明采用以下技術(shù)方案予以實(shí)現(xiàn)：
[0007] -種壓性斷裂帶結(jié)構(gòu)測(cè)井識(shí)別方法，包括如下步驟：
[0008] (1)判斷工區(qū)應(yīng)力環(huán)境是否為壓性應(yīng)力環(huán)境，若為壓性應(yīng)力環(huán)境，進(jìn)入步驟(2);
[0009] (2)遴選并計(jì)算裂縫敏感性參數(shù):①深淺電阻率差比RTC、②三孔隙度測(cè)井、③曲線 變化率、④井徑增大率CALd
[0010] ①深淺電阻率差比RTC:
[0011] 其中，RD為深側(cè)向測(cè)井值或深電阻率測(cè)井值，RS為淺側(cè)向測(cè)井值或淺電阻率測(cè)井 值;得到曲線F(Rtc)
[0012] ②三孔隙度測(cè)井
[0013] 三孔隙度測(cè)井曲線包括聲波時(shí)差(AC)、密度(DEN)及補(bǔ)償中子(CNL);其中，利用聲 波時(shí)差計(jì)算的孔隙度為原生孔隙度，利用密度和補(bǔ)償中子計(jì)算的孔隙度為至少包括裂縫、 溶蝕孔的次生孔隙度在內(nèi)的總孔隙度；
[0014] 聲波時(shí)差計(jì)算原生孔隙度公式為：
[0015] 其中，△ t為從聲波時(shí)差曲線讀出的地層聲波時(shí)差值，△ tf為孔隙中流體的聲波時(shí) 差值，A tma為巖石骨架的聲波時(shí)差值；
[0016] 密度測(cè)井計(jì)算總孔隙度公式為：
[0017] 其中，pb為密度測(cè)井曲線讀出的測(cè)量值，pf為孔隙中流體的密度值，pma為巖石骨架 密度值；
[0018] 補(bǔ)償中子測(cè)井計(jì)算總孔隙度公式為：
[0019] 其中，變&為補(bǔ)償中子曲線上讀出的測(cè)量值，攀激^為巖石骨架中子值，熟貧為孔隙 中流體的中子值；

[0021]由裂縫等因素引起的次生孔隙度變化值 [0022] ③曲線變化率
[0023]選取對(duì)裂縫敏感的曲線聲波時(shí)差(AC)、補(bǔ)償中子(CNL)及密度(DEN)進(jìn)行曲線變化 率計(jì)算，計(jì)算公式為：
[0024] Δ Xi= ( | Xi-i-Xi | + | Xi+i-Xi | )/Xi
[0025] 其中，XiS當(dāng)前深度點(diǎn)的測(cè)井曲線值，Xh和X1+1為當(dāng)前深度點(diǎn)鄰近的兩點(diǎn)的測(cè)井 曲線值；
[0026] ④井徑增大率CALd
[0028] CAL為鉆頭直徑，CALj為實(shí)測(cè)井徑；
[0029] (3)針對(duì)各個(gè)參數(shù)進(jìn)行歸一化計(jì)算。
[0030] 曲線(AC、CNL)變化率、三孔隙度、井徑增大率采用如下公式歸一化：
[0031] Y=(X-Xmin)/(Xmax-Xmin)
[0032] 其中，X表示曲線(AC、CNL)變化率、三孔隙度或井徑增大率實(shí)際測(cè)量值，Xmax表示該 種測(cè)井曲線的曲線(AC、CNL)變化率、三孔隙度或井徑增大率最大測(cè)量值，Xmin表示該種測(cè)井 曲線的曲線(AC、CNL)變化率、三孔隙度或井徑增大率最小測(cè)量值;得到曲線F( AAC)、F( Δ CNL)、F( ACALd)、F( Δ ΦΡ);
[0033] 密度(DEN)變化率采用如下公式歸一化：
[0034] Y=l-(X-Xmin)/(Xmax-Xmin)；
[0035] 其中，X表示密度測(cè)井測(cè)量值，Xmax表示密度測(cè)井的最大測(cè)量值，Xmin表示密度測(cè)井 的最小測(cè)量值;得到曲線F( ADEN);
[0036] (4)構(gòu)建斷裂帶內(nèi)部結(jié)構(gòu)綜合判別參數(shù)FZI
[0037] FZI = [F( AAC)+F( ACNL)+F( ADEN)+F(Rtc)+F( ACALd)+F( Δ ΦΡ)]/6
[0038] 其中，F(xiàn)為相對(duì)各種曲線的函數(shù)，滿足以下條件：
[0040] 式中，X代表 AAC、ACNL、ADEN、RTC、ACALd或 Δ (DP，a、b為異常臨界值，且有a〈b。
[0041] 本發(fā)明的優(yōu)點(diǎn)和積極效果是:本發(fā)明可以有效的利用測(cè)井資料對(duì)斷裂帶內(nèi)部結(jié)構(gòu) 進(jìn)行識(shí)別，并且選取的曲線類型為常規(guī)的標(biāo)準(zhǔn)測(cè)井，具有很強(qiáng)的可操作性。與現(xiàn)有技術(shù)(定 性的利用單條曲線變化趨勢(shì)進(jìn)行斷裂帶結(jié)構(gòu)劃分)相比，將斷裂帶內(nèi)部結(jié)構(gòu)測(cè)井曲線識(shí)別 的過(guò)程定量化，更具有說(shuō)服力，因而本發(fā)明對(duì)后續(xù)開展斷裂封閉性及斷裂控藏作用研究具 有重大的意義。本發(fā)明利用巖心觀察及成像測(cè)井資料對(duì)該識(shí)別方案進(jìn)行了檢驗(yàn)，發(fā)現(xiàn)該劃 分與巖心及成像測(cè)井顯示相符，驗(yàn)證了該方法的準(zhǔn)確性。
[0042] 結(jié)合附圖閱讀本發(fā)明實(shí)施方式的詳細(xì)描述后，本發(fā)明的其他特點(diǎn)和優(yōu)點(diǎn)將變得更 加清楚。
【附圖說(shuō)明】
[0043]圖1為本發(fā)明具體實(shí)施例的實(shí)施路線圖。
[0044]圖2為排66井?dāng)嗔褍?nèi)部結(jié)構(gòu)測(cè)井曲線定量識(shí)別綜合圖。
【具體實(shí)施方式】
[0045]下面結(jié)合附圖對(duì)本發(fā)明的【具體實(shí)施方式】進(jìn)行詳細(xì)地描述。
[0046]如圖1所示，對(duì)構(gòu)造應(yīng)力場(chǎng)進(jìn)行分析，判斷工區(qū)應(yīng)力環(huán)境是否為壓性應(yīng)力環(huán)境，若 為壓性應(yīng)力環(huán)境，繼續(xù)如下步驟。
[0047] 對(duì)目標(biāo)鉆井進(jìn)行選取，對(duì)測(cè)井曲線進(jìn)行篩選，確定曲線的異常值a、b。
[0048] 遴選并計(jì)算裂縫敏感性參數(shù):①深淺電阻率差比Rtc、②三孔隙度測(cè)井、③曲線變化 率、④井徑增大率CALd
[0050]其中，Rd為深側(cè)向測(cè)井值或深電阻率測(cè)井值，Rs為淺側(cè)向測(cè)井值或淺電阻率測(cè)井 值;得到曲線F( Rtc);
[0051 ]②三孔隙度測(cè)井
[0052] 三孔隙度測(cè)井曲線包括聲波時(shí)差(AC)、密度(DEN)及補(bǔ)償中子(CNL);其中，利用聲 波時(shí)差計(jì)算的孔隙度為原生孔隙度，利用密度和補(bǔ)償中子計(jì)算的孔隙度為至少包括裂縫、 溶蝕孔的次生孔隙度在內(nèi)的總孔隙度；
[0053] 聲波時(shí)差計(jì)算原生孔隙度公式為：
[0054] 其中，△ t為從聲波時(shí)差曲線讀出的地層聲波時(shí)差值，△ tf為孔隙中流體的聲波時(shí) 差值，A tma為巖石骨架的聲波時(shí)差值；
[0055] 密度測(cè)井計(jì)算總孔隙度公式為：
[0056] 其中，pb為密度測(cè)井曲線讀出的測(cè)量值，pf為孔隙中流體的密度值，pma為巖石骨架 密度值；
[0057] 補(bǔ)償中子測(cè)井計(jì)算總孔隙度公式為：
[0058] 其中，Λ為補(bǔ)償中子曲線上讀出的測(cè)量值，購(gòu)_為巖石骨架中子值，變尊為孔隙 中流體的中子值；

[0060] 由裂縫等因素引起的次生孔隙度變化值
[0061] ③曲線變化率
[0062] 選取對(duì)裂縫敏感的曲線聲波時(shí)差(AC)、補(bǔ)償中子(CNL)及密度(DEN)進(jìn)行曲線變化 率計(jì)算，計(jì)算公式為：
[0063] Δ Xi= ( | Xi-i-Xi | + | Xi+i-Xi | )/Xi
[0064] 其中，Xi為當(dāng)前深度點(diǎn)的測(cè)井曲線值，Xi-i和Xi+i為當(dāng)前深度點(diǎn)鄰近的兩點(diǎn)的測(cè)井 曲線值；
[0065] ④井徑增大率CALd
[0067] CAL為鉆頭直徑，CALj為實(shí)測(cè)井徑；
[0068] (3)針對(duì)各個(gè)參數(shù)進(jìn)行歸一化計(jì)算。
[0069]曲線(AC、CNL)變化率、三孔隙度、井徑增大率采用如下公式歸一化：
[0070] Y=(X-Xmin)/(Xmax-Xmin)
[0071] 其中，X表示曲線(AC、CNL)變化率、三孔隙度或井徑增大率實(shí)際測(cè)量值，Xmax表示該 種測(cè)井曲線的曲線(AC、CNL)變化率、三孔隙度或井徑增大率最大測(cè)量值，X min表示該種測(cè)井 曲線的曲線(AC、CNL)變化率、三孔隙度或井徑增大率最小測(cè)量值;得到曲線F( AAC)、F( Δ CNL)、F( ACALd)、F( Δ ΦΡ);
[0072] 密度(DEN)變化率采用如下公式歸一化：
[0073] Y=l-(X-Xmin)/(Xmax-Xmin)；
[0074] 其中，X表示密度測(cè)井測(cè)量值，Xmax表示密度測(cè)井的最大測(cè)量值，Xmin表示密度測(cè)井 的最小測(cè)量值;得到曲線F( ADEN);
[0075] (4)構(gòu)建斷裂帶內(nèi)部結(jié)構(gòu)綜合判別參數(shù)FZI
[0076] FZI = [F( AAC)+F( ACNL)+F( ADEN)+F(Rtc)+F( ACALd)+F( Δ ΦΡ)]/6
[0077] 其中，F(xiàn)為相對(duì)各種曲線的函數(shù)，滿足以下條件：
[0079] 式中，X代表 AAC、ACNL、ADEN、RTC、ACALd或 Δ (DP，a、b為異常臨界值，且有a〈b。
[0080] 本方法在準(zhǔn)噶爾盆地西北緣實(shí)驗(yàn)時(shí)取得了良好的實(shí)驗(yàn)效果。準(zhǔn)噶爾盆地是我國(guó)西 部重要的含油氣盆地，經(jīng)歷海西期至喜馬拉雅期構(gòu)造運(yùn)動(dòng)，壓性斷裂發(fā)育，斷裂活動(dòng)時(shí)間 長(zhǎng)、平面延伸距離遠(yuǎn)、垂向斷距大，且部分井鉆穿斷裂帶。利用該方法對(duì)鉆穿紅車斷裂帶的 排66井進(jìn)行了斷裂帶內(nèi)部結(jié)構(gòu)識(shí)別，效果良好(圖2)。紅車斷裂帶在該井可以分為上盤誘導(dǎo) 裂縫帶、滑動(dòng)破碎帶和下盤誘導(dǎo)裂縫帶三部分。相對(duì)于兩盤誘導(dǎo)裂縫帶來(lái)說(shuō)，滑動(dòng)破碎帶Rtc 值為負(fù)或零，顯示裂縫不發(fā)育;CALd值明顯偏小且較穩(wěn)定，暗示井徑垮塌程度弱于誘導(dǎo)裂縫 帶；ADEN較為穩(wěn)定，但兩側(cè)誘導(dǎo)裂縫帶則出現(xiàn)明顯的變化；AAC、ACNL較小且變化較小，但 誘導(dǎo)裂縫帶則出現(xiàn)了較為明顯的跳動(dòng)；Φρ表現(xiàn)穩(wěn)定，與巖性整體為泥巖有關(guān);FZI則呈現(xiàn)明 顯低值，表明在滑動(dòng)破碎帶裂縫整體不發(fā)育，而在誘導(dǎo)裂縫帶中裂縫發(fā)育明顯。
[0081] 最后應(yīng)說(shuō)明的是：以上實(shí)施例僅用以說(shuō)明本發(fā)明的技術(shù)方案，而非對(duì)其限制;盡管 參照前述實(shí)施例對(duì)本發(fā)明進(jìn)行了詳細(xì)的說(shuō)明，本領(lǐng)域的普通技術(shù)人員應(yīng)當(dāng)理解:其依然可 以對(duì)前述各實(shí)施例所記載的技術(shù)方案進(jìn)行修改，或者對(duì)其中部分技術(shù)特征進(jìn)行等同替換； 而這些修改或者替換，并不使相應(yīng)技術(shù)方案的本質(zhì)脫離本發(fā)明各實(shí)施例技術(shù)方案的精神和 范圍。
【主權(quán)項(xiàng)】
1. 一種壓性斷裂帶結(jié)構(gòu)測(cè)井識(shí)別方法，其特征在于，包括如下步驟： (1) 判斷工區(qū)應(yīng)力環(huán)境是否為壓性應(yīng)力環(huán)境，若為壓性應(yīng)力環(huán)境，進(jìn)入步驟(2); (2) 遽選并計(jì)算裂縫敏感性參數(shù):①深淺電阻率差比Rtc、②Ξ孔隙度測(cè)井、③曲線變化 率、④井徑增大率CALd; ① 深淺電阻率差比化C:其中，Rd為深側(cè)向測(cè)井值或深電阻率測(cè)井值，Rs為淺側(cè)向測(cè)井值或淺電阻率測(cè)井值;得 到曲線F(Rtc); ② Ξ孔隙度測(cè)井 Ξ孔隙度測(cè)井曲線包括聲波時(shí)差(AC)、密度(DEN)及補(bǔ)償中子((m);其中，利用聲波時(shí) 差計(jì)算的孔隙度為原生孔隙度，利用密度和補(bǔ)償中子計(jì)算的孔隙度為至少包括裂縫、溶蝕 孔的次生孔隙度在內(nèi)的總孔隙度； 聲波時(shí)差計(jì)算原生孔隙度公式為：其中，A t為從聲波時(shí)差曲線讀出的地層聲波時(shí)差值，Δ tf為孔隙中流體的聲波時(shí)差值， A tma為巖石骨架的聲波時(shí)差值； 密度測(cè)井計(jì)算總孔隙度公式為：其中，Pb為密度測(cè)井曲線讀出的測(cè)量值，Pf為孔隙中流體的密度值，Pma為巖石骨架密度 值； 補(bǔ)償中子測(cè)井計(jì)算總孔隙度公式為：其中，鱗為補(bǔ)償中子曲線上讀出的巧慢值，輸胃為巖石骨架中子值，資游為孔隙中流 體的中子值； 總孔隙度值'7 由裂縫等因素引起的次生孔隙度變化值③ 曲線變化率 選取對(duì)裂縫敏感的曲線聲波時(shí)差(AC)、補(bǔ)償中子(CNL)及密度(DEN)進(jìn)行曲線變化率計(jì) 算，計(jì)算公式為： AXi = (|Xi-i-Xi| + |Xi+i-Xi|)/Xi； 其中，Xi為當(dāng)前深度點(diǎn)的測(cè)井曲線值，Xi-1和Xi+1為當(dāng)前深度點(diǎn)鄰近的兩點(diǎn)的測(cè)井曲線 值； ④井徑增大率CALd 計(jì)算公式為：Z： CAL為鉆頭直徑，CAb為實(shí)測(cè)井徑； (3) 針對(duì)各個(gè)參數(shù)進(jìn)行歸一化計(jì)算： 曲線(AC、CNL)變化率、Ξ孔隙度、井徑增大率采用如下公式歸一化： Y=(X-Xmin)/(Xmax-Xmin) 其中，X表示曲線(AC、CNL)變化率、Ξ孔隙度或井徑增大率實(shí)際測(cè)量值，Xmax表示該種測(cè) 井曲線的曲線(AC、C化)變化率、Ξ孔隙度或井徑增大率最大測(cè)量值，Xmin表示該種測(cè)井曲線 的曲線(AC、(m)變化率、Ξ孔隙度或井徑增大率最小測(cè)量值;得到曲線F( AAC)、F( ACNL)、 F(ACALd)、F(A 巫p); 密度化ΕΝ)變化率采用如下公式歸一化： Y=l-(X-Xmin)/(Xmax-Xmin)； 其中，X表示密度測(cè)井測(cè)量值，Xmax表示密度測(cè)井的最大測(cè)量值，Xmin表示密度測(cè)井的最 小測(cè)量值;得到曲線F( ADEN); (4) 構(gòu)建斷裂帶內(nèi)部結(jié)構(gòu)綜合判別參數(shù)FZI FZI = [F( AAC)+F( ACNL)+F( ADEN)+F(化c)+F( ACALd)+F( Δ Φρ)]/6 其中，F(xiàn)為相對(duì)各種曲線的函數(shù)，滿足W下條件：式中，X代表AAC、ACNL、ADEN、化C、ACALd或Δ Φρ，3、6為異常臨界值，且有a<b。
【文檔編號(hào)】E21B49/00GK106089191SQ201610411508
【公開日】2016年11月9日
【申請(qǐng)日】2016年6月12日 公開號(hào)201610411508.X, CN 106089191 A, CN 106089191A, CN 201610411508, CN-A-106089191, CN106089191 A, CN106089191A, CN201610411508, CN201610411508.X
【發(fā)明人】吳孔友, 劉寅, 吳曉菲, 劉波, 裴仰文
【申請(qǐng)人】中國(guó)石油大學(xué)(華東), 北京大學(xué)附屬中學(xué)