層和下層砂體的可能厚度，以及上層砂體的起始位置組合，即作為兩層疊置砂體疊置關(guān)系 的最終預測結(jié)果。
[0113] 上述各實施例僅用于說明本發(fā)明，其中各部件的結(jié)構(gòu)、設(shè)置位置及其連接方式等 都是可以有所變化的，凡是在本發(fā)明技術(shù)方案的基礎(chǔ)上進行的等同變換和改進，均不應排 除在本發(fā)明的保護范圍之外。
【主權(quán)項】
1. 一種基于反射系數(shù)分析的砂體疊置關(guān)系判別方法，包括以下步驟： 1) 對目的層進行測井和地震勘探，根據(jù)得到的測井資料和實際地震記錄，估計目的層 砂體和泥巖隔層的厚度范圍；在該厚度范圍內(nèi)，分別設(shè)定兩層疊置砂體中上層砂體、泥巖隔 層和下層砂體的可能厚度，以得到不同的兩層疊置砂體中上層砂體、泥巖隔層和下層砂體 的可能厚度組合;根據(jù)測井資料和實際地震記錄，標定出上層砂體的起始位置； 2) 根據(jù)步驟1)得到的不同的兩層疊置砂體中上層砂體、泥巖隔層和下層砂體的可能厚 度組合，以及上層砂體的起始位置，生成不同的兩層疊置砂體的反射系數(shù)向量； 3) 將步驟2)得到的不同的兩層疊置砂體的反射系數(shù)向量分別與地震子波進行褶積，生 產(chǎn)不同的正演地震記錄； 4) 將實際地震記錄與步驟3)得到的不同的正演地震記錄進行對比，判斷出最接近實際 地震記錄的正演地震記錄，該正演地震記錄對應的兩層疊置砂體的上層砂體、泥巖隔層和 下層砂體可能厚度組合，以及上層砂體起始位置，即為砂體疊置關(guān)系的最終預測結(jié)果。2. 如權(quán)利要求1所述的一種基于反射系數(shù)分析的砂體疊置關(guān)系判別方法，其特征在于， 所述步驟2)中生成不同兩層疊置砂體的反射系數(shù)向量，具體包括以下步驟： ① 根據(jù)兩層疊置砂體中上層砂體、泥巖隔層和下層砂體的可能厚度組合，以及上層砂 體的起始位置，得到兩層疊置砂體中四個分界面的位置； ② 通過測井資料得到目的層砂體和泥巖隔層的阻抗值，根據(jù)兩層疊置砂體的四個分界 面位置，以及砂體和泥巖隔層的阻抗值，生成兩層疊置砂體的反射系數(shù)向量； ③ 根據(jù)不同的兩層疊置砂體中上層砂體、泥巖隔層和下層砂體的可能厚度組合，以及 上層砂體的起始位置，得到不同的分界面位置組合，生成不同的兩層疊置砂體反射系數(shù)向 量。3. 如權(quán)利要求2所述的一種基于反射系數(shù)分析的砂體疊置關(guān)系判別方法，其特征在于， 所述步驟②中生成兩層疊置砂體的反射系數(shù)向量的公式為： r=[r(l)，· · ·，r(i)，· · ·，r(n)]TS= (Il+l2)/(Il_l2) 式中，r表示兩層疊置砂體的反射系數(shù)向量;n表示反射系數(shù)向量r的元素個數(shù);r(i)表 示反射系數(shù)向量r的第i個元素，1 < i <n;du、dM和dD分別表示兩層疊置砂體中上層砂體、泥 巖隔層和下層砂體的可能厚度;do表示上層砂體的起始位置;IjPI 2分別表示砂體和泥巖隔 層的阻抗值;8、￥1、￥2、￥3、；1^1和€2均為計算參數(shù)；當(111 = 0時，￥1 = 0，否則，￥1=1;當(^=0時，￥2 =〇,否貝1J，V2=1 ;當dD = 0時，V3 = 0,否貝lJ，V3=l ;當V2 = V3 = 0時，fl=l，否貝>J，fl = V2;當V2 = VI = 0時，f2 = 1，否貝lj，f2 = V2。4. 如權(quán)利要求1或2或3所述的一種基于反射系數(shù)分析的砂體疊置關(guān)系判別方法，其特 征在于，所述步驟3)中生產(chǎn)不同的正演地震記錄，具體包括以下步驟： I) 從實際地震記錄中提取一個包括nw個元素的列向量作為地震子波，并設(shè)定子波中心 為Cw; II) 分別將不同的兩層疊置砂體的反射系數(shù)向量與地震子波進行褶積，生產(chǎn)不同的正 演地震記錄； III) 將不同的正演地震記錄組合，得到正演地震記錄集合；正演地震記錄集合的每一 列均為一個正演地震記錄。5. 如權(quán)利要求4所述的一種基于反射系數(shù)分析的砂體疊置關(guān)系判別方法，其特征在于， 所述步驟II)中生產(chǎn)正演地震記錄的公式為： y=[y(l)，· · ·，y(i)，· · ·，y(n)]T當 i-j+Cw<l或者i-j+Cw>nw時，w(i-j+Cw) =0 式中，y表示正演地震記錄;η表示正演地震記錄y的元素個數(shù);y (i)表示正演地震記錄y 的第i個元素，1 < i j為計算參數(shù);r(j)表示反射系數(shù)向量r的第j個元素;w(i-j+cw)表 示地震子波w的第i-j+cw個元素。6. 如權(quán)利要求4所述的一種基于反射系數(shù)分析的砂體疊置關(guān)系判別方法，其特征在于， 所述步驟4)中將實際地震記錄與不同的正演地震記錄對比，包括以下步驟： a) 從實際地震記錄中，上下擾動頂層砂體估計位置地多次截取一段包括η個元素的列 向量，作為對照實際地震記錄;η也是正演地震記錄的元素個數(shù)； b) 將正演地震記錄集合的各個列進行歸一化，得到歸一化的正演地震記錄集合； c) 根據(jù)歸一化的正演地震記錄集合和不同的對照實際地震記錄，計算比對向量集合， 比對向量集合的每一列均為一個比對向量； d) 查找比對向量集合中元素的最大值，該元素索引對應的即為最接近實際地震記錄的 正演地震記錄;該正演地震記錄對應的兩層疊置砂體的上層砂體、泥巖隔層和下層砂體的 可能厚度組合，以及上層砂體的起始位置，即為兩層疊置砂體的砂體疊置關(guān)系的最終預測 結(jié)果。7. 如權(quán)利要求6所述的一種基于反射系數(shù)分析的砂體疊置關(guān)系判別方法，其特征在于， 所述步驟b)中正演地震記錄歸一化的公式為： z = [z(l)，...，z(i)，...，z(n)]T式中，z表示歸一化的正演地震記錄;η表示歸一化的正演地震記錄z的元素個數(shù)，也是 正演地震記錄y的元素個數(shù);z(i)表示歸一化的正演地震記錄ζ的第i個元素，1 < i <n;y(i) 表示正演地震記錄y的第i個元素。8. 如權(quán)利要求6或7所述的一種基于反射系數(shù)分析的砂體疊置關(guān)系判別方法，其特征在 于，所述步驟c)中比對向量的計算公式為： Cv - Z Xv ? v - 1，2，. . .，?χ 式中，Ζ表示歸一化的正演地震記錄集合;cv表示與第ν次截取得到的對照實際地震記錄 Xv相對應的比對向量;Xv表示第W欠截取的對照實際地震記錄;nx表示對照實際地震記錄的 截取次數(shù)。9.如權(quán)利要求1或2或3或5或6或7所述的一種基于反射系數(shù)分析的砂體疊置關(guān)系判別 方法，其特征在于，所述步驟1)中提到的厚度均指時間采樣厚度，即地質(zhì)體在時間域地震剖 面上對應的采樣點數(shù)。
【專利摘要】本發(fā)明涉及一種基于反射系數(shù)分析的砂體疊置關(guān)系判別方法，包括以下步驟：1)對目的層進行測井和地震勘探，根據(jù)測井資料和實際地震記錄，設(shè)定兩層疊置砂體中上層砂體、泥巖隔層和下層砂體的可能厚度，并標定上層砂體的起始位置；2)根據(jù)不同的兩層疊置砂體的上層砂體、泥巖隔層和下層砂體可能厚度組合，以及上層砂體起始位置，生成不同的兩層疊置砂體反射系數(shù)向量；3)將不同的兩層疊置砂體反射系數(shù)向量分別與地震子波進行褶積，生產(chǎn)不同的正演地震記錄；4)將實際地震記錄與不同的正演地震記錄進行對比，判斷最接近實際地震記錄的正演地震記錄，該正演地震記錄對應的兩層疊置砂體的砂體疊置關(guān)系，即為最終預測結(jié)果。
【IPC分類】G01V1/30
【公開號】CN105676288
【申請?zhí)枴緾N201610086619
【發(fā)明人】范廷恩, 胡光義, 馬淑芳, 王暉, 宋來明, 王海峰, 王宗俊, 張晶玉, 梁旭, 井涌泉, 肖大坤, 劉向南, 陳飛, 王帥, 解超
【申請人】中國海洋石油總公司, 中海油研究總院
【公開日】2016年6月15日
【申請日】2016年2月16日