本發(fā)明屬于地質(zhì)對(duì)象含油氣概率評(píng)價(jià)方法領(lǐng)域，具體地，涉及一種確定地質(zhì)風(fēng)險(xiǎn)空間分布的方法及系統(tǒng)。

背景技術(shù)：

隨著對(duì)油氣資源評(píng)價(jià)和勘探效益的要求不斷提高，同時(shí)為了充分利用勘探資料，作為資源評(píng)價(jià)重要組成部分的地質(zhì)風(fēng)險(xiǎn)空間分布引起了人們的重視，借助數(shù)學(xué)地質(zhì)方法和計(jì)算機(jī)技術(shù)，地質(zhì)風(fēng)險(xiǎn)空間分布的研究進(jìn)展迅速，出現(xiàn)了多種地質(zhì)風(fēng)險(xiǎn)空間分布預(yù)測(cè)方法。

目前地質(zhì)風(fēng)險(xiǎn)空間分布的預(yù)測(cè)方法主要有成因模型法、隨機(jī)過(guò)程的統(tǒng)計(jì)模型法和信息集成法三大類(lèi)。成因模型法主要考慮油氣成藏的地質(zhì)因素，適用于勘探初期地區(qū)；隨機(jī)過(guò)程的統(tǒng)計(jì)模型法主要以鉆探結(jié)果為基礎(chǔ)、以地質(zhì)統(tǒng)計(jì)方法為手段來(lái)預(yù)測(cè)油氣分布，適用于鉆探程度高、油氣地質(zhì)條件相對(duì)簡(jiǎn)單的地區(qū)；信息集成法與油氣空間分布有關(guān)的地質(zhì)、地球物理和勘探工程等信息為一體，能夠綜合反映油氣空間分布的特征。

技術(shù)實(shí)現(xiàn)要素：

然而，目前的地質(zhì)風(fēng)險(xiǎn)空間分布的預(yù)測(cè)方法，只給出了含油氣概率在不同地理位置的分布情況，并沒(méi)有進(jìn)一步與地理信息系統(tǒng)進(jìn)行結(jié)合，進(jìn)行地質(zhì)風(fēng)險(xiǎn)空間分布的地圖代數(shù)計(jì)算。對(duì)于一片地質(zhì)區(qū)域，影響其含油氣概率的因素有多個(gè)，對(duì)于區(qū)域中某個(gè)部分，需要對(duì)這些因素進(jìn)行綜合計(jì)算來(lái)得到此部分的含油氣概率預(yù)測(cè)情況。

另外，對(duì)于不同的地質(zhì)區(qū)域，含油氣概率相關(guān)因素對(duì)其影響的權(quán)重也不相同，為了能夠計(jì)算出不同地質(zhì)區(qū)域不同部分的含油氣概率分布情況，需要利用地圖代數(shù)計(jì)算功能，將每塊區(qū)域的含油氣概率因素進(jìn)行加權(quán)計(jì)算。在計(jì)算結(jié)果的展示上，也應(yīng)與地理信息系統(tǒng)技術(shù)相結(jié)合，從而將地質(zhì)風(fēng)險(xiǎn)的空間分布與區(qū)域的地 質(zhì)對(duì)象同時(shí)展現(xiàn)，更直觀(guān)地查看地質(zhì)風(fēng)險(xiǎn)的空間分布。

因此，本發(fā)明針對(duì)上述問(wèn)題，提供一種與地理信息系統(tǒng)結(jié)合更為緊密，具有更加完善的地圖代數(shù)功能，能夠更好展示含油氣概率在地理空間上分布的確定地質(zhì)風(fēng)險(xiǎn)空間分布的方法。

為了解決上述技術(shù)問(wèn)題，本申請(qǐng)的實(shí)施例首先提供了一種確定地質(zhì)風(fēng)險(xiǎn)空間分布的方法，包括：定義地質(zhì)風(fēng)險(xiǎn)評(píng)價(jià)所需的矢量圖層，其中，矢量圖層由一系列含有空間位置信息的圖元所組成，且每個(gè)圖元含有與之相應(yīng)的多項(xiàng)地質(zhì)風(fēng)險(xiǎn)評(píng)價(jià)值；對(duì)矢量圖層進(jìn)行柵格化處理得到柵格化矢量圖層；設(shè)置地圖代數(shù)公式，利用地圖代數(shù)公式對(duì)柵格化矢量圖層進(jìn)行地圖代數(shù)計(jì)算得到柵格化矢量圖層中每個(gè)格子的地質(zhì)風(fēng)險(xiǎn)評(píng)價(jià)值；基于柵格化矢量圖層中每個(gè)格子的地質(zhì)風(fēng)險(xiǎn)評(píng)價(jià)值進(jìn)行地質(zhì)風(fēng)險(xiǎn)空間分布的展示。

優(yōu)選地，在對(duì)矢量圖層進(jìn)行柵格化處理得到柵格化矢量圖層的步驟中，將矢量圖層的總區(qū)域劃分為網(wǎng)格形式，網(wǎng)格中每個(gè)格子的地質(zhì)風(fēng)險(xiǎn)評(píng)價(jià)值為落在這個(gè)格子內(nèi)的所有圖元的地質(zhì)風(fēng)險(xiǎn)評(píng)價(jià)值的均值。

優(yōu)選地，對(duì)于沒(méi)有圖元落在其中的格子，采用插值方法計(jì)算其地質(zhì)風(fēng)險(xiǎn)評(píng)價(jià)值，所述插值方法包括距離反平方方法和克里金插值方法。

優(yōu)選地，基于柵格化矢量圖層中每個(gè)格子的地質(zhì)風(fēng)險(xiǎn)評(píng)價(jià)值，采用等值線(xiàn)圖形進(jìn)行地質(zhì)風(fēng)險(xiǎn)空間分布的展示的步驟中。

優(yōu)選地，所述地質(zhì)風(fēng)險(xiǎn)評(píng)價(jià)值包括圈閉概率、充注概率、儲(chǔ)層概率和保存概率。

本申請(qǐng)的實(shí)施例還提供了一種確定地質(zhì)風(fēng)險(xiǎn)空間分布的系統(tǒng)，包括：定義模塊，其定義地質(zhì)風(fēng)險(xiǎn)評(píng)價(jià)所需的矢量圖層，其中，矢量圖層由一系列含有空間位置信息的圖元所組成，且每個(gè)圖元含有與之相應(yīng)的多項(xiàng)地質(zhì)風(fēng)險(xiǎn)評(píng)價(jià)值；柵格化模塊，其對(duì)矢量圖層進(jìn)行柵格化處理得到柵格化矢量圖層；計(jì)算模塊，其設(shè)置地圖代數(shù)公式，利用地圖代數(shù)公式對(duì)柵格化矢量圖層進(jìn)行地圖代數(shù)計(jì)算得到柵格化矢量圖層中每個(gè)格子的地質(zhì)風(fēng)險(xiǎn)評(píng)價(jià)值；展示模塊，其基于柵格化矢量圖層中每個(gè)格子的地質(zhì)風(fēng)險(xiǎn)評(píng)價(jià)值進(jìn)行地質(zhì)風(fēng)險(xiǎn)空間分布的展示。

優(yōu)選地，所述柵格化模塊將矢量圖層的總區(qū)域劃分為網(wǎng)格形式，網(wǎng)格中每個(gè)格子的地質(zhì)風(fēng)險(xiǎn)評(píng)價(jià)值為落在這個(gè)格子內(nèi)的所有圖元的地質(zhì)風(fēng)險(xiǎn)評(píng)價(jià)值的均值。

優(yōu)選地，對(duì)于沒(méi)有圖元落在其中的格子，所述柵格化模塊采用插值方法計(jì)算其地質(zhì)風(fēng)險(xiǎn)評(píng)價(jià)值，所述插值方法包括距離反平方方法和克里金插值方法。

優(yōu)選地，所述展示模塊基于柵格化矢量圖層中每個(gè)格子的地質(zhì)風(fēng)險(xiǎn)評(píng)價(jià)值，采用等值線(xiàn)圖形進(jìn)行地質(zhì)風(fēng)險(xiǎn)空間分布的展示的步驟中。

優(yōu)選地，所述地質(zhì)風(fēng)險(xiǎn)評(píng)價(jià)值包括圈閉概率、充注概率、儲(chǔ)層概率和保存概率。

與現(xiàn)有技術(shù)相比，上述方案中的一個(gè)或多個(gè)實(shí)施例可以具有如下優(yōu)點(diǎn)或有益效果。

本發(fā)明實(shí)施例的計(jì)算地質(zhì)風(fēng)險(xiǎn)空間分布的方法，將地質(zhì)風(fēng)險(xiǎn)的空間分布柵格化，與地理信息系統(tǒng)技術(shù)結(jié)合起來(lái)，直觀(guān)地展示了含油氣概率在地理空間上的分布情況，并采用插值方法對(duì)周邊地理區(qū)域進(jìn)行含油氣概率的近似估計(jì)，實(shí)現(xiàn)含油氣概率的預(yù)測(cè)，同時(shí)提供了完善的地圖代數(shù)公式設(shè)定功能，使得地質(zhì)風(fēng)險(xiǎn)空間分布的計(jì)算更加方便。

本發(fā)明的其它特征和優(yōu)點(diǎn)將在隨后的說(shuō)明書(shū)中闡述，并且，部分地從說(shuō)明書(shū)中變得顯而易見(jiàn)，或者通過(guò)實(shí)施本發(fā)明的技術(shù)方案而了解。本發(fā)明的目的和其他優(yōu)點(diǎn)可通過(guò)在說(shuō)明書(shū)、權(quán)利要求書(shū)以及附圖中所特別指出的結(jié)構(gòu)和/或流程來(lái)實(shí)現(xiàn)和獲得。

附圖說(shuō)明

附圖用來(lái)提供對(duì)本申請(qǐng)的技術(shù)方案或現(xiàn)有技術(shù)的進(jìn)一步理解，并且構(gòu)成說(shuō)明書(shū)的一部分。其中，表達(dá)本申請(qǐng)實(shí)施例的附圖與本申請(qǐng)的實(shí)施例一起用于解釋本申請(qǐng)的技術(shù)方案，但并不構(gòu)成對(duì)本申請(qǐng)技術(shù)方案的限制。

圖1是根據(jù)本發(fā)明實(shí)施例的確定地質(zhì)風(fēng)險(xiǎn)空間分布的方法流程圖；

圖2是根據(jù)本發(fā)明實(shí)施例的導(dǎo)入矢量圖層的顯示界面示意圖；

圖3是根據(jù)本發(fā)明實(shí)施例的輸入地圖代數(shù)的顯示界面示意圖；

圖4是根據(jù)本發(fā)明實(shí)施例的等值線(xiàn)展示地質(zhì)風(fēng)險(xiǎn)空間分布的顯示效果示意圖；

圖5是根據(jù)本發(fā)明實(shí)施例的確定地質(zhì)風(fēng)險(xiǎn)空間分布的系統(tǒng)結(jié)構(gòu)圖。

具體實(shí)施方式

以下將結(jié)合附圖及實(shí)施例來(lái)詳細(xì)說(shuō)明本發(fā)明的實(shí)施方式，借此對(duì)本發(fā)明如何應(yīng)用技術(shù)手段來(lái)解決技術(shù)問(wèn)題，并達(dá)成相應(yīng)技術(shù)效果的實(shí)現(xiàn)過(guò)程能充分理解并據(jù)以實(shí)施。本申請(qǐng)實(shí)施例以及實(shí)施例中的各個(gè)特征，在不相沖突前提下可以相互結(jié)合，所形成的技術(shù)方案均在本發(fā)明的保護(hù)范圍之內(nèi)。

另外，附圖的流程圖示出的步驟可以在諸如一組計(jì)算機(jī)可執(zhí)行指令的計(jì)算機(jī)系統(tǒng)中執(zhí)行。并且，雖然在流程圖中示出了邏輯順序，但是在某些情況下，可以以不同于此處的順序執(zhí)行所示出或描述的步驟。

(實(shí)施例一)

針對(duì)現(xiàn)有計(jì)算地質(zhì)風(fēng)險(xiǎn)空間分布的方法所存在的不足，本發(fā)明實(shí)施例提供了一種確定地質(zhì)風(fēng)險(xiǎn)空間分布的方法。概括來(lái)說(shuō)，本方法的計(jì)算流程就是導(dǎo)入地質(zhì)風(fēng)險(xiǎn)評(píng)價(jià)所需的矢量圖層，對(duì)導(dǎo)入的矢量圖層進(jìn)行柵格化處理，輸入地圖代數(shù)公式，解析地圖代數(shù)公式并利用地圖代數(shù)公式對(duì)柵格化處理后的矢量圖層進(jìn)行計(jì)算得到地質(zhì)風(fēng)險(xiǎn)空間分布，最后用等值線(xiàn)展示地質(zhì)風(fēng)險(xiǎn)空間分布。

圖1是根據(jù)本發(fā)明一實(shí)施例的確定地質(zhì)風(fēng)險(xiǎn)空間分布的方法流程圖。下面參考圖1來(lái)詳細(xì)說(shuō)明本方法的各個(gè)步驟。

在步驟S110中，定義地質(zhì)風(fēng)險(xiǎn)評(píng)價(jià)所需的矢量圖層。

需要說(shuō)明的是，矢量圖層由一系列含有空間位置信息的圖元所組成，圖元的類(lèi)型有點(diǎn)、線(xiàn)、多邊形、多點(diǎn)、多線(xiàn)、多多邊形六種，圖元具體的形狀和位置信息都是通過(guò)坐標(biāo)來(lái)存儲(chǔ)和表示的，通過(guò)這些坐標(biāo)中X和Y的最值可以得到圖元所占的范圍大小。如圖2中所示，本例中導(dǎo)入了3個(gè)圖層，分別為“川西J2SX控制儲(chǔ)量”，“川西J2SX控制儲(chǔ)量3”，“川西J2SX控制儲(chǔ)量4”，均為多邊形類(lèi)型的矢量圖層。

在地質(zhì)學(xué)中，對(duì)于用于地質(zhì)風(fēng)險(xiǎn)評(píng)價(jià)的矢量圖層，矢量圖層代表一片地理區(qū)域，矢量圖層中每個(gè)圖元代表這片區(qū)域中一個(gè)地質(zhì)對(duì)象，并在其中存儲(chǔ)了這個(gè)地質(zhì)對(duì)象地理上的形狀和位置坐標(biāo)信息。

本例中三個(gè)矢量圖層為川西坳陷沙溪廟組的勘探程度圖，圖2中不同顏色深度的圖層分別代表了這三個(gè)矢量圖層對(duì)應(yīng)的地理區(qū)域的形狀、位置和范圍。

在定義地質(zhì)風(fēng)險(xiǎn)評(píng)價(jià)所需的矢量圖層時(shí)，會(huì)選取若干個(gè)矢量圖層文件及其屬性文件，矢量圖層文件中存儲(chǔ)了一片區(qū)域的所有地質(zhì)對(duì)象的形狀和位置信息，屬性文件中存儲(chǔ)了區(qū)域中各個(gè)地質(zhì)對(duì)象的配套屬性。每項(xiàng)配套屬性是一個(gè)數(shù)據(jù)庫(kù)字 段，對(duì)于用于地質(zhì)風(fēng)險(xiǎn)評(píng)價(jià)的矢量圖層，其屬性文件中應(yīng)有評(píng)價(jià)其地質(zhì)風(fēng)險(xiǎn)的相關(guān)字段。地質(zhì)風(fēng)險(xiǎn)評(píng)價(jià)即含油氣概率的評(píng)價(jià)，對(duì)于用于地質(zhì)風(fēng)險(xiǎn)評(píng)價(jià)的矢量圖層，圖層中每個(gè)圖元會(huì)含有與之相應(yīng)的多項(xiàng)地質(zhì)風(fēng)險(xiǎn)評(píng)價(jià)值作為屬性信息，代表這個(gè)圖元對(duì)應(yīng)的地質(zhì)對(duì)象所在區(qū)域含油氣概率的大小，可以利用這些值來(lái)得到這個(gè)矢量圖層所對(duì)應(yīng)區(qū)域的含油氣概率分布情況。

在本實(shí)施例中，對(duì)導(dǎo)入的三個(gè)矢量圖層添加其配套的屬性數(shù)據(jù)，使得圖層中不同圖元有不同的地質(zhì)風(fēng)險(xiǎn)評(píng)價(jià)值，可以利用這些值來(lái)得到這個(gè)矢量圖層所對(duì)應(yīng)區(qū)域的含油氣概率分布情況。如圖2中下方輸出空間中所示，給出了“川西J2SX控制儲(chǔ)量”這一圖層對(duì)應(yīng)的屬性數(shù)據(jù)，其中包含五個(gè)字段，分別為“geometryid”，“圈閉概率”，“充注概率”，“儲(chǔ)層概率”，“保存概率”，第一個(gè)字段為圖元的編號(hào)，后面4個(gè)字段代表此圖元的4項(xiàng)地質(zhì)風(fēng)險(xiǎn)評(píng)價(jià)參數(shù)，圖元的含油氣概率為這4個(gè)參數(shù)的加權(quán)平均。

在步驟S120中，對(duì)矢量圖層進(jìn)行柵格化處理得到柵格化矢量圖層。

具體地，在對(duì)導(dǎo)入的矢量圖層進(jìn)行柵格化處理時(shí)，需要將所有矢量圖層的總區(qū)域按X和Y兩個(gè)方向進(jìn)行切分，劃分為網(wǎng)格形式，網(wǎng)格中每個(gè)格子的地質(zhì)風(fēng)險(xiǎn)評(píng)價(jià)值為落在這個(gè)格子內(nèi)的所有圖元的地質(zhì)風(fēng)險(xiǎn)評(píng)價(jià)值的均值。

所有圖層都按照上述同樣的劃分方式進(jìn)行劃分，由于矢量圖層中圖元的形狀是不規(guī)則的，將其網(wǎng)格化就可以得到規(guī)則的區(qū)域所對(duì)應(yīng)的地質(zhì)風(fēng)險(xiǎn)評(píng)價(jià)值，使得圖層與圖層之間可以進(jìn)行地圖代數(shù)計(jì)算及地質(zhì)風(fēng)險(xiǎn)空間分布的展示。例如對(duì)于“川西J2SX控制儲(chǔ)量”圖層，其柵格化后的一個(gè)格子中的“圈閉概率”的值為落在這個(gè)格子中的所有圖元的“圈閉概率”的值的平均數(shù)。

對(duì)于沒(méi)有圖元落在其中的格子，其地質(zhì)風(fēng)險(xiǎn)評(píng)價(jià)值用插值的方法進(jìn)行計(jì)算，插值方法分為距離反平方和克里金插值兩種，用戶(hù)可以自行選擇。

距離反平方方法是根據(jù)格子周?chē)鷪D元的地質(zhì)風(fēng)險(xiǎn)評(píng)價(jià)值，加權(quán)平均計(jì)算出格子的地質(zhì)風(fēng)險(xiǎn)評(píng)價(jià)值，其中權(quán)重是圖元到格子距離的平方的倒數(shù)。

克里金插值又稱(chēng)空間自協(xié)方差最佳插值法，其廣泛應(yīng)用于地下水模擬，土壤制圖等領(lǐng)域，是一種很有用的地質(zhì)統(tǒng)計(jì)網(wǎng)格化方法，關(guān)于其具體的算法實(shí)現(xiàn)已有較多資料介紹，在此不再贅述。本實(shí)施例中所采用的為克里金插值，對(duì)沒(méi)有圖元落在其中的格子給出插值計(jì)算得到的地質(zhì)風(fēng)險(xiǎn)評(píng)價(jià)相關(guān)參數(shù)值。

在本步驟中，由于將與地質(zhì)風(fēng)險(xiǎn)的空間分布相關(guān)的矢量圖層進(jìn)行了柵格化， 與地理信息系統(tǒng)技術(shù)結(jié)合起來(lái)，因此能夠直觀(guān)地展示含油氣概率在地理空間上的分布情況。另外，采用插值方法對(duì)周邊地理區(qū)域進(jìn)行含油氣概率的近似估計(jì)，實(shí)現(xiàn)了含油氣概率的預(yù)測(cè)。

在步驟S130中，設(shè)置地圖代數(shù)公式，利用地圖代數(shù)公式對(duì)柵格化矢量圖層進(jìn)行地圖代數(shù)計(jì)算得到柵格化矢量圖層中每個(gè)格子的地質(zhì)風(fēng)險(xiǎn)評(píng)價(jià)值。

需要說(shuō)明的是，地圖代數(shù)計(jì)算即利用用戶(hù)自定義的公式對(duì)柵格化后的各個(gè)圖層進(jìn)行計(jì)算，計(jì)算結(jié)果也是柵格化的矢量圖層，區(qū)域及格子的劃分都與之前未經(jīng)計(jì)算的矢量圖層相同。計(jì)算結(jié)果中每個(gè)格子的地質(zhì)風(fēng)險(xiǎn)評(píng)價(jià)值為各個(gè)矢量圖層的格子的地質(zhì)風(fēng)險(xiǎn)評(píng)價(jià)值經(jīng)過(guò)公式運(yùn)算得到的值。

在地圖代數(shù)公式中，最簡(jiǎn)單的如兩個(gè)圖層相加，表示將這兩個(gè)圖層中的地質(zhì)風(fēng)險(xiǎn)評(píng)價(jià)值相加，計(jì)算結(jié)果中每個(gè)格子的地質(zhì)風(fēng)險(xiǎn)評(píng)價(jià)值為這兩個(gè)圖層在此格子對(duì)應(yīng)的地質(zhì)風(fēng)險(xiǎn)評(píng)價(jià)值之和。

如圖3所示，本實(shí)施例中設(shè)置并輸入的地圖代數(shù)公式為：

{[Layer(川西J2SX控制儲(chǔ)量.圈閉概率)*0.3+Layer(川西J2SX控制儲(chǔ)量.充注概率)*0.3+Layer(川西J2SX控制儲(chǔ)量.儲(chǔ)層概率)*0.2+Layer(川西J2SX控制儲(chǔ)量.保存概率)*0.2]+[Layer(川西J2SX控制儲(chǔ)量3.圈閉概率)*0.25+Layer(川西J2SX控制儲(chǔ)量3.充注概率)*0.25+Layer(川西J2SX控制儲(chǔ)量3.儲(chǔ)層概率)*0.3+Layer(川西J2SX控制儲(chǔ)量3.保存概率)*0.2]+[Layer(川西J2SX控制儲(chǔ)量4.圈閉概率)*0.25+Layer(川西J2SX控制儲(chǔ)量4.充注概率)*0.3+Layer(川西J2SX控制儲(chǔ)量4.儲(chǔ)層概率)*0.3+Layer(川西J2SX控制儲(chǔ)量4.保存概率)*0.15]}/3。

其含義為，首先對(duì)每個(gè)柵格化矢量圖層計(jì)算其地質(zhì)風(fēng)險(xiǎn)評(píng)價(jià)分布，方法為對(duì)地質(zhì)風(fēng)險(xiǎn)評(píng)價(jià)相關(guān)的4個(gè)字段的值進(jìn)行加權(quán)平均，公式中不同圖層的4個(gè)字段的權(quán)重有所不同，如“川西J2SX控制儲(chǔ)量”4個(gè)字段的權(quán)重分別為0.3、0.3、0.2、0.2，而“川西J2SX控制儲(chǔ)量3”4個(gè)字段的權(quán)重分別為0.25、0.25、0.3、0.2表示不同評(píng)價(jià)參數(shù)對(duì)于不同的地理區(qū)域?qū)ο笤诤蜌飧怕噬嫌绊懥^(qū)別。得到每個(gè)圖層的地質(zhì)風(fēng)險(xiǎn)評(píng)價(jià)分布后，將這些圖層的地質(zhì)風(fēng)險(xiǎn)評(píng)價(jià)值取均值，得到了總體區(qū)域的地質(zhì)風(fēng)險(xiǎn)評(píng)價(jià)分布。

另外，如圖3所示，本實(shí)施例中網(wǎng)格的劃分方式為10行×10列，表示將矢量圖層按照10×10的方式劃分來(lái)實(shí)現(xiàn)柵格化。

為了更好地利用設(shè)置的地圖代數(shù)公式，在計(jì)算機(jī)處理的過(guò)程中，計(jì)算機(jī)會(huì)對(duì) 用戶(hù)輸入的公式進(jìn)行解析。在本實(shí)施例中，公式中的“LAYER(圖層名.字段名)”表示一個(gè)柵格化矢量圖層的字段，運(yùn)算符除了四則運(yùn)算外，還有指數(shù)、對(duì)數(shù)、三角函數(shù)等等，都可以設(shè)置到公式中，公式在解譯過(guò)程中會(huì)按照這些運(yùn)算符的優(yōu)先級(jí)，逐步解譯。

對(duì)于“LAYER(圖層名.字段名)”這樣的公式片段，會(huì)解譯為具體的柵格化矢量圖層；對(duì)于數(shù)值型的公式片段，也會(huì)解譯為柵格化矢量圖層，其中每個(gè)格子的值即片段代表的數(shù)值，例如本實(shí)施例公式中含有數(shù)值0.3，則會(huì)解譯為一個(gè)柵格化矢量圖層，其網(wǎng)格劃分方式與圖層的劃分方式相同，都為10×10的網(wǎng)格，每個(gè)格子對(duì)應(yīng)的屬性值均為0.3。

柵格化矢量圖層的計(jì)算是對(duì)柵格化矢量圖層每個(gè)格子的地質(zhì)風(fēng)險(xiǎn)評(píng)價(jià)值進(jìn)行計(jì)算。如“l(fā)og(LAYER(圖層名.字段名))”的計(jì)算結(jié)果為與這個(gè)柵格化矢量圖層范圍及格子劃分相同的柵格化矢量圖層，不同的是其中每個(gè)格子的地質(zhì)風(fēng)險(xiǎn)評(píng)價(jià)值為之前的值取對(duì)數(shù)；“l(fā)og(LAYER(圖層名.字段名))*0.3”的計(jì)算結(jié)果為與這個(gè)柵格化矢量圖層范圍及格子劃分相同的柵格化矢量圖層，不同的是其中每個(gè)格子的地質(zhì)風(fēng)險(xiǎn)評(píng)價(jià)值為之前的0.3倍；“l(fā)og(LAYER(圖層名1.字段名1))+log(LAYER(圖層名2.字段名2))”的計(jì)算結(jié)果為與這個(gè)兩個(gè)柵格化矢量圖層范圍及格子劃分相同的柵格化矢量圖層，其中每個(gè)格子的地質(zhì)風(fēng)險(xiǎn)評(píng)價(jià)值為這個(gè)兩個(gè)柵格化矢量圖層對(duì)應(yīng)格子的地質(zhì)風(fēng)險(xiǎn)評(píng)價(jià)值之和。經(jīng)過(guò)用戶(hù)自定義公式的計(jì)算，最終得到的柵格化矢量圖層中每個(gè)格子的地質(zhì)風(fēng)險(xiǎn)評(píng)價(jià)值即經(jīng)過(guò)同樣公式運(yùn)算得到的結(jié)果。

在本實(shí)施例中，計(jì)算結(jié)果為一個(gè)柵格化矢量圖層，圖層的網(wǎng)格劃分方式為10行×10列，對(duì)于每個(gè)格子，其格子對(duì)應(yīng)的屬性值為各個(gè)圖層在此格子的對(duì)應(yīng)的屬性值經(jīng)過(guò)實(shí)施例中設(shè)定的公式計(jì)算得到的結(jié)果。

由于本步驟提供了完善的地圖代數(shù)公式設(shè)定功能，使得地質(zhì)風(fēng)險(xiǎn)空間分布的計(jì)算更加方便。

另外，為了更直觀(guān)地看到計(jì)算得到的地質(zhì)風(fēng)險(xiǎn)空間分布，本實(shí)施例還包括步驟S140。

在步驟S140中，基于柵格化矢量圖層中每個(gè)格子的地質(zhì)風(fēng)險(xiǎn)評(píng)價(jià)值進(jìn)行地質(zhì)風(fēng)險(xiǎn)空間分布的展示。

具體地，在本步驟中，對(duì)于計(jì)算得到的柵格化矢量圖層，用等值線(xiàn)圖形來(lái)展 示這個(gè)圖層的地質(zhì)風(fēng)險(xiǎn)空間分布。等值線(xiàn)圖層傳入的信息為各個(gè)點(diǎn)及各個(gè)點(diǎn)對(duì)應(yīng)的值，即計(jì)算結(jié)果中各個(gè)網(wǎng)格的中心點(diǎn)坐標(biāo)及各個(gè)網(wǎng)格的地質(zhì)風(fēng)險(xiǎn)評(píng)價(jià)值。

如圖4所示，展示了地質(zhì)風(fēng)險(xiǎn)空間分布的結(jié)果。圖中不同位置有不同的色彩深淺度，右側(cè)的柱狀圖例展示了顏色深淺度和數(shù)值的對(duì)應(yīng)關(guān)系，數(shù)值即計(jì)算得到的地質(zhì)風(fēng)險(xiǎn)評(píng)價(jià)值?？梢钥闯鰷\色代表的值比較大，即相應(yīng)位置的地質(zhì)風(fēng)險(xiǎn)評(píng)價(jià)值較高，深色代表的值比較小，表示相應(yīng)位置的地質(zhì)風(fēng)險(xiǎn)評(píng)價(jià)值較低。

至此，地質(zhì)風(fēng)險(xiǎn)空間分布的計(jì)算及展示工作進(jìn)行完畢。

綜上所述，本發(fā)明實(shí)施例的計(jì)算地質(zhì)風(fēng)險(xiǎn)空間分布的方法，將地質(zhì)風(fēng)險(xiǎn)的空間分布柵格化，與地理信息系統(tǒng)技術(shù)結(jié)合起來(lái)，直觀(guān)地展示了含油氣概率在地理空間上的分布情況，并采用插值方法對(duì)周邊地理區(qū)域進(jìn)行含油氣概率的近似估計(jì)，實(shí)現(xiàn)含油氣概率的預(yù)測(cè)，同時(shí)提供了完善的地圖代數(shù)公式設(shè)定功能，使得地質(zhì)風(fēng)險(xiǎn)空間分布的計(jì)算更加方便。

(實(shí)施例二)

圖5是根據(jù)本發(fā)明實(shí)施例的確定地質(zhì)風(fēng)險(xiǎn)空間分布的系統(tǒng)結(jié)構(gòu)圖。下面參考圖5來(lái)說(shuō)明本系統(tǒng)的各部分組成。

如圖5所示，該系統(tǒng)包括定義模塊51、柵格化模塊52、計(jì)算模塊53和展示模塊54。定義模塊51、柵格化模塊52、計(jì)算模塊53和展示模塊54分別執(zhí)行實(shí)施例一的步驟S110、步驟S120、步驟S130和步驟S140的操作，在此不再詳細(xì)展開(kāi)。

本領(lǐng)域的技術(shù)人員應(yīng)該明白，上述的本發(fā)明實(shí)施例所提供的裝置和/或系統(tǒng)的各組成部分，以及方法中的各步驟，它們可以集中在單個(gè)的計(jì)算裝置上，或者分布在多個(gè)計(jì)算裝置所組成的網(wǎng)絡(luò)上。可選地，它們可以用計(jì)算裝置可執(zhí)行的程序代碼來(lái)實(shí)現(xiàn)。從而，可以將它們存儲(chǔ)在存儲(chǔ)裝置中由計(jì)算裝置來(lái)執(zhí)行，或者將它們分別制作成各個(gè)集成電路模塊，或者將它們中的多個(gè)模塊或步驟制作成單個(gè)集成電路模塊來(lái)實(shí)現(xiàn)。這樣，本發(fā)明不限制于任何特定的硬件和軟件結(jié)合。

雖然本發(fā)明所揭露的實(shí)施方式如上，但所述的內(nèi)容僅為便于理解本發(fā)明技術(shù)方案而采用的實(shí)施方式，并非用以限定本發(fā)明。任何本發(fā)明所屬領(lǐng)域內(nèi)的技術(shù)人員，在不脫離本發(fā)明所揭露的精神和范圍的前提下，可以在實(shí)施的形式及細(xì)節(jié)上進(jìn)行任何的修改與變化，但本發(fā)明的專(zhuān)利保護(hù)范圍，仍須以所附的權(quán)利要求書(shū)所 界定的范圍為準(zhǔn)。