技術(shù)特征：

1.一種地質(zhì)剖面圖的生成方法，其特征在于，包括以下步驟：

基于地質(zhì)知識(shí)建立協(xié)同克立格方程組矩陣；

根據(jù)所述協(xié)同克立格方程組矩陣建立二維地質(zhì)剖面圖；

根據(jù)所述協(xié)同克立格方程組矩陣建立三維地質(zhì)模型剖面圖。

2.根據(jù)權(quán)利要求1所述的地質(zhì)剖面圖的生成方法，其特征在于，所述基于地質(zhì)知識(shí)建立協(xié)同克立格方程組矩陣具體步驟為：

設(shè)有N個(gè)變量構(gòu)成協(xié)同區(qū)域化變量集合{Z_k(x),k＝1,2,…N}；假設(shè)一階矩和二階矩存在且平穩(wěn)：

$E\left(Z_k\right(x\left)\right)=m_k,\∀x,\∀k---\left(1\right);$

其中：E(Z_k(x))為表示區(qū)域化變量Z_k(x)的數(shù)學(xué)期望；m_k為為常數(shù)；

設(shè)k₀為k＝0,1,2…N中某一數(shù)值，計(jì)算中心點(diǎn)在x₀處的承載上估計(jì)變量的平均值

$Z_{V_{k_0}}=\frac{1}{v_{k_0}}\underset{v_{k_0}}{\∫}{Z}_{k_0}\left(x\right)dx---\left(2\right);$

表示變量x₀的變化域；

滿足二階平穩(wěn)的協(xié)同區(qū)域化變量必須每一個(gè)區(qū)域化變量Z_k(x)的期望存在且平穩(wěn)，交叉協(xié)方差函數(shù)為：

C_kk′(h)＝E[Z_k′(x+h)·Z_k(x)] (3)；

根據(jù)權(quán)重系數(shù)和交叉協(xié)方差函數(shù)C_kk′(h)確定協(xié)同克立格方程組，具體公式如下：

$\left\{\begin{array}{c}{\mathrm{\Σ}}_{k^{\′}=1}^N{\mathrm{\Σ}}_{{\mathrm{\β}}_{k^{\′}}}^{n_{k^{\′}}}{\mathrm{\λ}}_{{\mathrm{\β}}_{k^{\′}}}{c}_{k^{\′}k}(v_{{\mathrm{\β}}_{k^{\′}}},v_{{\mathrm{\α}}_k})-{\mathrm{\μ}}_k=C_{k_{0}k}\left(V_{k_0^{\′}},v_{{\mathrm{\α}}_k}\right),\\ \∀{\mathrm{\α}}_k=1,2,...n_k;k^{\′}=1,2,...N,\\ {\mathrm{\Σ}}_{{\mathrm{\α}}_{k_0}=1}^{n_{k_0}}{\mathrm{\λ}}_{{\mathrm{\α}}_{k_0}}=1,\\ {\mathrm{\Σ}}_{{\mathrm{\α}}_{k=1}}^{n_k}{\mathrm{\λ}}_{{\mathrm{\α}}_k}=0,\∀k\≠k_0\end{array}\right.---\left(4\right);$

其中，β_k′為樣本標(biāo)號(hào)；n_k′為樣品個(gè)數(shù)；為權(quán)系數(shù)；為點(diǎn)承載；μ_k為拉格朗日參數(shù)；為協(xié)方差；為變量的變化域；

對(duì)公式(4)進(jìn)行線性變換，得到協(xié)同克立格方程組矩陣為：

kλ＝M (5)；

其中：

$\mathrm{\λ}=\left[\begin{array}{c}{\mathrm{\λ}}_{\left(1\right)\left(1\right)}\\ {\mathrm{\λ}}_{\left(1\right)\left(2\right)}\\ \begin{array}{c}.\\ .\\ .\end{array}\\ {\mathrm{\λ}}_{\left(1\right)\left(n_1\right)}\\ {\mathrm{\λ}}_{\left(2\right)\left(1\right)}\\ \begin{array}{c}.\\ .\\ .\end{array}\\ {\mathrm{\λ}}_{\left(n\right)\left(n_N\right)}\\ -{\mathrm{\μ}}_1\\ \begin{array}{c}.\\ .\\ .\end{array}\\ -{\mathrm{\μ}}_N\end{array}\right],M=\left[\begin{array}{c}C\left(v,v_{\left(1\right)\left(1\right)}\right)\\ C\left(v,v_{\left(1\right)\left(2\right)}\right)\\ \begin{array}{c}.\\ .\\ .\end{array}\\ C\left(v,v_{\left(1\right)\left(n_1\right)}\right)\\ C(v,v_{\left(2\right)\left(1\right)}\\ \begin{array}{c}.\\ .\\ .\end{array}\\ C(v,v_{\left(N\right)\left(n_N\right)}\\ 1\\ \begin{array}{c}.\\ .\\ .\end{array}\\ 0\end{array}\right],$

$K=\left[\begin{array}{cccccc}C(v_{\left(1\right)\left(1\right)},v_{\left(1\right)\left(1\right)})& ...& ...& C(v_{\left(1\right)\left(1\right)},v_{\left(N\right)\left(n_N\right)})& 1& 1\\ \begin{array}{c}.\\ .\\ .\end{array}& \begin{array}{c}.\\ .\\ .\end{array}& \begin{array}{c}.\\ .\\ .\end{array}& \begin{array}{c}.\\ .\\ .\end{array}& \begin{array}{c}.\\ .\\ .\end{array}& \begin{array}{c}.\\ .\\ .\end{array}\\ v_{\left(N\right)\left(n_N\right)}& ...& ...& C(v_{\left(N\right)\left(n_N\right)},v_{\left(N\right)\left(n_N\right)})& 1& 1\\ 1& 1& 0& 0& & \\ 0& 0& 1& 1& & \end{array}\right];$

其中，λ矩陣權(quán)系數(shù)向量；M矩陣為樣本數(shù)據(jù)和待觀察點(diǎn)的協(xié)方差矩陣；K矩陣為樣本數(shù)據(jù)間的協(xié)方差矩陣。

3.根據(jù)權(quán)利要求1所述的地質(zhì)剖面圖的生成方法，其特征在于，所述根據(jù)所述協(xié)同克立格方程組矩陣建立二維地質(zhì)剖面圖的具體步驟為：

基于地質(zhì)數(shù)據(jù)建立地質(zhì)屬性數(shù)據(jù)庫；

定義地質(zhì)對(duì)象之間的時(shí)空關(guān)系；

根據(jù)所述地質(zhì)屬性數(shù)據(jù)庫和所述時(shí)空關(guān)系采用協(xié)同克立格方程組矩陣構(gòu)建第一勢場；

根據(jù)二維地質(zhì)剖面圖的位置和所述第一勢場構(gòu)建第一地質(zhì)邊界線；

根據(jù)所述第一地質(zhì)邊界線的地質(zhì)屬性確定二維地質(zhì)剖面圖。

4.根據(jù)權(quán)利要求3所述的地質(zhì)剖面圖的生成方法，其特征在于，所述地質(zhì)數(shù)據(jù)包括：產(chǎn)狀、地層邊界、構(gòu)造、驗(yàn)證鉆孔、實(shí)測剖面和露頭。

5.根據(jù)權(quán)利要求3所述的地質(zhì)剖面圖的生成方法，其特征在于，所述第一勢場中任意兩點(diǎn)應(yīng)滿足以下公式：

$T^{*}\left(p\right)-T^{*}\left(p_0\right)={\mathrm{\Σ}}_{\mathrm{\α}=1}^M{\mathrm{\μ}}_{\mathrm{\α}}\left(T\right(p_{\mathrm{\α}})-T(p_{\mathrm{\α}}^{\′}\left)\right)+{\mathrm{\Σ}}_{\mathrm{\β}=1}^N{v}_{\mathrm{\β}}\frac{\∂T}{\∂u_{\mathrm{\β}}}\left(p_{\mathrm{\β}}\right)---\left(6\right);$

其中，T^*(p)表示任意點(diǎn)的勢值；T^*(p₀)表示已知點(diǎn)的勢值；M表示測量點(diǎn)的個(gè)數(shù)；N表示測量產(chǎn)狀的個(gè)數(shù)；μ_α表示測量點(diǎn)的權(quán)值；v_β表示產(chǎn)狀點(diǎn)的權(quán)值；T(p_α)表示序號(hào)為α的一個(gè)勢值；T(p′_α)表示序號(hào)為α的另一個(gè)勢值。

6.根據(jù)權(quán)利要求1所述的地質(zhì)剖面圖的生成方法，其特征在于，所述根據(jù)所述協(xié)同克立格方程組矩陣建立三維地質(zhì)模型剖面圖的具體步驟為：

采用網(wǎng)格索引與包圍盒碰撞檢測算法優(yōu)化歐式距離；

根據(jù)優(yōu)化后的所述歐式距離，采用協(xié)同克立格方程組矩陣構(gòu)建第二勢場；

根據(jù)所述第二勢場和三維地質(zhì)模型剖面的位置構(gòu)建第二地質(zhì)邊界線；

根據(jù)所述第二地質(zhì)邊界線的地質(zhì)屬性確定三維地質(zhì)模型剖面圖。

7.根據(jù)權(quán)利要求6所述的地質(zhì)剖面圖的生成方法，其特征在于，所述根據(jù)所述第二勢場和三維地質(zhì)剖面的位置構(gòu)建第二地質(zhì)邊界線具體包括：

針對(duì)實(shí)體模型時(shí)，根據(jù)所述第二勢場和三維地質(zhì)剖面的位置直接追蹤獲得第二地質(zhì)邊界線；

針對(duì)表面體模型時(shí)，先建三角網(wǎng)格，然后在所述三角網(wǎng)格上根據(jù)所述第二勢場和三維地質(zhì)剖面的位置追蹤獲得第二地質(zhì)邊界線。

8.根據(jù)權(quán)利要求7所述的地質(zhì)剖面圖的生成方法，其特征在于，所述實(shí)體模型包括：四面體模型或者六面體模型。

9.根據(jù)權(quán)利要求6所述的地質(zhì)剖面圖的生成方法，其特征在于，所述歐式距離的公式為：

$d=\frac{|{\mathrm{Ax}}_0+{\mathrm{By}}_0+{\mathrm{Cz}}_0+D|}{\sqrt{A^2+B^2+C^2}}---\left(7\right);$

其中，(x₀,y₀,z₀)為空間坐標(biāo)一點(diǎn)，A、B、C、D為系數(shù)。