專利名稱：預測儲層的方法
預測儲層的方法技術(shù)領(lǐng)域
本申請涉及地球物理勘探領(lǐng)域，更具體地講，涉及一種儲層預測方法。
背景技術(shù)：
目前，鮞灘儲層和生物礁儲層是四川盆地的勘探熱點。鮞灘儲層是發(fā)育在三疊系 飛仙關(guān)組的鮞粒狀灘相儲層，礁儲層是發(fā)育在二疊系長興組的與生物礁相關(guān)的灘相儲層。 由于礁灘儲層具有埋藏深、較大非均質(zhì)性、各向異性較強的特點，因此預測鮞灘儲層和礁灘 儲層的勘探風險大、成本高。除位于臺緣區(qū)域附近的鮞灘儲層與生物礁儲層具有較好的物 性，能獲得較高的產(chǎn)能外，臺內(nèi)礁灘儲層受泥巖的影響，預測準確率較低。例如，在四川盆 地，對礁灘儲層的量化預測還僅處在利用常規(guī)的疊后速度反演(波阻抗)和伽瑪反演對儲 層厚度進行預測，不能指示真正的礁灘開發(fā)有利區(qū)。因此，需要一種既能有效預測鮞灘儲層 和礁灘儲層，又能排除其它影響因素的方法。發(fā)明內(nèi)容
根據(jù)本發(fā)明的示例性實施例，提供了一種預測儲層的方法，所述方法包括(a)采 集待測地區(qū)的地震數(shù)據(jù)以及待測地區(qū)的測井數(shù)據(jù)；(b)對測井數(shù)據(jù)進行頻譜分析，以確定 對儲層敏感的儲層敏感頻率和對泥巖層敏感的泥巖層敏感頻率；(C)對采集的地震數(shù)據(jù)進 行分頻處理，以獲取頻域地震數(shù)據(jù)；(d)從所述頻域地震數(shù)據(jù)獲取待測地區(qū)的特定地質(zhì)層 位處的與不同頻率對應的分頻振幅切片；(e)將所述特定地質(zhì)層位的與儲層敏感頻率對應 的分頻振幅切片的強振幅區(qū)域的分布確定為所述特定地質(zhì)層位的儲層分布，并將所述特定 地質(zhì)層位的與泥巖層敏感頻率對應的分頻振幅切片的強振幅區(qū)域的分布確定為所述特定 地質(zhì)層位的泥巖層分布。
所述的方法還可包括在對采集的地震數(shù)據(jù)進行分頻處理之前，對所述地震數(shù)據(jù) 進行地震資料處理，其中，可在上述步驟(C)中對經(jīng)過地震資料處理的地震數(shù)據(jù)進行分頻 處理來獲取待測地區(qū)的頻域地震數(shù)據(jù)。
所述地震資料處理可包括對采集的地震數(shù)據(jù)進行動校正和靜校正中的至少一 個。
所述分頻處理的步驟可包括對采集的地震數(shù)據(jù)進行離散傅里葉變換或小波變 換。
在上述步驟(b)中對測井數(shù)據(jù)進行頻譜分析的步驟可包括對測井數(shù)據(jù)中的聲波 測井曲線和電阻率曲線進行廣義S變換頻譜分析。
有益效果
通過利用根據(jù)本發(fā)明的預測儲層的方法，能夠有效預測礁灘儲層，同時排除泥巖 對儲層預測的影響，從而能夠預測真正的儲層發(fā)育區(qū)。


通過下面對結(jié)合示例性實施例示出的附圖進行的描述，本發(fā)明的上述和其他目的 和特點將會變得更加清楚，其中
圖1是示出根據(jù)本發(fā)明的示例性實施例的預測儲層的方法的流程圖2是示出根據(jù)本發(fā)明的方法預測的三疊系飛仙關(guān)組鮞灘儲層分布的示例；
圖3是示出根據(jù)本發(fā)明的方法預測的三疊系飛仙關(guān)組泥巖層分布的示例。
以下通過參考附圖描述實施例以解釋本發(fā)明。
具體實施方式
在下文中，除非另有定義，否則這里使用的所有術(shù)語(包括技術(shù)術(shù)語和科學術(shù)語) 具有與本發(fā)明所屬領(lǐng)域的普通技術(shù)人員通常理解的含義相同的含義。還將理解，除非這里 確切地定義，否則術(shù)語(諸如在常用詞典中定義的那些術(shù)語)應被解釋為具有與所述術(shù)語 在相關(guān)領(lǐng)域的上下文中的含義一致的含義，而不應被理想化或過于正式地解釋。
圖1是示出根據(jù)本發(fā)明的示例性實施例的預測儲層的方法的流程圖。
參照圖1，在步驟101，首先需要采集待測地區(qū)的地震數(shù)據(jù)以及待測地區(qū)的測井數(shù) 據(jù)。具體地講，通常為了獲得某一地區(qū)的目的層儲層分布情況，可在該地區(qū)進行地震勘探， 人工激發(fā)地震波，接收反射縱波信息，并通過對采集的地震信息進行處理解釋來獲得與不 同巖石的彈性或波阻抗差異有關(guān)的地震數(shù)據(jù)。同時，對于該區(qū)已有的鉆井，利用井下儀器 沿井身測量地層的各物理參數(shù)，包括聲波、密度、電阻率等參數(shù)，從而確定儲層等地質(zhì)層的 深度或厚度(例如，確定儲層分布在地表以下多少米處)并獲得相應的地球物理參數(shù)數(shù)據(jù) (如儲層的厚度、孔隙度等)。然而，通過測井的方式來獲得地表下的地質(zhì)分布情況(即，地 表以下不同地質(zhì)層的深度或厚度)需要花費較大的時間、人力物力成本，并且還只是一孔 之見，不能查明待測地區(qū)的目的層的平面展布。因此，在本發(fā)明中，可通過在待測地區(qū)進行 地震勘探來獲得該地區(qū)的地震數(shù)據(jù)，并基于通過測井獲得的有效測井結(jié)果，使用在待測地 區(qū)采集的地震數(shù)據(jù)來預測整個待測地區(qū)的地質(zhì)分布情況，尤其是目的層段的儲層和泥巖層 的分布情況。以下將具體說明預測的方法。
在步驟103，可通過對在步驟101中采集的測井數(shù)據(jù)進行頻譜分析，從而確定對儲 層敏感的儲層敏感頻率和對泥巖層敏感的泥巖層敏感頻率。具體地講，可對測井數(shù)據(jù)中的 聲波測井曲線和電阻率曲線進行廣義S變換頻譜分析。然后，在頻域中，可將井區(qū)中探測到 儲層的地質(zhì)層位的聲波測井曲線和電阻率曲線的主頻確定為對儲層敏感的儲層敏感頻率。 此外，還可將井區(qū)中探測到泥巖層的地質(zhì)層位的聲波測井曲線和電阻率曲線的主頻確定為 對泥巖層敏感的泥巖層敏感頻率。應該理解，在某一井區(qū)的同一地質(zhì)層位，相同巖性的電阻 率曲線與聲波測井曲線的主頻通常是一致的。當電阻率曲線與聲波測井曲線的主頻不一致 時，還可使用其他測井數(shù)據(jù)來確定對儲層敏感的頻率和對泥巖層敏感的頻率，例如，還可使 用伽瑪曲線或橫波曲線來確定對儲層敏感的頻率和對泥巖層敏感的頻率。此外，應該理解， 可在待測地區(qū)鉆多個鉆井，并通過對所述多個鉆井的測井數(shù)據(jù)進行上述分析，以更準確地 確定對儲層敏感的儲層敏感頻率和對泥巖層敏感的泥巖層敏感頻率。此外，對測井數(shù)據(jù)進 行頻譜分析的方法對于本領(lǐng)域技術(shù)人員而言是已知的，因此為了簡明，在此將不再進行詳 細描述。
然后，在步驟105，可對在該待測地區(qū)采集的地震數(shù)據(jù)進行分頻處理，以獲取待測地區(qū)的頻域地震數(shù)據(jù)。具體地講，可對采集的待測地區(qū)的時域地震數(shù)據(jù)進行離散傅里葉變 換或小波變換，從而將時域的地震數(shù)據(jù)轉(zhuǎn)換為頻域地震數(shù)據(jù)。此外，在對采集的地震數(shù)據(jù)進 行分頻處理之前，必須對所述地震數(shù)據(jù)進行地震資料處理(例如，對所述地震數(shù)據(jù)進行動 校正、靜校正等中的至少一項處理)，從而使得經(jīng)過處理的地震數(shù)據(jù)能夠更真實的反映待測 地區(qū)的地質(zhì)構(gòu)造。在對采集的地震數(shù)據(jù)進行地震資料處理之后，可對經(jīng)過地震資料處理的 地震數(shù)據(jù)進行分頻處理來獲取待測地區(qū)的頻域地震數(shù)據(jù)，從而有利于獲得更為準確的預測 結(jié)果。應該理解，對地震數(shù)據(jù)進行地震資料處理的方法以及對地震數(shù)據(jù)進行分頻處理的方 法對于本領(lǐng)域技術(shù)人員而言是已知的，因此為了簡明，在此將不進行詳細描述。此外，應該 理解，盡管以上描述中在步驟103之后執(zhí)行步驟105，但步驟105也可在步驟103之前執(zhí)行。
接下來。在步驟107，可從待測地區(qū)的頻域地震數(shù)據(jù)獲取該待測地區(qū)的特定地質(zhì)層 位處的與不同頻率對應的分頻振幅切片。所述分頻振幅切片是指在特定頻率下，待測地區(qū) 的特定地質(zhì)層位的頻域地震數(shù)據(jù)的振幅分布。例如，假設(shè)在步驟105中獲得了某一地區(qū)的 O-1OOHz的頻域地震數(shù)據(jù)，則對于該地區(qū)的特定地質(zhì)層位平面，可將O-1OOHz的頻域地震數(shù)據(jù)中與該地質(zhì)層位平面對應的lHz、2Hz.....99Hz和IOOHz的頻域地震數(shù)據(jù)的振幅分布情況分別確定為該地質(zhì)層位平面的lHz、2Hz.....99Hz和IOOHz的分頻振幅切片(圖2和圖3中示出分頻振幅切片的示例，稍后將進行解釋)。
在步驟109，在待測地區(qū)的特定地質(zhì)層位處(即，待測地區(qū)地表下的處于該特定地 質(zhì)層位的平面)，可將該特定地質(zhì)層位的與儲層敏感頻率對應的分頻振幅切片的強振幅區(qū) 域的分布確定為該特定地質(zhì)層位的儲層分布，并將該特定地質(zhì)層位的與泥巖層敏感頻率對 應的分頻振幅切片的強振幅區(qū)域的分布確定為該特定地質(zhì)層位的泥巖層分布。
圖2至圖3詳細示出了通過分頻振幅切片來預測儲層和泥巖層分布的示例。圖2 是示出根據(jù)本發(fā)明的方法預測的三疊系飛仙關(guān)組鮞灘儲層分布的示例。圖3是示出根據(jù)本 發(fā)明的方法預測的三疊系飛仙關(guān)組泥巖層分布的示例。假設(shè)在步驟103中確定在待測地區(qū) 中對儲層敏感的儲層敏感頻率為18Hz，對泥巖層敏感的泥巖層敏感頻率為40Hz。在圖2中， 標定為龍崗I的井區(qū)的實際測井結(jié)果表明井區(qū)的該地質(zhì)層位的位置為儲層，并且在圖3中 標定為龍崗37的井區(qū)的實際測井結(jié)果表明井區(qū)的該地質(zhì)層位的位置為泥巖層。
圖2示出待測地區(qū)三疊系飛仙關(guān)組內(nèi)部的18Hz的分頻振幅切片中的振幅分布情 況，其中，圖2中標注為210的淺色區(qū)域是強振幅區(qū)，根據(jù)本發(fā)明的預測方法，所述強振幅區(qū) (即，標注為210的淺色區(qū)域)可被確定為分布有儲層的區(qū)域。應該理解，龍崗I處于圖2 中的強振幅區(qū)，這與預測結(jié)果相符。
圖3示出待測地區(qū)三疊系飛仙關(guān)組內(nèi)部的40Hz的分頻振幅切片中的振幅分布情 況，其中，圖3中標注為310的淺色區(qū)域是強振幅區(qū)，根據(jù)本發(fā)明的預測方法，所述強振幅區(qū) (即，標注為310的淺色區(qū)域)可被確定為分布有泥巖層的區(qū)域。應該理解，龍崗37處于圖 3中的強振幅區(qū)，這與預測結(jié)果相符。
通過圖2和圖3的預測和測井結(jié)果可知，使用本發(fā)明的儲層預測方法能夠準確地 預測礁、灘儲層的分布，并可為礁灘勘探開發(fā)提供更好的地震技術(shù)支撐。
盡管已經(jīng)參照其示例性實施例具體顯示和描述了本發(fā)明，但是本領(lǐng)域的技術(shù)人員 應該理解，在不脫離權(quán)利要求所限定的本發(fā)明的精神和范圍的情況下，可以對其進行形式 和細節(jié)上的各種改變。
權(quán)利要求
1.一種預測儲層的方法，所述方法包括(a)采集待測地區(qū)的地震數(shù)據(jù)以及待測地區(qū)的測井數(shù)據(jù)；(b)對測井數(shù)據(jù)進行頻譜分析，以確定對儲層敏感的儲層敏感頻率和對泥巖層敏感的泥巖層敏感頻率；(C)對采集的地震數(shù)據(jù)進行分頻處理，以獲取頻域地震數(shù)據(jù)；(d)從所述頻域地震數(shù)據(jù)獲取待測地區(qū)的特定地質(zhì)層位處的與不同頻率對應的分頻振幅切片；(e)將所述特定地質(zhì)層位的與儲層敏感頻率對應的分頻振幅切片的強振幅區(qū)域的分布確定為所述特定地質(zhì)層位的儲層分布，并將所述特定地質(zhì)層位的與泥巖層敏感頻率對應的分頻振幅切片的強振幅區(qū)域的分布確定為所述特定地質(zhì)層位的泥巖層分布。
2.如權(quán)利要求1所述的方法，還包括在對采集的地震數(shù)據(jù)進行分頻處理之前，對所述地震數(shù)據(jù)進行地震資料處理，其中，在步驟(C)中對經(jīng)過地震資料處理的地震數(shù)據(jù)進行分頻處理來獲取待測地區(qū)的頻域地震數(shù)據(jù)。
3.如權(quán)利要求2所述的方法，其中，所述地震資料處理包括對采集的地震數(shù)據(jù)進行動校正和靜校正中的至少一個。
4.如權(quán)利要求1所述的方法，其中，所述分頻處理的步驟包括對采集的地震數(shù)據(jù)進行離散傅里葉變換或小波變換。
5.如權(quán)利要求1所述的方法，其中，在步驟(b)中對測井數(shù)據(jù)進行頻譜分析的步驟包括對測井數(shù)據(jù)中的聲波測井曲線和電阻率曲線進行廣義S變換頻譜分析。
全文摘要
提供了一種預測儲層的方法，所述方法包括(a)采集待測地區(qū)的地震數(shù)據(jù)以及待測地區(qū)的測井數(shù)據(jù)；(b)對測井數(shù)據(jù)進行頻譜分析，以確定對儲層敏感的儲層敏感頻率和對泥巖層敏感的泥巖層敏感頻率；(c)對采集的地震數(shù)據(jù)進行分頻處理，以獲取頻域地震數(shù)據(jù)；(d)從所述頻域地震數(shù)據(jù)獲取待測地區(qū)的特定地質(zhì)層位處的與不同頻率對應的分頻振幅切片；(e)將所述特定地質(zhì)層位的與儲層敏感頻率對應的分頻振幅切片的強振幅區(qū)域的分布確定為所述特定地質(zhì)層位的儲層分布，并將所述特定地質(zhì)層位的與泥巖層敏感頻率對應的分頻振幅切片的強振幅區(qū)域的分布確定為所述特定地質(zhì)層位的泥巖層分布。
文檔編號G01V1/40GK103048678SQ20121057869
公開日2013年4月17日 申請日期2012年12月27日 優(yōu)先權(quán)日2012年12月27日
發(fā)明者鄧瑛, 徐敏, 王身建, 王萍, 彭才 申請人:中國石油集團川慶鉆探工程有限公司地球物理勘探公司