轉(zhuǎn)：公眾號(hào)，專業(yè)網(wǎng)站原創(chuàng)：劉彥鋒智能油田開發(fā)

前一陣去捷克參加了數(shù)學(xué)地質(zhì)年會(huì)（IAMG2018）,聽完匯報(bào)，發(fā)現(xiàn)在數(shù)學(xué)地質(zhì)領(lǐng)域用python的越來越多了，他們的很多成果都是用python做出來的。不管喜不喜歡，這個(gè)語言已經(jīng)成了主流，特別是在人工智能領(lǐng)域，幾乎所有的框架都有python接口，用起來很方便。秉承“人生苦短，我用python”的理念，也踏上了python學(xué)習(xí)之路，盡管感覺C#比python更優(yōu)雅。

學(xué)習(xí)python的目的是幫助油氣開發(fā)的智能化，特別是智能油田管理相關(guān)領(lǐng)域，為了避免重復(fù)造輪子，有必要了解目前有哪些開源的python包可以用。下面從數(shù)值巖心、測(cè)井解釋、地質(zhì)建模、地震處理解釋、數(shù)值模擬、裂縫模擬幾個(gè)方面介紹。

1、數(shù)字巖心

數(shù)字巖心是表征巖石微觀動(dòng)態(tài)和靜態(tài)特征的關(guān)鍵，對(duì)掃描得到的數(shù)字巖心進(jìn)行建模和各種模擬需要巖心重建，孔隙網(wǎng)絡(luò)模型是常見的巖心重建方法，該方法的python包有OpenPNM和Poreflow等。

http://www.gslib.com/

https://github.com/PMEAL

https://github.com/pmeal/openpnm

https://github.com/ggorman/poreflow

2測(cè)井解釋

包括測(cè)井?dāng)?shù)據(jù)的讀寫、各類解釋方法、圖形顯示等，比如標(biāo)準(zhǔn)的las格式數(shù)據(jù)讀寫的python包lasio，各種深度轉(zhuǎn)換的包las，各種濾波處理的logging，基本的測(cè)井解釋，用excel和簡(jiǎn)單的幾行python腳本就可以，復(fù)雜的解釋可以借助專門的解釋包，比如petropy等。

https://lasio.readthedocs.io/en/latest/

https://github.com/petroGG/Basic-Well-Log-Interpretation

https://github.com/petroGG/Basic-Well-Log-Interpretation/blob/master/Basic%20well%20log%20interpretation.ipynb

http://www.cwls.org/las/

https://github.com/WarrenWeckesser/las

https://github.com/WalBytes/PyLogFinder_Package

https://github.com/toddheitmann/petropy

3地質(zhì)建模

地質(zhì)建模方面，主要是地質(zhì)統(tǒng)計(jì)學(xué)方面，以前的Gslib是地質(zhì)統(tǒng)計(jì)學(xué)中使用最廣泛的fortran開源代碼，目前也已經(jīng)封裝成為python包，有PyGSLIB，GeostatsPy，除此之外還有基于gslib的擁有更高級(jí)接口的Pygeostat，pyGeoStatistics，也有不依賴GSLIB python包，比如pykriging，krige，geostatsmodels，HPGL，scikit-gstat，利用這些python包可以進(jìn)行變差函數(shù)求取，各種克里金插值，序貫高斯模擬、序貫指示模擬等。近兩年來大名鼎鼎的MichaelPyrcz（代表作是Geostatistical Reservoir Modeling2ndEdition）在twitter上一直很活躍地在分享它從GSLIB轉(zhuǎn)向python的學(xué)習(xí)過程(@GeostatsGuy)。

https://github.com/opengeostat/pygslib

https://github.com/GeostatsGuy/GeostatsPy

https://opengeostat.github.io/pygslib/

http://www.ccgalberta.com/pygeostat/index.html

https://github.com/whimian/pyGeoStatistics

https://github.com/cjohnson318/geostatsmodels/

https://github.com/hpgl/hpgl

http://hpgl.mit-ufa.com/

https://github.com/mmaelicke/scikit-gstat/tree/v0.2.1

http://pykriging.com/

https://twitter.com/GeostatsGuy

下圖是測(cè)試kriging包的代碼及運(yùn)算結(jié)果

4、地震資料解釋處理

地震資料讀寫、處理、解釋和反演方面也有不少包，比如地震資料讀寫的segpy、segyio，巖石物理的pyrocko，合成地震記錄的ricker，處理解釋的ObsPy，屬性解釋的interpies、geoprobe、d2geo，地震反演與成像的pysit，nessi，時(shí)深轉(zhuǎn)換的TotalDepth，地震層次成像的pytomo3d，地震數(shù)據(jù)3D顯示的seisplot等，以及三維解釋Java包3d-painting。

https://github.com/sixty-north/segpy

https://github.com/Statoil/segyio

https://pyrocko.org/

https://github.com/pyrocko

https://github.com/lijuhttps://github.com/computational-seismology/pytomo3dnzh/ricker

https://github.com/obspy/obspy/wiki

https://github.com/jobar8/interpies

https://github.com/joferkington/python-geoprobe

https://github.com/dfitzgerald3/d2geo

http://pysit.org/

https://github.com/pysit/pysit

https://github.com/PageotD/nessi

https://github.com/paulross/TotalDepth

https://github.com/computational-seismology/pytomo3d

https://github.com/agile-geosci培養(yǎng)桃紅ence/seisplot

https://github.com/chrisengelsma/3d-painting

5、數(shù)值模擬

油藏?cái)?shù)值模擬的難度相對(duì)較大，相應(yīng)的包也比較少，能找到的有單相裂縫-孔隙介質(zhì)流動(dòng)模擬的PorePy，單相2D模擬的SimplePetro，單相3D模擬的PRESTO。

https://pypi.org/project/porepy/#porepy-features

https://github.com/ishovkun/SimplePetro

https://github.com/padmec-reservoir/PRESTO

6、裂縫表征

裂縫特征統(tǒng)計(jì)的geo_frac_analysis，裂縫生成模擬的fracture。

https://github.com/TerminusEst/geo_frac_analysis

https://github.com/inconvergent/fracture

總結(jié)

以上是生活随笔為你收集整理的【油藏描述】分享一波油藏描述方面的Python开源包（转）的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。