亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        基于條件生成對抗網(wǎng)絡(luò)的成像測井圖像裂縫計算機識別

        2020-03-30 02:00:14魏伯陽潘保芝殷秋麗田超國郭曉雅
        石油物探 2020年2期
        關(guān)鍵詞:正弦曲線標(biāo)簽卷積

        魏伯陽,潘保芝,殷秋麗,田超國,郭曉雅

        (1.吉林大學(xué)地球探測科學(xué)與技術(shù)學(xué)院,吉林長春130026;2.冀東油田公司勘探開發(fā)研究院,河北唐山063000;3.鄭州中核巖土工程有限公司,河南鄭州450008)

        裂縫由構(gòu)造變形作用或物理成巖作用形成,是在巖石中天然存在的宏觀面狀不連續(xù)[1]。裂縫既是油氣的儲集空間,又是油氣的流動通道。裂縫與油氣的儲量和產(chǎn)能息息相關(guān),裂縫識別是油氣儲量評價和產(chǎn)能預(yù)測的關(guān)鍵。由于地下情況復(fù)雜,盡管目前對于地層微電阻率掃描成像(formation microscanner image,FMI)圖像中的裂縫已有成熟的識別技術(shù),但大多是基于人機交互的方法,人力資源消耗大(耗功耗時)且易受主觀因素影響,因而對裂縫判斷不夠精確。國內(nèi)外多位專家和學(xué)者進(jìn)行了裂縫自動識別技術(shù)研究。秦巍等[2]、李雪英等[3]引入圖像處理手段,利用裂縫灰度特性識別裂縫,實現(xiàn)了對地層中簡單裂縫信息的自動識別。張程恩[4]使用蟻群聚類算法提取成像圖裂縫形態(tài)。唐佳偉[5]使用二維Otsu算法,消除成像圖的噪聲干擾,并且對提取的裂縫圖像進(jìn)行細(xì)化。張群會等[6]提出了基于元胞自動機模型的裂縫自動識別算法,確定FMI圖像中裂縫的灰度范圍,利用免疫規(guī)則識別裂縫。張曉峰等[7]對FMI圖像進(jìn)行垂直和橫向的小波變換,剔除噪聲的同時重構(gòu)小波變換圖,最終得到裂縫圖像。潘保芝等[8]提出綜合常規(guī)測井?dāng)?shù)據(jù)裂縫信息與FMI圖像中裂縫響應(yīng)識別裂縫區(qū)域。低角度斜交裂縫的形態(tài)在FMI圖像中表現(xiàn)為一個單周期的正弦曲線,但因為地下結(jié)構(gòu)復(fù)雜(多條裂縫交叉、井壁坍塌等)及儀器自身原因(提升速度變化、壞極板等),造成裂縫面表現(xiàn)為一個不是完全成正弦曲線的形態(tài),因此給計算機自動識別FMI圖像中的裂縫帶來嚴(yán)重干擾。利用成像測井資料識別裂縫的方法雖然識別效果較好,但都存在共同的缺點,即裂縫識別費時費力,且識別的裂縫中包括無效縫[9]。

        2006年,HINTON等[10]提出了有效提取數(shù)據(jù)(包括圖像、音頻、文本等)的抽象特征的深度學(xué)習(xí)方法,并對數(shù)據(jù)進(jìn)行分類與回歸。卷積神經(jīng)網(wǎng)絡(luò)是深度學(xué)習(xí)中的一種重要方法。卷積神經(jīng)網(wǎng)絡(luò)對圖像與圖像間的非線性關(guān)系有良好的表征能力,可以對圖像抽象出人工難以理解的特征進(jìn)而對圖像進(jìn)行識別分類。目前該方法已經(jīng)在基于常規(guī)測井曲線的測井?dāng)?shù)據(jù)分類領(lǐng)域取得很高的成功率,如安鵬等[11]應(yīng)用深度學(xué)習(xí)方法劃分9種巖性的準(zhǔn)確率達(dá)90.02%,吳正陽等[12]應(yīng)用卷積神經(jīng)網(wǎng)絡(luò)劃分裂縫性儲層類別與人工解釋結(jié)果也有相當(dāng)高的符合率。由于FMI圖像受井眼、地層等條件的影響,包含的噪聲較多,因此,普通的卷積神經(jīng)網(wǎng)絡(luò)模型難以識別FMI圖像中的裂縫信息。

        2014年,GOODFELLOW等[13]在卷積神經(jīng)網(wǎng)絡(luò)的基礎(chǔ)上提出生成對抗網(wǎng)絡(luò)(generative adversarial nets,GAN),通過對實際圖片進(jìn)行無監(jiān)督學(xué)習(xí),然后利用噪聲生成與實際圖片特征相似的圖片。MIRZA等[14]在GAN的基礎(chǔ)上提出條件生成對抗網(wǎng)絡(luò)(conditional generative adversarial nets,CGAN)。CGAN將條件變量加入到GAN模型的輸入中,將無監(jiān)督的GAN變?yōu)橛斜O(jiān)督模型[14]。本文借助CGAN思想識別計算機模擬的裂縫,介紹了CGAN方法的抗干擾能力,并采用CGAN進(jìn)行了實際測井圖像的裂縫識別。

        1 方法技術(shù)

        CGAN是在GAN的基礎(chǔ)上添加條件標(biāo)簽圖像,利用該標(biāo)簽圖像約束生成器從實際圖像中提取信息,將GAN無監(jiān)督的訓(xùn)練過程變?yōu)橛斜O(jiān)督的訓(xùn)練過程。

        1.1 生成對抗網(wǎng)絡(luò)(GAN)

        GAN由生成器(Generator)和判別器(Discriminator)兩部分組成。GAN損失函數(shù)V(D,G)由生成器最小化損失函數(shù)(minG)和判別器最大化損失函數(shù)(maxD)構(gòu)成:

        minGmaxDV(D,G)=Ex~Pdata(x)[log(D(x))]+
        Ez~Pz(z)[log(1-D(G(z)))]

        (1)

        式中:x表示實際圖像;z表示噪聲圖像;D(x)表示判別器對實際圖像的判別值;D(G(z))表示判別器對生成圖像G(z)的判別值;Pdata(x)表示實際圖像的分布;Pz(z)表示噪聲圖像分布;Ex~Pdata(x)表示x服從實際圖像分布取樣的期望;Ez~Pz(z)表示z服從噪聲圖像分布取樣的期望。

        生成器是以z作為輸入,對x通過卷積提取圖像特征信息,將z生成為服從實際圖像分布的圖像G(z)。生成器的損失函數(shù)V(G)越小,G(z)與x越相似。判別器將G(z)與x同時作為輸入,判斷G(z)與x是否一致。判別器的損失函數(shù)V(D)在[0,1]之間,當(dāng)V(D)為1時,生成的圖像G(z)與實際圖像x一致;當(dāng)V(D)為0時,G(z)與x完全不一致。V(D)越大,G(z)與x越相似。判別器使G(z)盡可能相似于x,最終得到D(x)。生成器和判別器不斷博弈,最終使生成器生成的G(z)無限接近于x。

        1.2 條件生成對抗網(wǎng)絡(luò)

        條件生成對抗網(wǎng)絡(luò)CGAN的損失函數(shù)為:

        minGmaxDV(D,G)=Ex~Pdata(x)[log(D(x|y))]+
        Ez~Pz(z)[log(1-D(G(z|y)))]

        (2)

        式中:y表示條件即標(biāo)簽圖像;D(x|y)表示判別器對實際圖像x與標(biāo)簽圖像y拼接后的判別值;G(z|y)表示噪聲圖像z與標(biāo)簽圖像y拼接后生成器生成的圖像;D(G(z|y))表示判別器對于G(z|y)與y拼接后的判別值。從損失函數(shù)((2)式)可以看出,y不僅參與了判別器的輸入,同時也參與了生成器的輸入。

        CGAN也由生成器和判別器兩部分組成(圖1)。生成器以隨機噪聲圖像z以及裂縫標(biāo)簽圖像y為輸入,生成服從實際圖像數(shù)據(jù)分布的圖像G(z|y),判別器將G(z|y)與x、y作為輸入,判斷G(z|y)分別與x和y拼接后的判別值是否一致,最終提取出只包含標(biāo)簽特征信息的圖像。

        圖2為CGAN單次迭代(單張圖片訓(xùn)練一次)識別裂縫流程。生成過程首先輸入y和z,將二者數(shù)據(jù)拼接之后變成像素為256*256、層數(shù)為3的圖像(用256*256@3表示)。輸入到生成器,生成器中包含8個卷積層和8個反卷積層,輸出通過激發(fā)函數(shù)tanh生成G(z|y)。判別過程將x和y與G(z|y)分別拼接,改變圖像尺寸為256*256@3后輸入到判別器。判別器中包含4個卷積層,將256*256@3的輸入圖像卷積為32*32@1的特征圖像。通過激發(fā)函數(shù)sigmoid得到D(x|y)。其中輸入x和y到判別器,得到判別器對實際圖像的判別結(jié)果,輸入x與G(z|y)到判別器,得到判別器對生成圖像的判別結(jié)果。

        圖1 CGAN識別裂縫模型

        圖2 CGAN單次迭代(單張圖片訓(xùn)練一次)識別裂縫流程

        2 模擬圖像裂縫識別

        天然裂縫按角度分為水平縫(傾角<5°)、低角度縫(5°~30°)、中角度縫(30°~75°)、高角度縫(傾角>75°)和豎直縫(傾角=90°)[15]。其中低角度縫、中角度縫及高角度縫在FMI圖像中響應(yīng)為黑色正弦形態(tài),水平縫表現(xiàn)為水平方向黑色直線,豎直縫響應(yīng)為垂向黑色貫穿直線(本文暫不考慮豎直縫)。

        2.1 CGAN訓(xùn)練

        使用計算機模擬裂縫圖像,模擬不同方位、傾角和寬度的單條裂縫圖片共216張,其中水平縫、低角度縫、中角度縫、高角度縫的圖片分別為36,36,72,72張,如圖3左圖所示。圖3右圖為相應(yīng)的標(biāo)簽圖像。將所有訓(xùn)練圖像訓(xùn)練一次的過程稱為1個循環(huán)(epoch),本次訓(xùn)練過程經(jīng)過100個epoch,每個epoch迭代216次(訓(xùn)練一張圖片視為迭代一次),共迭代21600次,耗時30min。圖4為CGAN計算機模擬裂縫圖像訓(xùn)練結(jié)果,可以看出,訓(xùn)練得到的結(jié)果與標(biāo)簽一致。

        2.2 模擬裂縫圖像的CGAN識別

        利用訓(xùn)練好的CGAN識別未參與訓(xùn)練的圖像,并觀察CGAN模型的裂縫識別效果及對圖像噪聲的抗干擾能力。構(gòu)造的4張裂縫圖像分別為:沿正弦曲線添加了隨機噪聲、正弦曲線中間斷開、對整個圖像添加了隨機噪聲和2條正弦曲線疊加圖像,如圖5 所示。測試發(fā)現(xiàn):CGAN對有噪聲的不規(guī)則正弦曲線(圖5a)和包含脈沖噪聲的圖像(圖5c)有較強的抗干擾能力,識別結(jié)果和裂縫標(biāo)簽形態(tài)一致;對疊加裂縫也有良好的識別結(jié)果(圖5d);在識別中間斷開的正弦曲線形態(tài)時,得到中間斷開的正弦曲線(圖5b)。分析其原因,認(rèn)為:CGAN強烈依賴訓(xùn)練樣本和標(biāo)簽,而上述訓(xùn)練樣本只包含單條連續(xù)正弦曲線樣本,故不能將中間斷開的正弦曲線識別為連續(xù)正弦曲線,需要在訓(xùn)練過程中加入斷開的正弦曲線樣本進(jìn)行訓(xùn)練。最后設(shè)置100張各種干擾(沿曲線添加隨機噪聲、曲線斷開、多條曲線交叉等)下計算機模擬裂縫圖像,結(jié)果顯示,CGAN識別這些圖像中裂縫的準(zhǔn)確率達(dá)93.4%。

        圖3 計算機模擬裂縫圖像(左)及標(biāo)簽(右)a 水平縫; b 低角度縫; c 中角度縫; d 高角度縫

        圖4 CGAN計算機模擬裂縫圖像訓(xùn)練結(jié)果a 水平縫; b 低角度縫; c 中角度縫; d 高角度縫

        圖5 CGAN計算機模擬裂縫識別結(jié)果a 正弦曲線添加隨機噪聲; b 正弦曲線中間斷開; c 圖像添加隨機噪聲; d 多條曲線疊加

        3 FMI圖像裂縫識別

        3.1 CGAN訓(xùn)練

        將CGAN應(yīng)用于實際FMI圖像識別裂縫訓(xùn)練中,并將實際FMI圖像和相應(yīng)的裂縫標(biāo)簽作為輸入直接進(jìn)行訓(xùn)練。首先利用人機交互方法提取裂縫標(biāo)簽,共得到162個訓(xùn)練樣本,其中,單條裂縫圖像124張(水平縫49張,低角度縫62張,中角度縫13張),含兩條及兩條以上裂縫的圖片38張。圖6a、圖6b和圖6c分別代表其水平縫、低角度縫、中角度縫的FMI圖像及標(biāo)簽圖像。

        訓(xùn)練過程經(jīng)過500個epoch,每個epoch迭代162次,共迭代81000次,耗時2.5h。圖7展示了其中4張具有代表性的裂縫圖像及其訓(xùn)練結(jié)果,可以看出,CGAN對水平縫(圖7a)和低角度縫(圖7b)有較好的訓(xùn)練結(jié)果。訓(xùn)練過程中發(fā)現(xiàn)有3處異常,如圖8所示。分析圖8a和圖8c發(fā)現(xiàn),由于裂縫形狀的不規(guī)則或者周圍有微小裂縫與主裂縫連通造成生成的裂縫圖像與標(biāo)簽設(shè)定的裂縫有差異;圖8b中FMI圖像上部分可能存在裂縫但標(biāo)簽中沒有標(biāo)注出該條裂縫。

        圖6 FMI裂縫圖像(左)及標(biāo)簽(右)a 水平縫; b 低角度縫; c 中角度縫

        圖7 CGAN FMI圖像裂縫訓(xùn)練結(jié)果a 水平縫; b 低角度縫; c 中角度縫; d 多條曲線疊加

        圖8 CGAN訓(xùn)練過程判別器異常點a 異常1; b 異常2; c 異常3

        可見,設(shè)置標(biāo)簽的準(zhǔn)確性會影響CGAN的訓(xùn)練結(jié)果和識別效果。為避免此類異常出現(xiàn),設(shè)置裂縫標(biāo)簽時應(yīng)盡可能準(zhǔn)確地包含圖像中所有裂縫信息。

        3.2 FMI圖像裂縫CGAN識別

        用訓(xùn)練好的CGAN識別10張包含裂縫的圖像,識別準(zhǔn)確率達(dá)90%,識別速度為0.19s/張。圖9為訓(xùn)練模型識別結(jié)果,可見,當(dāng)水平縫與高角度縫交叉時,CGAN對水平縫也有良好的識別效果(圖9b);低角度縫周圍的破碎使裂縫變?yōu)椴灰?guī)則的正弦形態(tài),CGAN也可以避開干擾準(zhǔn)確識別裂縫(圖9c);當(dāng)中角度縫不是正弦形態(tài)的裂縫時,CGAN識別結(jié)果與標(biāo)簽之間出現(xiàn)一定的差異(圖9d)。分析原因,首先,受客觀條件限制,中角度縫訓(xùn)練樣本相對較少,其次,中角度縫破碎后裂縫形態(tài)不是正弦形態(tài)且易與其它裂縫形成貫通,因此對中角度縫的識別會有較大誤差。

        圖9 FMI圖像CGAN網(wǎng)絡(luò)識別裂縫結(jié)果a 水平縫; b 水平縫受高角度縫交叉干擾; c 不規(guī)則水平縫; d 不規(guī)則中角度縫

        目前采用的FMI圖像資料中截取到的中角度縫圖像較少,因而參與訓(xùn)練的中角度縫圖像較少,使得對該資料中少量的中角度縫尚不能有良好的識別效果,需要進(jìn)一步深入研究。由計算機模擬裂縫識別效果可知,在樣本數(shù)據(jù)充足的情況下,采用CGAN對此類裂縫有良好的識別效果。

        3.3 裂縫識別效果對比

        分別采用CGAN和蟻群算法對FMI圖像進(jìn)行裂縫識別,結(jié)果如圖10所示。蟻群算法將FMI圖像中黑色區(qū)域在一定范圍內(nèi)聚集為白點,然后連接圖像中白點的位置識別裂縫形態(tài)。CGAN識別得到的裂縫更連貫,裂縫形態(tài)更直觀。表1為CGAN與蟻群算法識別單張圖像耗時對比結(jié)果(其中,蟻群算法識別裂縫耗時研究來自文獻(xiàn)[16],CGAN識別裂縫耗時研究來自本文)。蟻群算法識別裂縫耗時與螞蟻數(shù)和迭代次數(shù)成正相關(guān)。螞蟻數(shù)和蟻群迭代次數(shù)太少時不能有效提取裂縫信息,因而限制了有效識別圖像中裂縫的最短耗時。在螞蟻數(shù)為3468,迭代次數(shù)為600時,耗時為28.74s/張。CGAN識別單張圖像過程耗時短,平均速度達(dá)到0.19s/張。

        圖10 采用本文方法和蟻群算法的裂縫識別結(jié)果a FMI圖像; b 本文方法識別裂縫結(jié)果; c 蟻群算法識別裂縫結(jié)果

        表1 CGAN與蟻群算法識別單張圖像耗時對比

        4 結(jié)論

        本文研究形成了CGAN自動識別裂縫技術(shù),研究結(jié)果表明:

        1) CGAN抗噪聲干擾能力強。即使裂縫圖像包含隨機噪聲,采用CGAN也可以準(zhǔn)確識別裂縫信息。

        2) CGAN識別FMI圖像低角度縫和水平裂縫準(zhǔn)確率高,實際成像測井圖像中裂縫識別準(zhǔn)確率為90%。識別圖像速度為0.19s/張,比蟻群算法識別裂縫速度有大幅度提高,且識別得到的裂縫形態(tài)更直觀。

        3) CGAN的缺點是對訓(xùn)練樣本和標(biāo)簽的依賴性較高。當(dāng)識別復(fù)雜形態(tài)裂縫時,設(shè)置的樣本標(biāo)簽對裂縫識別的干擾較大。各目標(biāo)訓(xùn)練樣本數(shù)量的均勻程度對識別結(jié)果也有影響。

        綜合計算機模擬裂縫識別結(jié)果、實際FMI圖像裂縫識別結(jié)果及與蟻群算法識別裂縫結(jié)果對比可知,CGAN是一種可以自動識別FMI圖像中裂縫的有效方法。目前CGAN只能識別裂縫,尚不能提取裂縫參數(shù)。下一步要完善訓(xùn)練過程中各形態(tài)的裂縫樣本,再結(jié)合其它深度學(xué)習(xí)方法,提取裂縫各項參數(shù),實現(xiàn)對裂縫的自動追蹤。

        猜你喜歡
        正弦曲線標(biāo)簽卷積
        基于3D-Winograd的快速卷積算法設(shè)計及FPGA實現(xiàn)
        基于FANUC系統(tǒng)的車削數(shù)控加工正弦曲線宏程序設(shè)計
        從濾波器理解卷積
        電子制作(2019年11期)2019-07-04 00:34:38
        圓環(huán)正弦曲線數(shù)控銑床編程加工探究
        無懼標(biāo)簽 Alfa Romeo Giulia 200HP
        車迷(2018年11期)2018-08-30 03:20:32
        不害怕撕掉標(biāo)簽的人,都活出了真正的漂亮
        海峽姐妹(2018年3期)2018-05-09 08:21:02
        基于傅里葉域卷積表示的目標(biāo)跟蹤算法
        基于宏程序的正弦曲線零件數(shù)控銑削加工
        標(biāo)簽化傷害了誰
        基于多進(jìn)制查詢樹的多標(biāo)簽識別方法
        計算機工程(2015年8期)2015-07-03 12:20:27
        久久综合久久综合久久| 伊人网在线视频观看| 亚洲乱在线播放| 九一精品少妇一区二区三区 | 国产一区不卡视频在线| 精品视频在线观看日韩| 天天躁日日躁狠狠躁| 国产91中文| 精品日产一区2区三区| 国产av天堂亚洲av刚刚碰| 无码欧美毛片一区二区三| 无码人妻一区二区三区在线视频| 欧美在线Aⅴ性色| 国产一区二区黄色的网站| 67194熟妇人妻欧美日韩| 99久久久无码国产精品9| 久久婷婷免费综合色啪| 开心久久婷婷综合中文字幕| 无码av一区二区大桥久未 | 色婷婷久久免费网站| 国产愉拍91九色国产愉拍| 挺进邻居丰满少妇的身体| 人人做人人妻人人精| 久草视频华人在线观看| 国产精品一区二区黄色| 天天色影网| 成年奭片免费观看视频天天看| 在线观看一区二区三区国产| 领导边摸边吃奶边做爽在线观看| 亚洲男人天堂| 在线观看视频日本一区二区三区| 亚洲国产精品中文字幕久久| 精品欧洲av无码一区二区| 亚洲伊人久久成人综合网| 精品一区2区3区4区| 国产av国片精品jk制服| 国产人妖视频一区二区| 亚洲一区日本一区二区| 欧美xxxxx高潮喷水| 女人大荫蒂毛茸茸视频| 亚洲国产精品成人久久av|