【文章內(nèi)容簡(jiǎn)介】 這樣在垂向上共形成了若干節(jié)點(diǎn)（ 6個(gè)節(jié)點(diǎn) ），求解每個(gè)節(jié)點(diǎn)上的超壓 Pai， i=1,2,…6 。若節(jié)點(diǎn) i 處于正常壓實(shí)帶，取 Pai=0。 取第二個(gè)時(shí)間步長(zhǎng) Δt ，該階段內(nèi) 又沉積 有一定厚度的地層，根據(jù)空間步長(zhǎng) Δz ，將 新沉積 的地層在垂向上 從下到上 劃分為若干小段 (假設(shè) 4小段 )，這樣在垂向上增加了若干節(jié)點(diǎn) (5個(gè)節(jié)點(diǎn) )，求解 已有的 每個(gè)節(jié)點(diǎn)上的超壓 Pai，i=1,2,…11 。若節(jié)點(diǎn) i 處于正常壓實(shí)帶，取 Pai=0。 取第三個(gè)時(shí)間步長(zhǎng) Δt ，該階段內(nèi) 又沉積 有一定厚度的地層，根據(jù)空間步長(zhǎng) Δz ，將新沉積的地層在垂向上 從下到上 劃分為若干小段 (假設(shè) 3小段 )，在垂向上又增加了若干節(jié)點(diǎn) (4個(gè)節(jié)點(diǎn) )，求解所有每個(gè)節(jié)點(diǎn)上的超壓 Pai，i=1,2,…15 。若節(jié)點(diǎn) i處于正常壓實(shí)帶，取 Pai=0。 … … … 取最后一個(gè)時(shí)間步長(zhǎng) Δt ，該階段內(nèi) 又沉積 有一定厚度的地層，根據(jù)空間步長(zhǎng) Δz ，將新沉積的地層在垂向上從下到上 劃分為若干小段 (假設(shè) 5小段 )，在垂向上又增加了若干節(jié)點(diǎn) (6個(gè)節(jié)點(diǎn) )，求解所有每個(gè)節(jié)點(diǎn)上的超壓 Pai，i=1,2,… n。若節(jié)點(diǎn) i處于正常壓實(shí)帶，取 Pai=0。其中 n為節(jié)點(diǎn)的總數(shù)。 超壓史求解完畢 。 在模擬過程中，正常壓實(shí)帶和欠壓實(shí)帶的分界線的確定方法： 在 模擬之前 先確定盆地 出現(xiàn)欠壓實(shí)的層位 ，由此可以確定平面上每個(gè)模擬點(diǎn)欠壓實(shí)出現(xiàn)的 深度 ，當(dāng)節(jié)點(diǎn)處于 該深度之下 時(shí)，認(rèn)為該節(jié)點(diǎn) 處于欠壓實(shí)帶 ，否則認(rèn)為該節(jié)點(diǎn)處于正常壓實(shí)帶。 求解超壓方程使用有限差分方法 六、 超壓方程的數(shù)值解法 (一維 ) 求解超壓方程一般使用有限差分等數(shù)值解法，由于二維和三維模擬在解方程時(shí)方法繁瑣，而且就一般的微機(jī)而言，進(jìn)行多維模擬的能力也不夠。因此，僅就一維超壓方程討論其有限差分?jǐn)?shù)值解法。 ? ? ? ??? ? ? ? ? ? ? ? ? ? ? ? ?? ? ? ? ? ?z K P z P t St P tw a w a w w n( ) ( ) ( ) ( )? ? ? ? ? ? ? ? ?1 1 1? ? ? ? ?一維超壓方程可表示為： 超壓方程的邊界條件可考慮為： 考慮盆地的基底為致密巖石組成的不滲透地層 ， 故下邊界條件為： (Z = 0 ) 上邊界為正常壓實(shí)和欠壓實(shí)的分界線 ， 故上邊界 條件為： (Z = L) ??Pza ? 0Pa ? 0因此，一維超壓方程的定解問題為： ? ?? ? ? ? ? ?? ???????? ? ???? ???? ? ?????zK PZPtStPtPtZP Z Lw awaw wnaa?????? ? ? ? ? ? ? ? ?? ?? ????????????1 1 10 00? ? ? ? ?( )( )求解上述超壓方程使用有限差分方法。 最終求解出每個(gè)縱向節(jié)點(diǎn)上的超壓值。 第三節(jié) 反演法模型（地層回剝技術(shù)） 主要思路： 各地層在保持其骨架厚度不變的條件下 ， 從今天盆地的分層現(xiàn)狀出發(fā) ， 按地質(zhì)年代 從新到老 (由今至古 )逐層剝?nèi)?， 直至全部剝完 。 回剝技術(shù)適用于 正常壓實(shí)帶 ，即超壓為零的情況。它所用的關(guān)鍵參數(shù)是 孔隙度 深度曲線 。而孔隙度 深度曲線應(yīng)滿足以下三個(gè)條件： (1)、 各地層應(yīng)有各自的孔隙度 深度曲線 ， 主要是因?yàn)椴煌貙拥膸r性和結(jié)構(gòu)存在差異 。 (2)、 用今天實(shí)測(cè)的孔隙度 深度曲線代替古孔隙度 深度曲線 ， 但今測(cè)的孔隙度 深度曲線可能受到各種地質(zhì)事件 (如剝蝕 、 斷層等 )的影響 ， 不能真實(shí)地反映正常壓實(shí)下的孔隙度隨深度的變化規(guī)律 ， 因此必須盡量消除這些影響 ， 獲得正常壓實(shí)狀況下的孔隙度 深度曲線 。 (3)、各地層應(yīng)有三種巖性的孔隙度 深度曲線 (砂巖、泥巖、灰?guī)r )，其中砂巖假定為除泥巖和灰?guī)r外的所有其它巖性的總和。 一、骨架厚度公式和地層底界公式的提出 隨著埋藏深度的增加，地層的上覆負(fù)載也增加，導(dǎo)致孔隙度變小，體積變小。若假定地層的橫向在沉積過程中保持不變，則地層被壓實(shí)僅僅是縱向上變化，地層體積的變小就歸結(jié)為地層厚度的變小 。另外，根據(jù)巖石骨架不可壓縮的假設(shè)，地層的骨架厚度 (孔隙度為零時(shí)的地層厚度 )始終不變，除非發(fā)生剝蝕等地質(zhì)事件。地層的骨架厚度始終小于實(shí)際的地層厚度。計(jì)算地層骨架厚度的公式為： ? ?? ???21)(1zzs dzzh其中： hs ～地層的骨架厚度 (m) ； z1 ～地層的頂界深度 (m)； z2 ～地層的底界深度 (m)； Φ(z)～地層的孔隙度 深度曲線函數(shù) (小數(shù) )。 Φ(z)的通式可表示為： ? ? ? ?( ) ( ) ( ) ( )z P z P z P zs s m m l l? ? ?其中： Ps ～地層的砂巖含量 (小數(shù) )； Pm ～地層的泥巖含量 (小數(shù) )； Pl ～地層的灰?guī)r含量 (小數(shù) )； Φs(z)～ 地層砂巖孔隙度 深度曲線函數(shù)； Φm(z)～ 地層泥巖孔隙度 深度曲線函數(shù)； Φl (z) ～ 地層灰?guī)r孔隙度 深度曲線函數(shù)； Ps + Pm + Pl =1 Φ(z)表示為砂巖、泥巖、灰?guī)r三種巖性孔 深曲線函數(shù)的 加權(quán)平均值 。 一般，砂巖、泥巖、灰?guī)r三種巖性的孔隙度 深度曲線函數(shù)形式如下： 砂巖： Φ0s ～砂巖的初始 (z=0)孔隙度 (小數(shù) )， cs ～常數(shù) (m1 ) 。 泥巖： Φ0 m ～泥巖的初始 (z=0)孔隙度 (小數(shù) )， cm ～常數(shù) (m1 ) 。 灰?guī)r： Φ0l ～灰?guī)r的初始 (z=0)孔隙度 (小數(shù) )， cl ～常數(shù) (m1 ) 。 按三種巖性含量加權(quán)平均得： ? ?s s c zz e s( ) ? ?0? ?m m c zz e m( ) ? ?0? ?l l c zz e l( ) ? ?0? ? ? ?( )z P e P e P es s c z m m c z l l c zs m l? ? ?? ? ?0 0 0將上式代入骨架厚度公式，并整理得： 上式是考慮三種巖性 (砂巖、泥巖、灰?guī)r )的骨架厚度公式。是 回剝技術(shù)中重要的公式之一 。如果地層沒有發(fā)生剝蝕，則根據(jù)目前資料求出的各地層的骨架厚度在以前的各地質(zhì)時(shí)刻保持不變。 ? ?? ? ? ?12121200012zczclllzczcmmmzczcssssllmmsseecPeecPeecPzzh?????????????????將骨架厚度公式變化位置可得：