(19)國家知識產(chǎn)權(quán)局(12)發(fā)明專利(65)同一申請的已公布的文獻號道163號賀婉慧高子涵(74)專利代理機構(gòu)南京佰騰智信知識產(chǎn)權(quán)代理事務(wù)所(普通合伙)32509專利代理師呂丹一種基于離散元的逆沖推覆體系下三維斷層/裂縫網(wǎng)絡(luò)建模方法及裝置本發(fā)明公開了一種基于離散元的逆沖推覆體系下三維斷層/裂縫網(wǎng)絡(luò)建模方法及裝置，涉及三維斷層建模的技術(shù)領(lǐng)域，包括基于目標(biāo)區(qū)域的地震測井資料參數(shù)，構(gòu)建初始化斷層面模型；提取目標(biāo)區(qū)域的地層參數(shù)，結(jié)合平衡剖面恢復(fù)方法，確定變形過程；根據(jù)地層參數(shù)標(biāo)定離散元的細(xì)觀參數(shù)，沉積第一離散元模型，優(yōu)化生成第二離散元模型，設(shè)置邊界條件進行逆沖推覆模擬，并通過第二離散元模型的輸出結(jié)果與實際地質(zhì)三維斷層和裂縫網(wǎng)絡(luò)，分析和建立目標(biāo)區(qū)域的空間展布特征，并對三維斷層的活動性進行動態(tài)評估與趨勢預(yù)測。本發(fā)明融合地質(zhì)力學(xué)模型，克服并通過第二離散元模型的輸出結(jié)果與實際地質(zhì)剖面的對比，21.一種基于離散元的逆沖推覆體系下三維斷層/裂縫網(wǎng)絡(luò)建模方法，其特征在于：包基于目標(biāo)區(qū)域的地震測井資料參數(shù)，構(gòu)建初始化斷層面模型，并采用加密網(wǎng)格化方法沿所述初始化斷層面模型的走向方向提取基本斷層點集數(shù)據(jù)；提取目標(biāo)區(qū)域的地層參數(shù)，結(jié)合平衡剖面恢復(fù)方法，確定目標(biāo)區(qū)域在歷史時期的變形過程；根據(jù)所述地層參數(shù)標(biāo)定離散元的細(xì)觀參數(shù)，沉積第一離散元模型，構(gòu)建斷層函數(shù)，并優(yōu)化生成第二離散元模型；基于所述第二離散元模型，設(shè)置邊界條件進行逆沖推覆模擬，并通過所述第二離散元模型的輸出結(jié)果與實際地質(zhì)剖面的對比，調(diào)整模型參數(shù)，直至二者匹配；基于匹配后的模型結(jié)果，建立并分析目標(biāo)區(qū)域的三維斷層和裂縫網(wǎng)絡(luò)的空間展布特征，并對三維斷層的活動性進行動態(tài)評估與趨勢預(yù)測；基于第二離散元模型演化結(jié)果的顆粒坐標(biāo)，根據(jù)顆粒間相對位置向量，推導(dǎo)出不同顆粒位移梯度張量；根據(jù)第二離散元模型演化結(jié)果的顆粒位移梯度張量，繪制不同剖面的應(yīng)變張量，得到空間展布的三維斷層；基于第二離散元模型演化結(jié)果中各顆粒間的粘結(jié)關(guān)系，判斷顆粒間粘結(jié)是否發(fā)生破根據(jù)判斷結(jié)果記錄顆粒間粘結(jié)破壞的時間、位置和類型，其中所述類型包括張破裂和剪破裂；將顆粒間粘結(jié)狀態(tài)標(biāo)記為裂縫產(chǎn)生點，其中所述裂縫產(chǎn)生點為顆粒間粘結(jié)狀態(tài)由粘結(jié)壓縮狀態(tài)或粘結(jié)拉伸狀態(tài)轉(zhuǎn)變?yōu)闊o粘結(jié)接觸狀態(tài)或失去相互作用狀態(tài)的情況；對張破裂類型和剪破裂類型進行特征分析，其中所述張破裂類型表現(xiàn)為垂直于最大主應(yīng)力方向上的破壞；所述剪破裂類型表現(xiàn)為沿著最小主應(yīng)力方向發(fā)展的破壞；基于加載過程的推進，持續(xù)更新顆粒間粘結(jié)破壞的信息，建立裂縫網(wǎng)絡(luò)數(shù)據(jù)庫，其中所述裂縫網(wǎng)絡(luò)數(shù)據(jù)庫包括已發(fā)生粘結(jié)破壞的顆粒坐標(biāo)、破壞時間和破壞類型；基于所述裂縫網(wǎng)絡(luò)數(shù)據(jù)庫，進行裂縫網(wǎng)絡(luò)定量分析，計算統(tǒng)計指標(biāo)，其中所述統(tǒng)計指標(biāo)包括裂縫密度、平均裂縫長度和裂縫方向分布；通過所述統(tǒng)計指標(biāo)繪制不同時空演化下裂縫網(wǎng)絡(luò)的空間展布特征圖。2.如權(quán)利要求1所述的基于離散元的逆沖推覆體系下三維斷層/裂縫網(wǎng)絡(luò)建模方法，其特征在于：基于所述第二離散元模型，設(shè)置邊界條件進行逆沖推覆模擬，并通過所述第二離將第二離散元模型的后緣設(shè)置為活動邊界，其余邊界設(shè)置為固定邊界；參考目標(biāo)區(qū)域的斷層位移量和平衡剖面恢復(fù)縮短量，對所述活動邊界施加運動；控制所述活動邊界的運動，當(dāng)達到與平衡剖面恢復(fù)縮短量相等的擠壓量時，則停止所述第二離散元模型的運算；提取所述第二離散元模型的模擬結(jié)果，將所述模擬結(jié)果與實際地質(zhì)剖面進行一致性對若所述模擬結(jié)果與實際地質(zhì)剖面不匹配，則微調(diào)離散元顆粒參數(shù)或邊界運動情況，重3復(fù)上述步驟，直至所述模擬結(jié)果與實際地質(zhì)剖面匹配；若所述模擬結(jié)果與所述實際地質(zhì)剖面匹配，則執(zhí)行三維斷層/裂縫網(wǎng)絡(luò)構(gòu)建。3.如權(quán)利要求2所述的基于離散元的逆沖推覆體系下三維斷層/裂縫網(wǎng)絡(luò)建模方法，其特征在于：根據(jù)所述地層參數(shù)標(biāo)定離散元的細(xì)觀參數(shù)，沉積第一離散元模型，構(gòu)建斷層函建立離散元的細(xì)觀參數(shù)體系，其中所述細(xì)觀參數(shù)體系包括顆粒參數(shù)和顆粒粘結(jié)參數(shù)；確定離散元相鄰顆粒間的狀態(tài)關(guān)系，并通過三軸壓縮實驗測試校準(zhǔn)各個地層的顆粒材料細(xì)觀參數(shù)，其中所述狀態(tài)關(guān)系包括粘結(jié)壓縮狀態(tài)、粘結(jié)拉伸狀態(tài)、無粘結(jié)接觸狀態(tài)以及顆粒失去相互作用狀態(tài)；基于目標(biāo)區(qū)域相應(yīng)的空間，隨機分布半徑呈正態(tài)分布的顆粒，并建立顆粒細(xì)觀參數(shù)與地層宏觀性質(zhì)的對應(yīng)關(guān)系，賦予每個顆粒物理屬性和相互作用參數(shù)，其中所述物理屬性包括半徑和密度；所述相互作用參數(shù)包括法向剛度、剪切剛度和摩擦系數(shù)；將賦予后的顆粒自由沉降至系統(tǒng)穩(wěn)定，記錄各個地層顆粒的位置參數(shù)，完成第一離散元模型的構(gòu)建；基于標(biāo)定的顆粒粘結(jié)參數(shù)和目標(biāo)區(qū)域的地層信息，在所述第一離散元模型的顆粒之間設(shè)定粘結(jié)關(guān)系；對所述第一離散元模型進行細(xì)化，建立目標(biāo)區(qū)域不同層位顆粒的細(xì)觀參數(shù)和力學(xué)關(guān)根據(jù)斷層點集數(shù)據(jù)Fault,建立主要斷層模型函數(shù)isInsideFault,其中所述主要斷層模型函數(shù)用于判斷任一顆粒坐標(biāo)pt是否位于所述斷層點集Fault內(nèi)部；基于判斷結(jié)果，取消斷層面上顆粒間的粘結(jié)，同時降低斷層顆粒間摩擦系數(shù)，形成第二離散元模型。4.如權(quán)利要求3所述的基于離散元的逆沖推覆體系下三維斷層/裂縫網(wǎng)絡(luò)建模方法，其特征在于：提取目標(biāo)區(qū)域的地層參數(shù)，結(jié)合平衡剖面恢復(fù)方法，確定目標(biāo)區(qū)域在歷史時期的基于初始化斷層面模型，收集目標(biāo)區(qū)域的鉆井資料，從所述鉆井資料中獲取巖芯樣本；對所述巖芯樣本進行實驗室分析，測定不同地層的物理性質(zhì)參數(shù)，其中所述物理性質(zhì)采集并分析目標(biāo)區(qū)域的測井?dāng)?shù)據(jù)，獲取地下地層的詳細(xì)信息，其中所述測井?dāng)?shù)據(jù)包括電阻率測井?dāng)?shù)據(jù)和自然伽馬射線測井?dāng)?shù)據(jù)；解析地震剖面數(shù)據(jù)，提取目標(biāo)區(qū)域不同地層的信息，其中所述不同地層的信息包括地層巖性、力學(xué)參數(shù)和空間展布特征；所述力學(xué)參數(shù)包括彈性模量和泊松比；將所述不同地層的信息、所述物理性質(zhì)參數(shù)和所述詳細(xì)信息進行綜合分析，得到目標(biāo)區(qū)域的完整地層參數(shù)；基于所述完整地層參數(shù)，采用平衡剖面恢復(fù)方法對所述目標(biāo)區(qū)域進行建模；通過所述平衡剖面恢復(fù)方法確定目標(biāo)區(qū)域在歷史時期的變形過程，識別斷層的多期活對所述平衡剖面恢復(fù)方法生成的恢復(fù)后平衡剖面與原始剖面進行量化對比分析，計算每期斷層活動引起的擠壓變形位移。45.如權(quán)利要求4所述的基于離散元的逆沖推覆體系下三維斷層/裂縫網(wǎng)絡(luò)建模方法，其特征在于：所述初始化斷層面模型的構(gòu)建方法為，獲取目標(biāo)區(qū)域的地震測井資料參數(shù)，并對所述地震測井資料參數(shù)進行預(yù)處理；基于所述預(yù)處理后的地震測井資料參數(shù)，挑選出清晰顯示目標(biāo)斷層特征的地震剖面，優(yōu)先關(guān)注顯示出位移的區(qū)域；對所述地震剖面進行分析，識別并提取位移大且方向為擠壓方向的主要斷層的相關(guān)參數(shù)，其中所述相關(guān)參數(shù)包括斷層面的傾向、斷層面的傾角、沿斷層面的滑動距離、沿斷層面的滑動方向以及斷層兩側(cè)地層的相對位移量；基于所述相關(guān)參數(shù)，提取主要斷層的空間三維幾何形狀；對所述空間三維幾何形狀采用加密網(wǎng)格化方法進行處理；基于處理后的空間三維幾何形狀，沿主要斷層的走向方向提取基本斷層點集數(shù)據(jù)，設(shè)所述基本斷層點集數(shù)據(jù)為Fault,包括n個三維點S={p?,p?,…,p},其中，p,為斷層面上的第n個三維點；根據(jù)所述基本斷層點集數(shù)據(jù)構(gòu)建初始化斷層面模型。6.一種基于離散元的逆沖推覆體系下三維斷層/裂縫網(wǎng)絡(luò)建模裝置，基于權(quán)利要求1~5任一所述的基于離散元的逆沖推覆體系下三維斷層/裂縫網(wǎng)絡(luò)建模方法，其特征在于：包斷層模塊，基于目標(biāo)區(qū)域的地震測井資料參數(shù)，構(gòu)建初始化斷層面模型，并采用加密網(wǎng)格化方法沿所述初始化斷層面模型的走向方向提取基本斷層點集數(shù)據(jù)；反演模塊，用于通過觀測目標(biāo)區(qū)域的地質(zhì)并解析該地區(qū)地震資料，提取地層信息與變形過程，并對地震解析圖進行反演，斷層的空模型模塊，根據(jù)所述地層參數(shù)標(biāo)定離散元細(xì)觀參數(shù)，沉積第一離散元模型，構(gòu)建斷層函變形模塊，基于所述第二離散元模型，設(shè)置邊界條件進行逆沖推覆模擬，并通過所述第二離散元模型的輸出結(jié)果與實際地質(zhì)剖面的對比，調(diào)整模型參數(shù)，直至二者匹配；建模模塊，基于匹配后的模型結(jié)果，建立并分析目標(biāo)區(qū)域的三維斷層和裂縫網(wǎng)絡(luò)的空間展布特征，并對三維斷層的活動性進行動態(tài)評估與趨勢預(yù)測。7.一種計算機設(shè)備，包括存儲器和處理器，所述存儲器存儲有計算機程序，其特征在于：所述處理器執(zhí)行所述計算機程序時實現(xiàn)權(quán)利要求1~5任一所述的基于離散元的逆沖推覆體系下三維斷層/裂縫網(wǎng)絡(luò)建模方法的步驟。8.一種計算機可讀存儲介質(zhì)，其上存儲有計算機程序，其特征在于：所述計算機程序被處理器執(zhí)行時實現(xiàn)權(quán)利要求1~5任一所述的基于離散元的逆沖推覆體系下三維斷層/裂縫網(wǎng)絡(luò)建模方法的步驟。5一種基于離散元的逆沖推覆體系下三維斷層/裂縫網(wǎng)絡(luò)建模方法及裝置技術(shù)領(lǐng)域[0001]本發(fā)明涉及三維斷層建模研究技術(shù)領(lǐng)域，特別是一種基于離散元的逆沖推覆體系下三維斷層/裂縫網(wǎng)絡(luò)建模方法及裝置。背景技術(shù)[0002]斷層與裂縫作為地質(zhì)構(gòu)造活動的重要產(chǎn)物，是控制油氣運移、聚集及成藏的關(guān)鍵地質(zhì)要素。斷層指巖層受構(gòu)造應(yīng)力作用發(fā)生斷裂并產(chǎn)生顯著位移的結(jié)構(gòu)，而裂縫則表征巖石破裂形成的非連續(xù)界面。二者共同構(gòu)成的復(fù)雜網(wǎng)絡(luò)系統(tǒng)不僅直接影響油氣儲集空間的分布規(guī)律，更通過動態(tài)演化過程對圈閉完整性、流體輸導(dǎo)能力產(chǎn)生持續(xù)性改造作用。因此，構(gòu)建三維斷層/裂縫網(wǎng)絡(luò)模型已成為油氣勘探開發(fā)中破解構(gòu)造控藏機制的核心技術(shù)環(huán)節(jié)。[0003]當(dāng)前三維地質(zhì)建模面臨雙重挑戰(zhàn)：一方面，地表與深部斷層/裂縫系統(tǒng)具有顯著的空間非均質(zhì)性與幾何不連續(xù)性，其原始三維觀測數(shù)據(jù)往往呈現(xiàn)碎片化特征，復(fù)雜斷層構(gòu)造相互切割導(dǎo)致地層的連續(xù)性遭到嚴(yán)重破壞，各類空間插值算法及面模型、體模型的斷層構(gòu)建方法不再適用；另一方面，傳統(tǒng)建模方法多局限于靜態(tài)層狀結(jié)構(gòu)重建，在構(gòu)建復(fù)雜且層面較多的斷層模型時會極大增加建模難度和復(fù)雜程度，難以有效刻畫多期次構(gòu)造運動疊加形成的逆沖推覆體系。尤其當(dāng)涉及斷層尖端應(yīng)力奇異性、多斷層交互擴展及裂縫網(wǎng)絡(luò)動態(tài)重組等復(fù)雜過程時，常規(guī)建模技術(shù)常因網(wǎng)格重構(gòu)困難、演化過程缺失等問題導(dǎo)致預(yù)測結(jié)果存在顯著偏差。這種局限性在具有強構(gòu)造改造特征的含油氣盆地中尤為突出，使得二維構(gòu)造解釋的多解性問題無法通過簡單三維可視化得到根本解決。[0004]離散元模擬(DiscreteElementMethod)其基本思想是將材料內(nèi)部細(xì)觀尺度的單個離散顆粒視為一個離散單元，將顆粒集合體視為若干離散單元的集合，通過一系列離散單元來模擬材料科學(xué)的力學(xué)行為，在研究非連續(xù)介質(zhì)破裂問題、大變形問題時，有很大的優(yōu)勢。基于離散元方法的逆沖推覆體系三維斷層/裂縫建模技術(shù)改進傳統(tǒng)建模范式，該方法將地質(zhì)歷史中的應(yīng)力場演變、斷層成核擴展規(guī)律等動力學(xué)要素融入建模過程，采用離散元數(shù)值模擬再現(xiàn)構(gòu)造變形序列。該研究方法及裝置構(gòu)建三維地質(zhì)建模正從靜態(tài)描述向動態(tài)預(yù)測跨越，為復(fù)雜構(gòu)造區(qū)油氣勘探提供了全新的技術(shù)方法。發(fā)明內(nèi)容[0005]鑒于現(xiàn)有的基于地震資料的三維斷層裂縫建模方法在面對推覆構(gòu)造帶等復(fù)雜變及傳統(tǒng)方法無法動態(tài)反演斷裂演化過程，導(dǎo)致模型的空間精度與地質(zhì)合理性難以兼顧，預(yù)測結(jié)果偏差較大。此外，常規(guī)網(wǎng)格重構(gòu)在處理斷層尖端破裂、斷層交匯擴展及裂縫動態(tài)演化等過程中存在難以描述的應(yīng)力奇異性和建模障礙的問題，提出了本發(fā)明。[0006]因此，本發(fā)明所要解決的問題在于如何構(gòu)建一種結(jié)合地質(zhì)力學(xué)機制、能夠動態(tài)模擬斷層裂縫形成演化過程、適用于逆沖推覆構(gòu)造帶的高精度三維建模方法，以實現(xiàn)復(fù)雜斷6裂系統(tǒng)的時空展布特征精準(zhǔn)刻畫、斷層活動性動態(tài)評估，并克服傳統(tǒng)方法在建模精度、演化過程重建以及多期構(gòu)造疊加表達上的不足。[0007]為解決上述技術(shù)問題，本發(fā)明提供如下技術(shù)方案：[0008]第一方面，本發(fā)明提供了一種基于離散元的逆沖推覆體系下三維斷層/裂縫網(wǎng)絡(luò)[0009]基于目標(biāo)區(qū)域的地震測井資料參數(shù)，構(gòu)建初始化斷層面模型，并采用加密網(wǎng)格化方法沿所述初始化斷層面模型的走向方向提取基本斷層點集數(shù)據(jù)；[0010]提取目標(biāo)區(qū)域的地層參數(shù)，結(jié)合平衡剖面恢復(fù)方法，確定目標(biāo)區(qū)域在歷史時期的變形過程；[0011]根據(jù)所述地層參數(shù)標(biāo)定離散元的細(xì)觀參數(shù)，沉積第一離散元模型，構(gòu)建斷層函數(shù)，并優(yōu)化生成第二離散元模型；[0012]基于所述第二離散元模型，設(shè)置邊界條件進行逆沖推覆模擬，并通過所述第二離散元模型的輸出結(jié)果與實際地質(zhì)剖面的對比，調(diào)整模型參數(shù)，直至二者匹配；[0013]基于模擬結(jié)果，建立三維斷層和裂縫網(wǎng)絡(luò)，分析和建立目標(biāo)區(qū)域三維斷層和裂縫網(wǎng)絡(luò)的空間展布特征，并對三維斷層的活動性進行動態(tài)評估與趨勢預(yù)測。[0014]作為本發(fā)明所述基于離散元的逆沖推覆體系下三維斷層/裂縫網(wǎng)絡(luò)建模方法的一種優(yōu)選方案，其中：基于匹配后的模型結(jié)果，建立并分析目標(biāo)區(qū)域的三維斷層和裂縫網(wǎng)絡(luò)的空間展布特征，并對三維斷層的活動性進行動態(tài)評估與趨勢預(yù)測，包括：[0015]基于第二離散元模型演化結(jié)果的顆粒坐標(biāo)，根據(jù)顆粒間相對位置向量，推導(dǎo)出不同顆粒位移梯度張量；[0016]根據(jù)第二離散元模型演化結(jié)果的顆粒位移梯度張量，繪制不同剖面的應(yīng)變張量，得到空間展布的三維斷層；[0017]基于第二離散元模型演化結(jié)果中各顆粒間的粘結(jié)關(guān)系，判斷顆粒間粘結(jié)是否發(fā)生破壞；[0018]根據(jù)判斷結(jié)果記錄顆粒間粘結(jié)破壞的時間、位置和類型，其中所述類型包括張破裂和剪破裂；[0019]將顆粒間粘結(jié)狀態(tài)標(biāo)記為裂縫產(chǎn)生點，其中所述裂縫產(chǎn)生點為顆粒間粘結(jié)狀態(tài)由粘結(jié)壓縮狀態(tài)或粘結(jié)拉伸狀態(tài)轉(zhuǎn)變?yōu)闊o粘結(jié)接觸狀態(tài)或失去相互作用狀態(tài)的情況；[0020]對張破裂類型和剪破裂類型進行特征分析，其中所述張破裂類型表現(xiàn)為垂直于最大主應(yīng)力方向上的破壞；所述剪破裂類型表現(xiàn)為沿著最小主應(yīng)力方向發(fā)展的破壞；[0021]基于加載過程的推進，持續(xù)更新顆粒間粘結(jié)破壞的信息，建立裂縫網(wǎng)絡(luò)數(shù)據(jù)庫，其中所述裂縫網(wǎng)絡(luò)數(shù)據(jù)庫包括已發(fā)生粘結(jié)破壞的顆粒坐標(biāo)、破壞時間和破壞類型；[0022]基于所述裂縫網(wǎng)絡(luò)數(shù)據(jù)庫，進行裂縫網(wǎng)絡(luò)定量分析，計算統(tǒng)計指標(biāo)，其中所述統(tǒng)計指標(biāo)包括裂縫密度、平均裂縫長度和裂縫方向分布；[0023]通過所述統(tǒng)計指標(biāo)繪制不同時空演化下裂縫網(wǎng)絡(luò)的空間展布特征圖。[0024]作為本發(fā)明所述基于離散元的逆沖推覆體系下三維斷層/裂縫網(wǎng)絡(luò)建模方法的一種優(yōu)選方案，其中：基于所述第二離散元模型，設(shè)置邊界條件進行逆沖推覆模擬，并通過所述第二離散元模型的輸出結(jié)果與實際地質(zhì)剖面的對比，調(diào)整模型參數(shù)，直至二者匹配，包7[0025]將第二離散元模型的后緣設(shè)置為活動邊界，其余邊界設(shè)置為固定邊界；[0026]參考目標(biāo)區(qū)域的斷層位移量和平衡剖面恢復(fù)縮短量，對所述活動邊界施加運動；[0027]控制所述活動邊界的運動，當(dāng)達到與平衡剖面恢復(fù)縮短量相等的擠壓量時，則停止所述第二離散元模型的運算；[0028]提取所述第二離散元模型的模擬結(jié)果，將所述模擬結(jié)果與實際地質(zhì)剖面進行一致性對比；[0029]若所述模擬結(jié)果與實際地質(zhì)剖面不匹配，則微調(diào)離散元顆粒參數(shù)或邊界運動情況，重復(fù)上述步驟，直至所述模擬結(jié)果與實際地[0030]若所述模擬結(jié)果與所述實際地質(zhì)剖面匹配，則執(zhí)行三維斷層/裂縫網(wǎng)絡(luò)構(gòu)建。[0031]作為本發(fā)明所述基于離散元的逆沖推覆體系下三維斷層/裂縫網(wǎng)絡(luò)建模方法的一種優(yōu)選方案，其中：根據(jù)所述地層參數(shù)標(biāo)定離散元的細(xì)觀參數(shù)，沉積第一離散元模型，構(gòu)建[0032]建立離散元的細(xì)觀參數(shù)體系，其中所述細(xì)觀參數(shù)體系包括顆粒參數(shù)和顆粒粘結(jié)參[0033]確定離散元相鄰顆粒間的狀態(tài)關(guān)系，并通過三軸壓縮實驗測試校準(zhǔn)各個地層的顆粒材料細(xì)觀參數(shù)，其中所述狀態(tài)關(guān)系包括粘結(jié)壓縮狀態(tài)、粘結(jié)拉伸狀態(tài)、無粘結(jié)接觸狀態(tài)以及顆粒失去相互作用狀態(tài)；[0034]基于目標(biāo)區(qū)域相應(yīng)的空間，隨機分布半徑呈正態(tài)分布的顆粒，并建立顆粒細(xì)觀參數(shù)與地層宏觀性質(zhì)的對應(yīng)關(guān)系，賦予每個顆粒物理屬性和相互作用參數(shù)，其中所述物理屬性包括半徑和密度；所述相互作用參數(shù)包括法向剛度、剪切剛度和摩擦系數(shù)；[0035]將賦予后的顆粒自由沉降至系統(tǒng)穩(wěn)定，記錄各個地層顆粒的位置參數(shù)，完成第一離散元模型的構(gòu)建；[0036]基于標(biāo)定的顆粒粘結(jié)參數(shù)和目標(biāo)區(qū)域的地層信息，在所述第一離散元模型的顆粒之間設(shè)定粘結(jié)關(guān)系；[0037]對所述第一離散元模型進行細(xì)化，建立目標(biāo)區(qū)域不同層位顆粒的細(xì)觀參數(shù)和力學(xué)[0038]根據(jù)斷層點集數(shù)據(jù)Fault,建立主要斷層模型函數(shù)isInsideFault,其中所述主要斷層模型函數(shù)用于判斷任一顆粒坐標(biāo)pt是否位于所述斷層點集Fault內(nèi)部；[0039]基于判斷結(jié)果，取消斷層面上顆粒間的粘結(jié)，同時降低斷層顆粒間摩擦系數(shù)，形成第二離散元模型。[0040]作為本發(fā)明所述基于離散元的逆沖推覆體系下三維斷層/裂縫網(wǎng)絡(luò)建模方法的一種優(yōu)選方案，其中：提取目標(biāo)區(qū)域的地層參數(shù)，結(jié)合平衡剖面恢復(fù)方法，確定目標(biāo)區(qū)域在歷[0041]基于初始化斷層面模型，收集目標(biāo)區(qū)域的鉆井資料，從所述鉆井資料中獲取巖芯樣本；[0042]對所述巖芯樣本進行實驗室分析，測定不同地層的物理性質(zhì)參數(shù)，其中所述物理性質(zhì)參數(shù)包括巖性、孔隙度和滲透率；[0043]采集并分析目標(biāo)區(qū)域的測井?dāng)?shù)據(jù)，獲取地下地層的詳細(xì)信息，其中所述測井?dāng)?shù)據(jù)包括電阻率測井?dāng)?shù)據(jù)和自然伽馬射線測井?dāng)?shù)據(jù)；8[0044]解析地震剖面數(shù)據(jù)，提取目標(biāo)區(qū)域不同地層的信息，其中所述不同地層的信息包括地層巖性、力學(xué)參數(shù)和空間展布特征；所述力學(xué)參數(shù)包括彈性模量和泊松比；[0045]將所述不同地層的信息、所述物理性質(zhì)參數(shù)和所述詳細(xì)信息進行綜合分析，得到目標(biāo)區(qū)域的完整地層參數(shù)；[0046]基于所述完整地層參數(shù)，采用平衡剖面恢復(fù)方法對所述目標(biāo)區(qū)域進行建模；[0047]通過所述平衡剖面恢復(fù)方法確定目標(biāo)區(qū)域在歷史時期的變形過程，識別斷層的多期活動；[0048]對所述平衡剖面恢復(fù)方法生成的恢復(fù)后平衡剖面與原始剖面進行量化對比分析，計算每期斷層活動引起的擠壓變形位移。[0049]作為本發(fā)明所述基于離散元的逆沖推覆體系下三維斷層/裂縫網(wǎng)絡(luò)建模方法的一[0050]獲取目標(biāo)區(qū)域的地震測井資料參數(shù)，并對所述地震測井資料參數(shù)進行預(yù)處理；[0051]基于所述預(yù)處理后的地震測井資料參數(shù)，挑選出清晰顯示目標(biāo)斷層特征的地震剖[0053]對所述地震剖面進行分析，識別并提取位移大且方向為擠壓方向的主要斷層的相關(guān)參數(shù)，其中所述相關(guān)參數(shù)包括斷層面的傾向、斷層面的傾角、沿斷層面的滑動距離、沿斷層面的滑動方向以及斷層兩側(cè)地層的相對位移量；[0054]基于所述相關(guān)參數(shù)，提取主要斷層的空間三維幾何形狀；[0055]對所述空間三維幾何形狀采用加密網(wǎng)格化方法進行處理；[0056]基于處理后的空間三維幾何形狀，沿主要斷層的走向方向提取基本斷層點集數(shù)據(jù)，設(shè)所述基本斷層點集數(shù)據(jù)為Fault,包括n個三維點S={p?,p?,…,p?},其中，p。為斷層面上的第n個三維點；[0057]根據(jù)所述基本斷層點集數(shù)據(jù)構(gòu)建初始化斷層面模型。[0058]第二方面，本發(fā)明提供了一種基于離散元的逆沖推覆體系下三維斷層/裂縫網(wǎng)絡(luò)[0059]斷層模塊，基于目標(biāo)區(qū)域的地震測井資料參數(shù)，構(gòu)建初始化斷層面模型，并采用加密網(wǎng)格化方法沿所述初始化斷層面模型的走向方向提取基本斷層點集數(shù)據(jù)；[0060]反演模塊，用于通過觀測目標(biāo)區(qū)域的地質(zhì)并解析該地區(qū)地震資料，提取地層信息[0061]模型模塊，根據(jù)所述地層參數(shù)標(biāo)定離散元細(xì)觀參數(shù)，沉積第一離散元模型，構(gòu)建斷[0062]變形模塊，基于所述第二離散元模型，設(shè)置邊界條件進行逆沖推覆模擬，并通過所述第二離散元模型的輸出結(jié)果與實際地質(zhì)剖面的對比，調(diào)整模型參數(shù)，直至二者匹配；[0063]建模模塊，基于匹配后的模型結(jié)果，建立并分析目標(biāo)區(qū)域的三維斷層和裂縫網(wǎng)絡(luò)的空間展布特征，并對三維斷層的活動性進行動態(tài)評估與趨勢預(yù)測。[0064]第三方面，本發(fā)明提供了一種計算機設(shè)備，包括存儲器和處理器，所述存儲器存儲有計算機程序，其中：所述處理器執(zhí)行所述計算機程序時實現(xiàn)上述的基于離散元的逆沖推覆體系下三維斷層/裂縫網(wǎng)絡(luò)建模方法的任一步驟。9[0065]第四方面，本發(fā)明提供了一種計算機可讀存儲介質(zhì)，其上存儲有計算機程序，其中：所述計算機程序被處理器執(zhí)行時實現(xiàn)上述的基于離散元的逆沖推覆體系下三維斷層/裂縫網(wǎng)絡(luò)建模方法的任一步驟。[0066]與現(xiàn)有技術(shù)對比，本發(fā)明有益效果為：通過引入離散元方法對斷層裂縫系統(tǒng)進行動態(tài)模擬，有效規(guī)避了傳統(tǒng)建模中因網(wǎng)格重構(gòu)困難導(dǎo)致的斷層尖端破裂擴展行為難以精準(zhǔn)刻畫的問題；通過在建模過程中融合地質(zhì)力學(xué)參數(shù)與構(gòu)造演化過程，建立了先存斷層與新生裂縫之間的時空耦合關(guān)系，實現(xiàn)了多期構(gòu)造疊加效應(yīng)的量化解析；同時建立了從構(gòu)造跡象反演到未來構(gòu)造演化預(yù)測的閉環(huán)建模流程，顯著提升了復(fù)雜斷裂系統(tǒng)三維空間展布建模的地質(zhì)合理性與預(yù)測能力，為油氣勘探中隱伏構(gòu)造識別、斷層封堵性評價及構(gòu)造控藏機制研究提供了科學(xué)、可靠的技術(shù)支撐。附圖說明[0067]為了更清楚地說明本發(fā)明實施例的技術(shù)方案，下面將對實施例描述中所需要使用的附圖作簡單地介紹，顯而易見地，下面描述中的附圖僅僅是本發(fā)明的一些實施例，對于本領(lǐng)域普通技術(shù)人員來講，在不付出創(chuàng)造性勞動的前提下，還可以根據(jù)這些附圖獲得其它的[0068]圖1為基于離散元的逆沖推覆體系下三維斷層/裂縫網(wǎng)絡(luò)建模方法的流程示意圖。[0069]圖2為基于離散元的逆沖推覆體系下三維斷層/裂縫網(wǎng)絡(luò)建模方法的沉積模型。[0070]圖3為基于離散元的逆沖推覆體系下三維斷層/裂縫網(wǎng)絡(luò)建模方法的顆粒接觸和粘結(jié)關(guān)系。[0071]圖4為基于離散元的逆沖推覆體系下三維斷層/裂縫網(wǎng)絡(luò)建模方法的示例平衡剖面演化過程。[0072]圖5為基于離散元的逆沖推覆體系下三維斷層/裂縫網(wǎng)絡(luò)建模方法的模擬裝置結(jié)構(gòu)示意圖。[0073]圖6為基于離散元的逆沖推覆體系下三維斷層/裂縫網(wǎng)絡(luò)建模方法的所示的示例區(qū)域中離散元地質(zhì)模型與實際地質(zhì)情況對比。[0074]圖7為基于離散元的逆沖推覆體系下三維斷層/裂縫網(wǎng)絡(luò)建模方法的模擬地層的應(yīng)力應(yīng)變曲線。[0075]圖8為基于離散元的逆沖推覆體系下三維斷層/裂縫網(wǎng)絡(luò)建模方法的示例瑪扎塔格構(gòu)造帶地震剖面。[0076]圖9為基于離散元的逆沖推覆體系下三維斷層/裂縫網(wǎng)絡(luò)建模方法的示例斷層組[0077]圖10為基于離散元的逆沖推覆體系下三維斷層/裂縫網(wǎng)絡(luò)建模方法的示例斷層縱[0078]圖11為基于離散元的逆沖推覆體系下三維斷層/裂縫網(wǎng)絡(luò)建模方法的示例斷層橫[0079]圖12為基于離散元的逆沖推覆體系下三維斷層/裂縫網(wǎng)絡(luò)建模方法的示例裂縫密度統(tǒng)計。具體實施方式[0080]為使本發(fā)明的上述目的、特征和優(yōu)點能夠更加明顯易懂，下面結(jié)合說明書附圖對本發(fā)明的具體實施方式做詳細(xì)的說明，顯然所描述的實施例是本發(fā)明的一部分實施例，而不是全部實施例?；诒景l(fā)明中的實施例，本領(lǐng)域普通人員在沒有做出創(chuàng)造性勞動前提下所獲得的所有其他實施例，都應(yīng)當(dāng)屬于本發(fā)明的保護的范圍。[0081]在下面的描述中闡述了很多具體細(xì)節(jié)以便于充分理解本發(fā)明，但是本發(fā)明還可以采用其他不同于在此描述的其它方式來實施，本領(lǐng)域技術(shù)人員可以在不違背本發(fā)明內(nèi)涵的情況下做類似推廣，因此本發(fā)明不受下面公開的具體實施例的限制。式中的特定特征、結(jié)構(gòu)或特性。在本說明書中不同地方出現(xiàn)的“在一個實施例中”并非均指同一個實施例，也不是單獨的或選擇性的與其他實施例互相排斥的實施例。[0083]本發(fā)明結(jié)合示意圖進行詳細(xì)描述，在詳述本發(fā)明實施例時，為便于說明，表示器件結(jié)構(gòu)的剖面圖會不依一般比例作局部放大，而且所述示意圖只是示例，其在此不應(yīng)限制本或位置關(guān)系為基于附圖所示的方位或位置關(guān)系，僅是為了便于描述本發(fā)明和簡化描述，而不是指示或暗示所指的裝置或元件必須具有特定的方位、以特定的方位構(gòu)造和操作，因此為指示或暗示相對重要性。也可以通過中間媒介間接相連，也可以是兩個元件內(nèi)部的連通。對于本領(lǐng)域的普通技術(shù)人員而言，可以具體情況理解上述術(shù)語在本發(fā)明中的具體含義。[0087]參照圖1~圖5,為本發(fā)明第一個實施例，該實施例提供了一種基于離散元的逆沖推覆體系下三維斷層/裂縫網(wǎng)絡(luò)建模方法，包括：[0088]S1:基于目標(biāo)區(qū)域的地震測井資料參數(shù)，構(gòu)建初始化斷層面模型，并采用加密網(wǎng)格化方法沿初始化斷層面模型的走向方向提取基本斷層點集數(shù)據(jù)。[0089]具體的，初始化斷層面模型的構(gòu)建方法為，獲取目標(biāo)區(qū)域的地震測井資料參數(shù)，并對地震測井資料參數(shù)進行預(yù)處理。[0090]在可選的實施方式中，通過地質(zhì)調(diào)查和地震資料解析獲取目標(biāo)區(qū)域的地震測井資料參數(shù)，并采用Petrol、LandMark軟件對地震測井資料參數(shù)進行預(yù)處理后，包括噪音去除和反射波提取。[0091]進一步的，基于預(yù)處理后的地震測井資料參數(shù)，挑選出清晰顯示目標(biāo)斷層特征的地震剖面，優(yōu)先關(guān)注顯示出位移的區(qū)域；基于地震剖面，提取坐標(biāo)數(shù)據(jù)、分層數(shù)據(jù)、斷層數(shù)據(jù)；對地震剖面進行分析，識別并提取位移大且方向為擠壓方向的主要斷層的相關(guān)參數(shù)，其中相關(guān)參數(shù)包括斷層面的傾向、斷層面的傾角、沿斷層面的滑動距離、沿斷層面的滑動方向以及斷層兩側(cè)地層的相對位移量；基于相關(guān)參數(shù)，提取主要斷層的空間三維幾何形狀；對空間三維幾何形狀采用加密網(wǎng)格化方法進行處理。11[0092]優(yōu)選的，基于處理后的空間三維幾何形狀，沿主要斷層的走向方向提取基本斷層點集數(shù)據(jù)，設(shè)基本斷層點集數(shù)據(jù)為Fault,包括n個三維點S={p?,p?,…,p?},其中，p為斷層面上的第n個三維點；[0093]更進一步的，根據(jù)基本斷層點集數(shù)據(jù)構(gòu)建初始化斷層面模型。[0094]S2:提取目標(biāo)區(qū)域的地層參數(shù)，結(jié)合平衡剖面恢復(fù)方法，確定目標(biāo)區(qū)域在歷史時期的變形過程。[0095]具體的，收集目標(biāo)區(qū)域的鉆井資料，從鉆井資料中獲取巖芯樣本；對巖芯樣本進行實驗室分析，測定不同地層的物理性質(zhì)參數(shù)，其中物理性質(zhì)參數(shù)包括巖性、孔隙度和滲透率；采集并分析目標(biāo)區(qū)域的測井?dāng)?shù)據(jù)，獲取地下地層的詳細(xì)信息，其中測井?dāng)?shù)據(jù)包括電阻率測井?dāng)?shù)據(jù)和自然伽馬射線測井?dāng)?shù)據(jù)。[0096]需要說明的是，地下地層的詳細(xì)信息包括地層厚度、巖性分界和流體飽[0097]進一步的，解析地震剖面數(shù)據(jù)，提取目標(biāo)區(qū)域不同地層的信息，不同地層的信息包括地層巖性、力學(xué)參數(shù)和空間展布特征；力學(xué)參數(shù)包括彈性模量和泊松比；將不同地層的信息、物理性質(zhì)參數(shù)和詳細(xì)信息進行綜合分析，得到目標(biāo)區(qū)域的完整地層參數(shù)；基于完整地層參數(shù)，采用平衡剖面恢復(fù)方法對目標(biāo)區(qū)域進行建模。[0098]需要說明的是，平衡剖面恢復(fù)方法是通過平衡剖面恢復(fù)軟件實現(xiàn)的，其中平衡剖面恢復(fù)軟件為3dmove。[0099]更進一步的，如圖4所示，通過平衡剖面恢復(fù)方法確定目標(biāo)區(qū)域在歷史時期的變形過程，識別斷層的多期活動；對平衡剖面恢復(fù)方法生成的恢復(fù)后平衡剖面與原始剖面進行量化對比分析，計算每期斷層活動引起的擠壓變形位移。[0100]S3:根據(jù)地層參數(shù)標(biāo)定離散元的細(xì)觀參數(shù)，沉積第一離散元模型，構(gòu)建斷層函數(shù)，并優(yōu)化生成第二離散元模型。[0101]具體的，如圖2所示，建立離散元細(xì)觀參數(shù)體系，其中離散元細(xì)觀參數(shù)體系包括顆粒參數(shù)和顆粒粘結(jié)參數(shù)。[0102]需要說明的是，顆粒參數(shù)描述顆粒本身的力學(xué)特性和顆粒間的接觸行為；顆粒粘結(jié)參數(shù)描述顆粒間的粘結(jié)效應(yīng)，適用于模擬膠結(jié)材料(如固結(jié)巖石),顆粒間接觸個粘結(jié)的關(guān)系如圖3所示。[0103]在可選的實施方式中，基于Hertz-Mindlin理論建立顆粒間接觸本構(gòu)模型，在本構(gòu)模型中，顆粒間被賦予彈性-摩擦接觸特性；定義顆粒接觸時的力學(xué)響應(yīng)，包括法向正應(yīng)力fn=kδ和顆粒間剪應(yīng)力fs=kδs,其中法向正應(yīng)力與顆粒間的法向剛度k和顆粒間的重疊量δ相關(guān)；顆粒間剪應(yīng)力與顆粒間的剪切剛度ks和顆粒中心的剪切偏移量δs相關(guān)；kn和ks值都是與重疊顆粒接觸面積有關(guān)的非線性量，與顆粒剪切模量Gp和泊松比p有關(guān)。[0104]進一步的，確定離散元相鄰顆粒間的狀態(tài)關(guān)系，并通過三軸壓縮實驗測試校準(zhǔn)各個地層的顆粒材料細(xì)觀參數(shù)，其中狀態(tài)關(guān)系包括粘結(jié)壓縮狀態(tài)、粘結(jié)拉伸狀態(tài)、無粘結(jié)接觸狀態(tài)以及顆粒失去相互作用狀態(tài)。[0105]示例說明的，三軸壓縮實驗測試采用半徑膨脹法建立1:1:2的初始模型，在初始模型中，頂?shù)變擅鎵σ韵嗤乃俣认蛟嚇又行囊苿?，周邊設(shè)置不同的圍壓；獲取三軸壓縮實驗測試試樣在應(yīng)變20%時的應(yīng)力應(yīng)變曲線，通過對應(yīng)力應(yīng)變曲線進行后處理，得到試樣的粘聚力、摩擦系數(shù)等參數(shù)。建立顆粒細(xì)觀參數(shù)與地層宏觀性質(zhì)的對應(yīng)關(guān)系。[0106]更進一步的，基于目標(biāo)區(qū)域相應(yīng)的空間，隨機分布半徑呈正態(tài)分布的顆粒，并建立顆粒細(xì)觀參數(shù)與地層宏觀性質(zhì)的對應(yīng)關(guān)系，賦予每個顆粒物理屬性和相互作用參數(shù)，其中物理屬性包括半徑和密度；相互作用參數(shù)包括法向剛度、剪切剛度和摩擦系數(shù)。將賦予后的顆粒自由沉降至系統(tǒng)穩(wěn)定，記錄各個地層顆粒的位置參數(shù)，完成第一離散元模型的構(gòu)建。[0107]具體的，基于標(biāo)定的顆粒粘結(jié)參數(shù)和目標(biāo)區(qū)域的地層信息，在第一離散元模型的顆粒之間設(shè)定粘結(jié)關(guān)系。[0108]需要說明的是，粘結(jié)關(guān)系通過粘結(jié)強度、抗拉強度和內(nèi)摩擦角的參數(shù)定義顆粒間粘結(jié)的失效準(zhǔn)則。[0109]進一步的，對第一離散元模型進行細(xì)化，建立目標(biāo)區(qū)域不同層位顆粒的細(xì)觀參數(shù)和力學(xué)關(guān)系；根據(jù)斷層點集數(shù)據(jù)Fault,建立主要斷層模型函數(shù)isInsideFault,其中主要斷層模型函數(shù)用于判斷任一顆粒坐標(biāo)pt是否位于斷層點集Fault內(nèi)部。[0110]需要說明的是，當(dāng)顆粒坐標(biāo)pt位于斷層點集Fault所定義的斷層面內(nèi)部時，則主要斷層模型函數(shù)isInsideFault返回真值，當(dāng)顆粒坐標(biāo)pt位于斷層點集Fault所定義的斷層面外部時，則主要斷層模型函數(shù)isInsideFault返回假值。[0111]更進一步的，基于判斷結(jié)果，取消斷層面上顆粒間的粘結(jié)，同時降低斷層顆粒間摩擦系數(shù)，形成第二離散元模型。[0112]S4:基于第二離散元模型，設(shè)置邊界條件進行逆沖推覆模擬，并通過第二離散元模型的輸出結(jié)果與實際地質(zhì)剖面的對比，調(diào)整模型參數(shù)，直至二者匹配。[0113]具體的，設(shè)置邊界條件進行逆沖推覆模擬，包括：將第二離散元模型的后緣設(shè)置為活動邊界，其余邊界設(shè)置為固定邊界；參考目標(biāo)區(qū)域的斷層位移量和平衡剖面恢復(fù)縮短量，對活動邊界施加運動；控制活動邊界的運動，當(dāng)達到與平衡剖面恢復(fù)縮短量相等的擠壓量時，則停止第二離散元模型的運算。[0114]進一步的，提取第二離散元模型的模擬結(jié)果，將模擬結(jié)果與實際地質(zhì)剖面進行一致性對比；若模擬結(jié)果與實際地質(zhì)剖面不匹配，則微調(diào)離散元顆粒參數(shù)或邊界運動情況，重復(fù)上述步驟，直至模擬結(jié)果與實際地質(zhì)剖面匹配；若模擬結(jié)果與實際地質(zhì)剖面匹配，則執(zhí)行三維斷層/裂縫網(wǎng)絡(luò)構(gòu)建。[0115]S5:基于匹配后的模型結(jié)果，建立并分析目標(biāo)區(qū)域的三維斷層和裂縫網(wǎng)絡(luò)的空間展布特征，并對三維斷層的活動性進行動態(tài)評估與趨勢預(yù)測。[0116]具體的，基于第二離散元模型演化結(jié)果的顆粒坐標(biāo)，根據(jù)顆粒間相對位置向量，推導(dǎo)出不同顆粒位移梯度張量；根據(jù)第二離散元模型演化結(jié)果的顆粒位移梯度張量，繪制不同剖面的應(yīng)變張量，得到空間展布的三維斷層。[0117]進一步的，三維斷層包括：將第二離散元模型的空間劃分為每個顆粒的控制體積，其中顆粒i在初始位置為x,當(dāng)前位置為x!,鄰居顆粒j的位置為x9→x;為每個顆粒分配一個單元，其中單元包括所有到此顆粒最近的空間點集合V,={x∈R3|Px-x,RPx-x,P,Vj≠i},滿足條件：對于所有顆粒點j∈N。[0118]更進一步的，基于顆粒在參考狀態(tài)和當(dāng)前狀態(tài)的位置，計算位移向量和相對位置[0121]其中，△uij為顆粒j相對于顆粒i的位移向量，△Xij為顆粒j相對于顆粒i的位置向為顆粒j在初始參考狀態(tài)下的位置坐標(biāo)，x(1)為顆粒i在當(dāng)前狀態(tài)下的位置坐標(biāo)，x(0)為顆粒i在初始參考狀態(tài)下的位置坐標(biāo)。[0125]優(yōu)選的，根據(jù)位移梯度張量計算應(yīng)變張量，繪制不同剖面的應(yīng)變張量分布，得到空間展布的三維斷層。[0128]其中，E為顆粒i的應(yīng)變張量，H?為顆粒i的位移梯度張量，H為顆粒i的位移梯度張量的轉(zhuǎn)置。[0129]具體的，基于第二離散元模型演化結(jié)果中各顆粒間的粘結(jié)關(guān)系，判斷顆粒間粘結(jié)是否發(fā)生破壞。[0130]需要說明的是，若應(yīng)力值大于預(yù)設(shè)粘結(jié)強度，則判定顆粒間粘結(jié)發(fā)生破壞；若應(yīng)力值小于或等于預(yù)設(shè)粘結(jié)強度，則判定顆粒間粘結(jié)未發(fā)生破壞。[0131]進一步的，根據(jù)判斷結(jié)果記錄顆粒間粘結(jié)破壞的時間、位置和類型，其中類型包括張破裂和剪破裂；將顆粒間粘結(jié)狀態(tài)標(biāo)記為裂縫產(chǎn)生點，其中裂縫產(chǎn)生點為顆粒間粘結(jié)狀態(tài)由粘結(jié)壓縮狀態(tài)或粘結(jié)拉伸狀態(tài)轉(zhuǎn)變?yōu)闊o粘結(jié)接觸狀態(tài)或失去相互作用狀態(tài)的情況。[0132]更進一步的，對張破裂類型和剪破裂類型進行特征分析，其中張破裂類型表現(xiàn)為垂直于最大主應(yīng)力方向上的破壞；剪破裂類型表現(xiàn)為沿著最小主應(yīng)力方向發(fā)展的破壞。[0133]具體的，基于加載過程的推進，持續(xù)更新顆粒間粘結(jié)破壞的信息，建立裂縫網(wǎng)絡(luò)數(shù)據(jù)庫，其中裂縫網(wǎng)絡(luò)數(shù)據(jù)庫包括已發(fā)生粘結(jié)破壞的顆粒坐標(biāo)、破壞時間和破壞類型；基于裂縫網(wǎng)絡(luò)數(shù)據(jù)庫，進行裂縫網(wǎng)絡(luò)定量分析，計算統(tǒng)計指標(biāo)，其中統(tǒng)計指標(biāo)包括裂縫密度、平均裂縫長度和裂縫方向分布；[0134]進一步的，通過統(tǒng)計指標(biāo)繪制不同時空演化下裂縫網(wǎng)絡(luò)的空間展布特征圖。[0135]更進一步的，參照圖5,本實施例還提供一種基于離散元的逆沖推覆體系下三維斷[0136]斷層模塊，基于目標(biāo)區(qū)域的地震測井資料參數(shù)，構(gòu)建初始化斷層面模型，并采用加密網(wǎng)格化方法沿初始化斷層面模型的走向方向提取基本斷層點集數(shù)據(jù)；[0137]反演模塊，用于通過觀測目標(biāo)區(qū)域的地質(zhì)并解析該地區(qū)地震資料，提取地層信息與變形過程，并對地震解析圖進行反演，斷層的空間演化過程；[0138]模型模塊，根據(jù)地層參數(shù)標(biāo)定離散元細(xì)觀參數(shù)，沉積第一離散元模型，構(gòu)建斷層函[0139]變形模塊，基于第二離散元模型，設(shè)置邊界條件進行逆沖推覆模擬，并通過第二離散元模型的輸出結(jié)果與實際地質(zhì)剖面的對比，調(diào)整模型參數(shù)，直至二者匹配；[0140]建模模塊，基于匹配后的模型結(jié)果，建立并分析目標(biāo)區(qū)域的三維斷層和裂縫網(wǎng)絡(luò)的空間展布特征，并對三維斷層的活動性進行動態(tài)評估與趨勢預(yù)測。[0141]本實施例還提供一種計算機設(shè)備，適用于基于離散元的逆沖推覆體系下三維斷層/裂縫網(wǎng)絡(luò)建模方法的情況，包括存儲器和處理器；存儲器用于存儲計算機可執(zhí)行指令，處理器用于執(zhí)行計算機可執(zhí)行指令，實現(xiàn)如上述實施例提出的基于離散元的逆沖推覆體系下三維斷層/裂縫網(wǎng)絡(luò)建模方法。[0142]該計算機設(shè)備可以是終端，該計算機設(shè)備包括通過系統(tǒng)總線連接的處理器、存儲該計算機設(shè)備的存儲器包括非易失性存儲介質(zhì)、內(nèi)存儲器。該非易失性存儲介質(zhì)存儲有操作系統(tǒng)和計算機程序。該內(nèi)存儲器為非易失性存儲介質(zhì)中的操作系統(tǒng)和計算機程序的運行提供環(huán)境。該計算機設(shè)備的通信接口用于與外部的終端進行有線或無線方式的通信，無線方式可通過WIFI、運營商網(wǎng)絡(luò)、NFC(近場通信)或其他技術(shù)實現(xiàn)。該計算機設(shè)備的顯示屏可以是液晶顯示屏或者電子墨水顯示屏，該計算機設(shè)備的輸入裝置可以是顯示屏上覆蓋的觸摸層，也可以是計算機設(shè)備外殼上設(shè)置的按鍵、軌跡球或觸控板，還可以是外接的鍵盤、觸控板或鼠標(biāo)等。[0143]本實施例還提供一種存儲介質(zhì)，其上存儲有計算機程序，該程序被處理器執(zhí)行時實現(xiàn)如上述實施例提出的實現(xiàn)基于離散元的逆沖推覆體系下三維斷層/裂縫網(wǎng)絡(luò)建模方[0144]本實施例提出的存儲介質(zhì)與上述實施例提出的數(shù)據(jù)存儲方法屬于同一發(fā)明構(gòu)思，未在本實施例中詳盡描述的技術(shù)細(xì)節(jié)可參見上述實施例，并且本實施例與上述實施例具有相同的有益效果。[0145]綜上，本發(fā)明通過引入離散元方法對斷層裂縫系統(tǒng)進行動態(tài)模擬，有效規(guī)避了傳統(tǒng)建模中因網(wǎng)格重構(gòu)困難導(dǎo)致的斷層尖端破裂擴展行為難以精準(zhǔn)刻畫的問題；通過在建模過程中融合地質(zhì)力學(xué)參數(shù)與構(gòu)造演化過程，建立了先存斷層與新生裂縫之間的時空耦合關(guān)系，實現(xiàn)了多期構(gòu)造疊加效應(yīng)的量化解析；同時建立了從構(gòu)造跡象反演到未來構(gòu)造演化預(yù)測的閉環(huán)建模流程，顯著提升了復(fù)雜斷裂系統(tǒng)三維空間展布建模的地質(zhì)合理性與預(yù)測能力，為油氣勘探中隱伏構(gòu)造識別、斷層封堵性評價及構(gòu)造控藏機制研究提供了科學(xué)、可靠的技術(shù)支撐。[0147]參照圖4、圖6~圖12,為本發(fā)明第二個實施例，該實施例提供了一種基于離散元的逆沖推覆體系下三維斷層/裂縫網(wǎng)絡(luò)建模方法，為了驗證本發(fā)明的有益效果，通過經(jīng)濟效益計算和仿真實驗進行科學(xué)論證。[0148]以塔里木盆地西部巴楚隆起南緣中段為研究對象，選取其內(nèi)部分布有鳥山、瑪扎塔格與羅斯塔格三大構(gòu)造帶交匯的典型逆沖推覆構(gòu)造區(qū)作為目標(biāo)區(qū)域。該區(qū)域因多期構(gòu)造疊加改造，斷裂系統(tǒng)密集且交錯復(fù)雜，具備典型的推覆變形特征。多個構(gòu)造帶之間的擠壓匯聚不僅導(dǎo)致斷裂系統(tǒng)呈現(xiàn)多樣性與非均質(zhì)性，同時也孕育了豐富的油氣資源，明確深部斷裂展布及構(gòu)造樣式對于提升油氣勘探效率具有重要意義。[0149]鳥山構(gòu)造帶的成因可視為東側(cè)瑪扎塔格與西側(cè)羅斯塔格構(gòu)造帶后期推進速率差異演化的結(jié)果，其結(jié)構(gòu)演化主要依據(jù)瑪扎塔格和羅斯塔格構(gòu)造帶的地震資料進行模擬推斷。為此，首先選取高分辨率地震剖面，對其進行預(yù)處理和斷層識別，提取清晰的斷層反射特征。如圖8所示，在瑪扎塔格構(gòu)造帶內(nèi)識別出一條向北逆沖的主斷層F1,并觀察到沿該主斷層發(fā)育的多個分支反沖斷層。針對每條斷層，在其上、下盤反射軸突然中斷或錯位的位置，提取空間坐標(biāo)信息(包括橫向、縱向位置及時間/深度),將二維地震信息轉(zhuǎn)換為三維坐標(biāo)點集。最終，分別構(gòu)建羅斯塔格斷層點集p?和瑪扎塔格斷層數(shù)據(jù)點p?,為后續(xù)離散元建模提供結(jié)構(gòu)邊界依據(jù)。[0150]為了準(zhǔn)確構(gòu)建地層力學(xué)特性與破裂演化過程，進一步采用三軸壓縮實驗標(biāo)定顆粒的細(xì)觀參數(shù)。如圖7所示，試驗中采用半徑膨脹法構(gòu)建樣本，并通過頂?shù)變擅鎵Φ人賶嚎s加載，同時在5~80MPa的圍壓下分別測試樣本應(yīng)力-應(yīng)變響應(yīng)，記錄20%應(yīng)變處的曲線，結(jié)果表明：模型顆粒的宏觀響應(yīng)粘聚力約為32.8M