一種獲取地層孔隙度的方法及裝置的制造方法
【技術(shù)領域】
[0001] 本發(fā)明屬于石油勘探資料處理領域，特別涉及獲取地層孔隙度的方法及裝置。
【背景技術(shù)】
[0002] 地層孔隙度是儲層評價的重要參數(shù)之一，地層孔隙度與測井電位、鉆井時間、油氣 儲量以及含油飽和度等石油勘探參數(shù)緊密相連，因此，獲取地層孔隙度是石油勘探的一個 重要工作步驟。現(xiàn)有技術(shù)中，通常根據(jù)電纜聲波測得的時差計算得到所述地層孔隙度。
[0003] 在實現(xiàn)本申請過程中，發(fā)明人發(fā)現(xiàn)現(xiàn)有技術(shù)中至少存在如下問題：在很多地層中， 尤其在大斜度井以及水平井中，存在時差各向異性現(xiàn)象。這時電纜聲波測得的時差是水平 時差和垂直時差的綜合值。而各向異性現(xiàn)象多數(shù)是由于地層內(nèi)夾雜著具有不同時差的其他 地層而產(chǎn)生的。計算所述孔隙度所需要的時差一般是水平時差，因此使用電纜聲波測得的 時差計算得到的地層孔隙度偏離真實地層孔隙度，進而根據(jù)所述地層孔隙度計算得到的含 油飽和度等參數(shù)值誤差也相對較大。

【發(fā)明內(nèi)容】

[0004] 針對現(xiàn)有技術(shù)中根據(jù)聲波測井儀器測得的時差計算得到的孔隙度誤差較大，本發(fā) 明提出一種獲取地層孔隙度的方法及裝置，所述獲取地層孔隙度的方法及裝置具體是這樣 實現(xiàn)的：
[0005] 一種獲取地層孔隙度的方法，所述方法包括：
[0006] 獲取目標地層的聲波測井資料、井斜角曲線資料以及層位時差各向異性比；
[0007] 根據(jù)所述目標地層的井斜角曲線資料以及層位時差各向異性比，按照擬合方程計 算所述目標地層在不同井斜角、不同層位時差各向異性比條件下的慢度差，所述擬合方程 根據(jù)已知的實驗地層訓練得到；
[0008] 根據(jù)所述目標地層聲波測井資料中的時差以及所述目標地層的慢度差計算所述 目標地層的水平時差；
[0009] 根據(jù)所述目標地層的水平時差計算所述目標地層的孔隙度。
[0010] 可選的，在本發(fā)明一實施例中，所述根據(jù)已知的實驗地層訓練得到擬合方程，包 括：
[0011] 獲取所述實驗地層的水平時差，計算所述實驗地層在不同井斜角、不同時差各向 異性比條件下的時差；
[0012] 根據(jù)所述實驗地層的水平時差，計算所述實驗地層的慢度差，所述慢度差為所述 時差與所述實驗地層的水平時差的差值；
[0013] 建立所述井斜角以及所述慢度差的坐標系，在所述坐標系中插入所述不同井斜 角、不同時差各向異性比條件下的所述慢度差的離散點；
[0014] 將所述離散點擬合成曲線，并獲取所述曲線的擬合方程。
[0015] 可選的，在本發(fā)明一實施例中，所述計算所述實驗地層在不同井斜角、不同時差各 向異性比條件下的時差，包括：
[0016] 利用聲波測井儀采用有限差分法計算所述實驗地層在不同井斜角、不同時差各向 異性比條件下的時差，其中，所述井斜角從0到90°范圍內(nèi)按15°步長遞增，所述時差各向 異性比從I. 0到1. 2范圍內(nèi)按0. 02步長遞增。
[0017] 可選的，在本發(fā)明一實施例中，所述聲波測井儀包括雙發(fā)雙收聲波測井儀。
[0018] 可選的，在本發(fā)明一實施例中，所述建立所述井斜角以及所述慢度差的坐標系，具 體包括：
[0019] 建立所述井斜角正弦值以及所述慢度差的直角坐標系。
[0020] 可選的，在本發(fā)明一實施例中，所述擬合方程包括一元二次方程。
[0021] 可選的，在本發(fā)明一實施例中，所述按照所述擬合方程，計算所述目標地層在不同 井斜角、不同層位時差各向異性比條件下的慢度差，具體包括：
[0022] 根據(jù)所述實驗地層井斜角以及慢度差的坐標系中的曲線，按照所述擬合方程，計 算與所述實驗地層相同的時差各向異性比條件下，所述目標地層在不同井斜角條件下的慢 度差，當所述井斜角是步長15°的倍數(shù)時，按照所述擬合方程求解所述目標地層的慢度差； 當所述井斜角不是15°的倍數(shù)時，采用插值法計算所述目標地層的慢度差。
[0023] 可選的，在本發(fā)明一實施例中，根據(jù)所述目標地層的水平時差計算所述目標地層 的孔隙度，其中所述孔隙度的計算公式為：
[0025] 式中，為孔隙度，Δ tM、Δ tf、Δ tsh分別為巖石骨架聲波水平方向時差、地層流 體聲波水平方向時差、地層流體聲波水平方向時差，V sh為地層泥質(zhì)含量，CpS聲波壓實校正 系數(shù)，At為所述目的地層的水平時差。
[0026] -種獲取地層孔隙度的裝置，包括：
[0027] 資料獲取模塊，用于獲取目標地層的聲波測井資料、井斜角曲線資料以及層位時 差各向異性比；
[0028] 目標慢度差計算模塊，用于根據(jù)所述目標地層的井斜角曲線資料以及層位時差各 向異性比，按照擬合方程計算所述目標地層在不同井斜角、不同層位時差各向異性比條件 下的慢度差，所述擬合方程根據(jù)擬合裝置計算得到；
[0029] 水平時差計算模塊，用于根據(jù)所述目標地層聲波測井資料中的時差以及所述目標 地層的慢度差計算所述目標地層的水平時差；
[0030] 孔隙度計算模塊，用于根據(jù)所述目標地層的水平時差計算所述目標地層的孔隙 度。
[0031] 可選的，在本發(fā)明一實施例中，所述目標慢度差計算模塊還包括擬合裝置，所述擬 合裝置包括：
[0032] 實驗時差計算模塊，用于獲取實驗地層的水平時差，計算所述實驗地層在不同井 斜角、不同時差各向異性比條件下的時差；
[0033] 慢度差計算模塊，用于根據(jù)所述實驗地層的水平時差，計算所述實驗地層的慢度 差，所述慢度差為所述時差與所述實驗地層的水平時差的差值；
[0034] 坐標系建立模塊，用于建立所述井斜角以及所述慢度差的坐標系，在所述坐標系 中插入所述不同井斜角、不同時差各向異性比條件下的所述慢度差的離散點；
[0035] 擬合模塊，用于將所述離散點擬合成曲線，并獲取所述曲線的擬合方程。
[0036] 可選的，在本發(fā)明一實施例中，所述實驗時差計算模塊，還包括：
[0037] 時差計算子模塊，用于利用聲波測井儀采用有限差分法計算所述實驗地層在不同 井斜角、不同時差各向異性比條件下的時差，其中，所述井斜角從〇到90°范圍內(nèi)按15°步 長遞增，所述時差各向異性比從1. 〇到1. 2范圍內(nèi)按0. 02步長遞增。
[0038] 可選的，在本發(fā)明一實施例中，所述聲波測井儀包括雙發(fā)雙收聲波測井儀。
[0039] 可選的，在本發(fā)明一實施例中，所述坐標系建立模塊，還包括：
[0040] 坐標系建立子模塊，用于建立所述井斜角正弦值以及所述慢度差的直角坐標系。
[0041] 可選的，在本發(fā)明一實施例中，所述擬合方程包括一元二次方程。
[0042] 可選的，在本發(fā)明一實施例中，所述目標慢度差計算模塊，還包括：
[0043] 目標慢度差計算子模塊，用于根據(jù)所述實驗地層井斜角以及慢度差的坐標系中的 曲線，按照所述擬合方程，計算與所述實驗地層相同的時差各向異性比條件下，所述目標地 層在不同井斜角條件下的慢度差，當所述井斜角是步長15°的倍數(shù)時，按照所述擬合方程 求解所述目標地層的慢度差；當所述井斜角不是15°的倍數(shù)時，采用插值法計算所述目標 地層的慢度差。
[0044] 可選的，在本發(fā)明一實施例中，所述孔隙度