本發(fā)明公開了一種頁(yè)巖氣藏產(chǎn)能數(shù)值模擬方法，涉及地質(zhì)勘探技術(shù)領(lǐng)域，包括建立非結(jié)構(gòu)四面體網(wǎng)絡(luò)模型；建立基質(zhì)系統(tǒng)流動(dòng)方程；建立天然裂縫骨架體積應(yīng)變和大尺度裂縫骨架體積應(yīng)變的求解方程組；建立頁(yè)巖氣非線性滲流模型，并確定頁(yè)巖氣非線性滲流模型的定解條件；獲得所述基質(zhì)系統(tǒng)流動(dòng)方程、天然裂縫流動(dòng)方程、大尺度裂縫流動(dòng)方程的半離散格式；獲得所述目標(biāo)區(qū)域的水平井產(chǎn)量分布。本發(fā)明公開的頁(yè)巖氣藏產(chǎn)能數(shù)值模擬方法，針對(duì)現(xiàn)有技術(shù)中儲(chǔ)層產(chǎn)量模擬方法中忽略地層骨架體積應(yīng)變的問(wèn)題，提供一種頁(yè)巖氣藏產(chǎn)能數(shù)值模擬方法，在氣藏產(chǎn)能模擬中引入骨架體積應(yīng)變，提高模擬結(jié)果的準(zhǔn)確性。