本發(fā)明公開了一種頁巖氣藏壓裂井組生產(chǎn)動態(tài)數(shù)值模擬方法，包括以下步驟：頁巖儲層和壓裂水平井三維地質(zhì)體生成及非結(jié)構(gòu)四面體網(wǎng)格離散；建立頁巖氣藏多級壓裂水平井綜合滲流數(shù)學(xué)模型；基于控制體積有限元法構(gòu)建綜合滲流數(shù)學(xué)模型的全隱式數(shù)值模型；建立基于四面體網(wǎng)格離散的全隱式多級壓裂水平井?dāng)?shù)值井模型，并嵌入全隱式數(shù)值模型，構(gòu)建完整的多系統(tǒng)同步聯(lián)立求解矩陣；對多系統(tǒng)同步聯(lián)立求解矩陣迭代求解，并分析朗格繆爾體積等對頁巖氣藏壓裂井組生產(chǎn)動態(tài)的影響。本發(fā)明彌補(bǔ)了現(xiàn)有頁巖氣藏二維模擬器的不足，綜合考慮頁巖儲層物性層間非均質(zhì)性和水力壓裂裂縫縱向延伸擴(kuò)展特征，實現(xiàn)了三維深層頁巖儲層中立體開發(fā)壓裂井組生產(chǎn)動態(tài)的準(zhǔn)確模擬和預(yù)測。