[發(fā)明專利]一種含水層結(jié)構(gòu)變異轉(zhuǎn)移概率的地質(zhì)統(tǒng)計模擬方法及裝置

有效

• 申請?zhí)枺?/td>: 202010760640.8
• 申請日: 2020-07-31
• 公開（公告）號: CN111859709B
• 公開（公告）日: 2022-10-18
• 發(fā)明（設計）人: 代鋒剛;張發(fā)旺;王濱;么紅超
• 申請（專利權(quán)）人: 河北地質(zhì)大學;中國地質(zhì)調(diào)查局水文地質(zhì)環(huán)境地質(zhì)調(diào)查中心
• 主分類號: G06F30/20
• 分類號: G06F30/20;G06T17/05;G06F111/08
• 代理公司: 石家莊國為知識產(chǎn)權(quán)事務所 13120
• 代理人: 魏笑

1.一種含水層結(jié)構(gòu)變異轉(zhuǎn)移概率的地質(zhì)統(tǒng)計模擬方法，其特征在于，包括：

根據(jù)目標區(qū)域的地層巖性數(shù)據(jù)計算所述目標區(qū)域地質(zhì)類型的垂向離散馬爾可夫鏈模型；

根據(jù)連續(xù)性馬爾可夫鏈對所述垂向離散馬爾可夫鏈模型進行擬合，得到垂向連續(xù)馬爾可夫鏈模型；

根據(jù)所述垂向連續(xù)馬爾可夫鏈模型計算側(cè)向連續(xù)馬爾可夫鏈模型；

基于所述垂向連續(xù)馬爾可夫鏈模型和所述側(cè)向連續(xù)馬爾可夫鏈模型，計算三維馬爾可夫鏈模型；

根據(jù)所述三維馬爾可夫鏈模型獲得所述目標區(qū)域的含水層結(jié)構(gòu)變異的轉(zhuǎn)移概率的地質(zhì)統(tǒng)計模型；

在所述根據(jù)目標區(qū)域的地層巖性數(shù)據(jù)計算所述目標區(qū)域地質(zhì)類型的垂向離散馬爾可夫鏈模型之前，所述方法還包括：

獲取對所述目標區(qū)域進行地質(zhì)鉆孔實驗后得到的初始地層巖性數(shù)據(jù)；所述初始地層巖性數(shù)據(jù)包括：鉆孔深度、地層厚度數(shù)據(jù)、地層分層巖性數(shù)據(jù)以及地層滲透系數(shù)數(shù)據(jù)；

去除所述初始地層巖性數(shù)據(jù)中的噪聲和異常數(shù)據(jù)，得到第二地層巖性數(shù)據(jù)；

補全所述目標區(qū)域的第二地層巖性數(shù)據(jù)，得到所述目標區(qū)域的地層巖性數(shù)據(jù)；

其中，使用數(shù)據(jù)驅(qū)動方法對所述初始地層巖性數(shù)據(jù)進行異常判定，所述數(shù)據(jù)驅(qū)動方法包括根據(jù)初始地層巖性數(shù)據(jù)的歷史數(shù)據(jù)建立訓練模型的數(shù)據(jù)測試集合；

在信息不完全的情況下，利用已有的所述第二地層巖性數(shù)據(jù)進行模擬預測，補全所述第二地層巖性數(shù)據(jù)中缺失的部分，得到目標區(qū)域的地層巖性數(shù)據(jù)。

2.如權(quán)利要求1所述的一種含水層結(jié)構(gòu)變異轉(zhuǎn)移概率的地質(zhì)統(tǒng)計模擬方法，其特征在于，所述根據(jù)目標區(qū)域的地層巖性數(shù)據(jù)計算所述目標區(qū)域地質(zhì)類型的垂向離散馬爾可夫鏈模型，包括：

根據(jù)預設距離對所述地層巖性數(shù)據(jù)進行垂向上的劃分，得到垂向上的離散地層巖性數(shù)據(jù)；

根據(jù)所述垂向上的離散地層巖性數(shù)據(jù)計算所述垂向離散馬爾可夫鏈模型。

3.如權(quán)利要求2所述的一種含水層結(jié)構(gòu)變異轉(zhuǎn)移概率的地質(zhì)統(tǒng)計模擬方法，其特征在于，所述根據(jù)所述垂向上的離散地層巖性數(shù)據(jù)計算所述垂向離散馬爾可夫鏈模型，包括：

根據(jù)所述垂向上的離散地層巖性數(shù)據(jù)計算第一地質(zhì)類型向第二地質(zhì)類型轉(zhuǎn)移的轉(zhuǎn)移概率，并得到轉(zhuǎn)移概率矩陣；所述第一地質(zhì)類型和所述第二地質(zhì)類型均為任一地質(zhì)類型；

根據(jù)所述轉(zhuǎn)移概率矩陣計算所述垂向離散馬爾可夫鏈模型。

4.如權(quán)利要求1所述的一種含水層結(jié)構(gòu)變異轉(zhuǎn)移概率的地質(zhì)統(tǒng)計模擬方法，其特征在于，所述根據(jù)連續(xù)性馬爾可夫鏈對所述垂向離散馬爾可夫鏈模型進行擬合，得到垂向連續(xù)馬爾可夫鏈模型，包括：

根據(jù)所述垂向離散馬爾可夫鏈模型計算矩陣指數(shù)擬合曲線；

根據(jù)所述矩陣指數(shù)擬合曲線對所述垂向離散馬爾可夫鏈模型進行擬合，得到所述垂向連續(xù)馬爾可夫鏈模型。

5.如權(quán)利要求1所述的一種含水層結(jié)構(gòu)變異轉(zhuǎn)移概率的地質(zhì)統(tǒng)計模擬方法，其特征在于，所述根據(jù)所述垂向連續(xù)馬爾可夫鏈模型計算側(cè)向連續(xù)馬爾可夫鏈模型，包括：

根據(jù)所述地層巖性數(shù)據(jù)確定所述目標區(qū)域內(nèi)第一地質(zhì)類型在側(cè)向和垂向上的平均長度比例；所述第一地質(zhì)類型為所述目標區(qū)域內(nèi)任一地質(zhì)類型；

根據(jù)所述目標區(qū)域內(nèi)各個地質(zhì)類型對應的平均長度比例和所述垂向連續(xù)馬爾可夫鏈模型，計算所述側(cè)向連續(xù)馬爾可夫鏈模型。