本申請涉及油藏地球物理技術(shù)領(lǐng)域，尤其是涉及一種預(yù)測油砂體分布的方法及裝置。

背景技術(shù)：

油砂體是指地下深處含油的砂體，通常由很多不規(guī)則的砂體組成。在油田地質(zhì)工作和油田開發(fā)工作中，如何有效、準(zhǔn)確預(yù)測油砂體和邊界，提高油層鉆遇率，是油藏開發(fā)的關(guān)鍵。

目前比較常用的預(yù)測油砂體分布的方法主要包括兩種，一種是基于測井資料的，井間插值算法，測井資料直觀地反映了儲層的巖性、物性和含油氣性特征及變化，但是在井?dāng)?shù)較少，儲層分別復(fù)雜的油田，利用井間插值預(yù)測油砂體分布可靠性和準(zhǔn)確性都不高；還有一種是基于測井約束下的地震阻抗反演進(jìn)行油砂體預(yù)測，地震資料攜帶了大量的儲層地質(zhì)信息，且在橫向上具有連續(xù)性的優(yōu)勢，但目前基于測井約束下的地震阻抗反演預(yù)測油砂體分布往往只利用了測井聲波和密度信息，且兩種數(shù)據(jù)主要反映儲層的孔隙度參數(shù)，對儲層巖性、含油性的表征能力較弱，影響了油砂體的預(yù)測效果。如何有效利用地震數(shù)據(jù)橫向連續(xù)性優(yōu)勢，進(jìn)一步提高油砂體的預(yù)測精度，一直是業(yè)界的技術(shù)難題。

技術(shù)實現(xiàn)要素：

本申請實施例的目的在于提供一種預(yù)測油砂體分布的方法及裝置，以提高預(yù)測油砂體分布的精度。

為達(dá)到上述目的，本申請實施例提供了一種預(yù)測油砂體分布的方法，所述方法包括：

獲取研究區(qū)內(nèi)的地震數(shù)據(jù)、測井?dāng)?shù)據(jù)和開發(fā)數(shù)據(jù)；

根據(jù)所述地震數(shù)據(jù)、測井?dāng)?shù)據(jù)和開發(fā)數(shù)據(jù)，獲取所述研究區(qū)內(nèi)滿足預(yù)設(shè)條件的樣本砂體，并根據(jù)所述樣本砂體構(gòu)建樣本砂體集合，所述樣本砂體集合中包括至少一個油砂體和至少一個非油砂體；

提取所述樣本砂體集合中各個樣本砂體的特定地震屬性和特定測井屬性，并根據(jù)所述特定地震屬性和特定測井屬性構(gòu)建樣本屬性集合；

根據(jù)所述樣本屬性集合，確定等效權(quán)重矩陣；

獲取所述地震資料中與所述特定地震屬性相對應(yīng)的地震屬性，并結(jié)合所述等效權(quán)重矩陣，預(yù)測所述研究區(qū)內(nèi)的油砂體分布。

本申請實施例還提供了一種預(yù)測油砂體分布的裝置，所述裝置包括：

獲取模塊，用于獲取研究區(qū)內(nèi)的地震數(shù)據(jù)、測井?dāng)?shù)據(jù)和開發(fā)數(shù)據(jù)；

樣本砂體確定模塊，用于根據(jù)所述地震數(shù)據(jù)、測井?dāng)?shù)據(jù)和開發(fā)數(shù)據(jù)，獲取所述研究區(qū)內(nèi)滿足預(yù)設(shè)條件的樣本砂體，并根據(jù)所述樣本砂體構(gòu)建樣本砂體集合，所述樣本砂體集合中包括至少一個油砂體和至少一個非油砂體；

屬性提取模塊，用于提取所述樣本砂體集合中各個樣本砂體的特定地震屬性和特定測井屬性，并根據(jù)所述特定地震屬性和特定測井屬性構(gòu)建樣本屬性集合；

矩陣確定模塊，用于根據(jù)所述樣本屬性集合，確定等效權(quán)重矩陣；

預(yù)測模塊，用于獲取所述地震資料中與所述特定地震屬性相對應(yīng)的地震屬性，并結(jié)合所述等效權(quán)重矩陣，預(yù)測所述研究區(qū)內(nèi)的油砂體分布。

由上述本申請實施例所提供的技術(shù)方案可知，本申請實施例通過研究區(qū)內(nèi)已知的樣本砂體數(shù)據(jù)，確定了可以在識別油砂體過程中代替測井?dāng)?shù)據(jù)作用的等效權(quán)重矩陣，再將該等效權(quán)重矩陣與該研究區(qū)內(nèi)地震數(shù)據(jù)相結(jié)合，就可以預(yù)測研究區(qū)內(nèi)的油砂體分布。本申請實施例所提供技術(shù)方案充分利用了地震數(shù)據(jù)橫向連續(xù)性，可以提高預(yù)測預(yù)測油砂體分布的精度。

附圖說明

此處所說明的附圖用來提供對本申請實施例的進(jìn)一步理解，構(gòu)成本申請實施例的一部分，并不構(gòu)成對本申請實施例的限定。在附圖中：

圖1為本申請實施例的一種預(yù)測油砂體分布的方法流程示意圖；

圖2為本申請實施例的井震標(biāo)定后的過井剖面圖；

圖3為本申請實施例的研究區(qū)具有代表性的單井地質(zhì)沉積柱狀圖；

圖4為本申請實施例的獲取等效權(quán)重矩陣的過程示意圖；

圖5為本申請實施例的研究區(qū)平面分布圖；

圖6為本申請實施例所選擇的地震屬性和測井屬性；

圖7為本申請實施例重復(fù)計算20次獲得的特征權(quán)重數(shù)值；

圖8為本申請實施例20次特征權(quán)重數(shù)值計算的算術(shù)平均值；

圖9為本申請實施例的特征項篩選結(jié)果示意圖；

圖10為本申請實施例的新的敏感屬性集合示意圖；

圖11為本申請實施例的權(quán)系數(shù)計算結(jié)果示意圖；

圖12為本申請實施例的各個敏感屬性項正負(fù)相關(guān)性劃分結(jié)果示意圖；

圖13為本申請實施例的油砂體樣本評價值直方圖；

圖14為本申請實施例的非油砂體樣本評價值直方圖；

圖15為本申請實施例的Y1井的測井和開發(fā)數(shù)據(jù)；

圖16為本申請實施例的Y2井的測井和開發(fā)數(shù)據(jù)；

圖17為本申請實施例利用圖5中所示樣本井中2砂組評價值計算結(jié)果，按照收斂插值算法，所得到的油砂體平面預(yù)測結(jié)果；

圖18為本申請實施例從圖17中提取等值線大于或等于0.35的平面區(qū)域，作為油砂體的平面分布范圍示意圖；

圖19為本申請實施例的利用圖5中所示樣本井和普通井中2砂組評價值計算結(jié)果，按照收斂插值算法，所得到的油砂體平面預(yù)測結(jié)果；

圖20為本申請實施例從圖19中提取等值線大于或等于0.35的平面區(qū)域，作為油砂體的平面分布范圍示意圖；

圖21為本申請實施例的利用圖5中所示樣本井中2砂組的評價值計算結(jié)果，求取等效權(quán)重矩陣，得到油砂體的平面預(yù)測結(jié)果；

圖22為本申請實施例從圖21中提取等值線大于或等于0.35的平面區(qū)域，作為油砂體的平面分布范圍示意圖；

圖23為本申請實施例的一種預(yù)測油砂體分布的裝置示意圖。

具體實施方式

為使本申請實施例的目的、技術(shù)方案和優(yōu)點更加清楚明白，下面結(jié)合實施例和附圖，對本申請實施例做進(jìn)一步詳細(xì)說明。在此，本申請實施例的示意性實施例及其說明用于解釋本申請實施例，但并不作為對本申請實施例的限定。

下面結(jié)合附圖，對本申請實施例的具體實施方式作進(jìn)一步的詳細(xì)說明。

圖1所示為本申請實施例所提供的一種預(yù)測油砂體分布的方法流程示意圖，如圖所示，該方法可以包括以下幾個步驟。該方法可以用于預(yù)測研究區(qū)內(nèi)油砂體的分布。所述油砂體是指地下深處含油的砂體，通常由很多不規(guī)則的砂體組成。

步驟S101，獲取研究區(qū)內(nèi)的地震數(shù)據(jù)、測井?dāng)?shù)據(jù)和開發(fā)數(shù)據(jù)。

在本實施例中，所述測井?dāng)?shù)據(jù)可以包括井徑、自然伽馬、自然電位、聲波時差、密度、中子、深電阻率、中電阻率和淺電阻率等基礎(chǔ)曲線數(shù)據(jù)。所述開發(fā)數(shù)據(jù)可以包括泥質(zhì)含量、孔隙度、滲透率、含水飽和度解釋的成果數(shù)據(jù)和射孔、試油等。所述地震數(shù)據(jù)可以為地震勘探獲得的三維地震數(shù)據(jù)。

步驟S102，根據(jù)所述地震數(shù)據(jù)、測井?dāng)?shù)據(jù)和開發(fā)數(shù)據(jù)，獲取所述研究區(qū)內(nèi)滿足預(yù)設(shè)條件的樣本砂體，并根據(jù)所述樣本砂體構(gòu)建樣本砂體集合，所述樣本砂體集合中包括至少一個油砂體和至少一個非油砂體。

在本實施例中，所述樣本砂體可以為位于所述研究區(qū)內(nèi)，并且已知是否為油砂體的樣本。判斷所述樣本砂體是否為油砂體可以根據(jù)開發(fā)數(shù)據(jù)，通過判斷所述樣本砂體是否含油來實現(xiàn)。在本實施例中按照測井含油氣解釋結(jié)論來判別是否為油砂體，通常認(rèn)為油層、含水油層、油水同層的砂體為油砂體，但在其他實施例中仍要根據(jù)不同區(qū)塊的研究需求更改判別油砂體的標(biāo)準(zhǔn)。另外，由于本實施例所提供的技術(shù)方案的目的在于不依賴井間插值技術(shù)，根據(jù)同時擁有測井?dāng)?shù)據(jù)和地震數(shù)據(jù)的樣本砂體中獲得一個等效權(quán)重矩陣，使得該矩陣在識別是否為油砂體時，可以起到替代井間差值的作用，這樣在其他沒有測井資料，只有地震資料的區(qū)域，就可以利用地震資料和該等效權(quán)重矩陣來預(yù)測油砂體，因此在本實施例中，所述預(yù)設(shè)條件可以包括，樣本砂體同時擁有測井資料和地震資料的砂體。優(yōu)選的，為了方便后續(xù)提取樣本砂體的特征，所述預(yù)設(shè)條件還可以包括，樣本砂體選用沉積環(huán)境穩(wěn)定的砂體。

在本申請的一個實施例中，可以根據(jù)以下幾個步驟確定樣本砂體。

(1)根據(jù)所述地震數(shù)據(jù)和所述測井?dāng)?shù)據(jù)，進(jìn)行井震標(biāo)定，獲得井震標(biāo)定的結(jié)果。

所述井震標(biāo)定可以為使深度域測井分層和時間域地震層位，無論在深度域還是時間域，兩者在空間上均閉合的一種手段。具體的，可以通過速度場將深度域測井分層轉(zhuǎn)換成時間域分層，再結(jié)合時間域地震層位，使其在測井分層和地震測井在時間域上閉合。使得深度域測井分層與時間域地震分層結(jié)果重合的目的之一是為了后續(xù)提取油砂體性質(zhì)時，方便提取對應(yīng)的同一油砂體的地震屬性和測井屬性。

(2)根據(jù)所述地震數(shù)據(jù)、所述測井?dāng)?shù)據(jù)和所述井震標(biāo)定的結(jié)果，在所述研究區(qū)中確定滿足第一預(yù)設(shè)條件的目標(biāo)層。

在本實施例中，所述第一預(yù)設(shè)條件可以為地層沉積環(huán)境穩(wěn)定。沉積環(huán)境穩(wěn)定可以指各層段內(nèi)的巖性特征和沉積相特征基本一致，可以通過測井曲線和地震數(shù)據(jù)獲得。具體的，可以通過測井(自然伽馬、自然電位)巖性曲線的形態(tài)特征，巖心數(shù)據(jù)的顏色、巖性、粒度、沉積構(gòu)造特征，地震數(shù)據(jù)的反射特征，以及前人對研究區(qū)沉積特征的認(rèn)識來綜合判斷巖性特征、沉積相特征是否基本一致。例如，可以量化自然伽馬和/或自然電位曲線的抖動，當(dāng)?shù)貙又卸秳映^一定閾值，可以認(rèn)為相應(yīng)地層沉積環(huán)境不穩(wěn)定。

在本申請的一個具體實施例中，井震標(biāo)定后的過井剖面圖如圖2所示，圖中W4和W32為井編號。根據(jù)測井?dāng)?shù)據(jù)的分層結(jié)果依次為圖2中所示的1_Top、1_Bot、2_Bot、3_Bot、4_Bot和5_Bot，與測井分層相對應(yīng)的地震層位分別為1砂組頂、2砂組頂、3砂組頂、4砂組頂、5砂組頂和6砂組頂。如圖2中所示，在井震標(biāo)定的情況下，深度域的測井分層和時間域的地震層位在空間上閉合，便于分別在時間域提取地震屬性和深度域提取測井屬性，且能保證兩種屬性反映的是同一地層。確定研究目的層時還需要根據(jù)測井和地震數(shù)據(jù)，確定沉積環(huán)境穩(wěn)定的地層。具體的，可以通過判斷地層內(nèi)巖性特性、沉積相特性是否一致判斷沉積環(huán)境是否穩(wěn)定。通常，在研究開始之前，會首先確定研究區(qū)內(nèi)的目的層，例如，在本申請的一個實施例中，研究之前確定研究區(qū)內(nèi)的目的層為深度3000～4000m的地層，經(jīng)過井震標(biāo)定使得深度域測井分層和時間域地震分層在時間域上閉合，再根據(jù)測井和地震數(shù)據(jù)發(fā)現(xiàn)，其中深度3000～3500內(nèi)的地層沉積環(huán)境穩(wěn)定，因此目標(biāo)層為該研究區(qū)內(nèi)深度3000～3500內(nèi)的地層。

在本申請的另一個具體實施例中，研究區(qū)具有代表性的單井地質(zhì)沉積柱狀圖如圖3所示，由于在考慮第一預(yù)設(shè)條件的情況下，只有當(dāng)各個層段內(nèi)的沉積環(huán)境穩(wěn)定即巖性特征和沉積特征基本一致的情況下，才可以將各個層段作為樣本數(shù)據(jù)的來源，否則只能在某一地地層段內(nèi)取樣。由圖3可知，圖3中5個砂組沉積環(huán)境都相對穩(wěn)定，因此，這個5個砂組所在地層可以作為研究目的層。

(3)根據(jù)所述測井資料、所述目標(biāo)層和所述開發(fā)資料，確定滿足第二預(yù)設(shè)條件的樣本井。

在本實施例中，所述第二預(yù)設(shè)條件可以為，樣本井在研究區(qū)內(nèi)盡可能的均勻分布，且樣本井需具有較完整的井徑、自然伽馬、自然電位、聲波時差、密度、中子、深電阻率、中電阻率、淺電阻率等基礎(chǔ)曲線數(shù)據(jù)，泥質(zhì)含量、孔隙度、滲透率、含水飽和度解釋的成果數(shù)據(jù)和射孔、試油等開發(fā)數(shù)據(jù)，能明顯區(qū)分開儲層是否為油砂體。所述第二預(yù)設(shè)條件還可以包括，樣本井需盡可能多的包括目標(biāo)層。例如，步驟(2)中確定的目標(biāo)層為該研究區(qū)內(nèi)深度3000～3500內(nèi)的地層，那么樣本井中要盡可能的包含深度3000～3500內(nèi)的地層，因為只有這樣才方便確定樣本砂體。

(4)從所述樣本井中獲取樣本砂體，并根據(jù)所述樣本砂體構(gòu)建樣本砂體集合。

在本實施例中，由于所述樣本砂體是從樣本井中獲取的，所以可以認(rèn)為所述樣本砂體具有測井?dāng)?shù)據(jù)，又由于研究區(qū)的地震數(shù)據(jù)覆蓋了整個研究區(qū)，所以該樣本砂體還同時擁有地震數(shù)據(jù)。所述樣本砂體可以為所述樣本井內(nèi)的目標(biāo)層中的砂體。

由于目前提取屬性特征對地層厚度的限制，在本申請的各個實施例中，可以將一個砂組近似看成是一個砂體，具體的，目標(biāo)層為研究區(qū)內(nèi)深度3000～3500內(nèi)的地層，某樣本井A中包含深度3000～3500的地層，根據(jù)井震標(biāo)定的結(jié)果，該深度3000～3500的地層中包含5組砂組，因此，在本實施例中，可以將這5組砂組認(rèn)為是樣本砂體。

上述步驟(1)至(4)提供了一種獲取樣本砂體的方法，根據(jù)該方法提取的樣本砂體既同時擁有地震數(shù)據(jù)和測井?dāng)?shù)據(jù)，又具有沉積環(huán)境穩(wěn)定的有點，方便后續(xù)研究。

步驟S103，提取所述樣本砂體集合中各個樣本砂體的特定地震屬性和特定測井屬性，并根據(jù)所述特定地震屬性和特定測井屬性構(gòu)建樣本屬性集合。

所述特定地震屬性和特定測井屬性均可以為若干個。具體提取的特定地震屬性的數(shù)量和特定測井屬性的數(shù)量根據(jù)實際確定。優(yōu)選的，為了盡可能準(zhǔn)確的預(yù)測研究區(qū)油砂體分布，地震屬性和測井屬性的數(shù)量可以盡可能的多一點。在本實施例中，提取所述樣本砂體的特定地震屬性，可以為以地震分層的頂、低分層為約束分別提取若干地震屬性數(shù)據(jù)。其中，提取地震屬性可以包括以下幾種方式：一是輸入整個地震數(shù)據(jù)體，提取體屬性，然后以各層段的頂部層位數(shù)據(jù)作為約束，向下開時窗(時窗的大小要根據(jù)層段的厚度選擇)提取體屬性相應(yīng)的層位屬性數(shù)據(jù)，或是以各層段的頂、底層位數(shù)據(jù)作為約束，提取體屬性的層間屬性數(shù)據(jù)；二是輸入整個地震數(shù)據(jù)體，直接以各層段的頂、底層位數(shù)據(jù)作為約束，提取層間屬性數(shù)據(jù)。通常，提取特定地震屬性時，優(yōu)先選擇體屬性。

提取所述測井屬性可以為，根據(jù)井震標(biāo)定結(jié)果，以與地震分層相對應(yīng)的測井層位數(shù)據(jù)作為各層段的深度域分界線，以各層段的頂、底分層作為約束分別提取各層段內(nèi)的測井?dāng)?shù)據(jù)，例如泥質(zhì)含量、孔隙度、滲透率、含水飽和度、砂體厚度等測井屬性的算術(shù)平均值和幾何平均值等。提取的特定地震屬性和特定測井屬性構(gòu)成了樣本屬性集合。

在本申請的一個實施例中，共有樣本砂體20個，這20個樣本砂體中既包含油砂體又包含非油砂體，對每一個樣本砂體的提取特定的10個特定地震屬性的數(shù)據(jù)，同時提取每一個樣本砂體的特定測井?dāng)?shù)據(jù)10個。這樣每一個樣本砂體就對應(yīng)了20個屬性數(shù)據(jù)，其中包括10個地震屬性的，10個測井屬性的，樣本屬性集合中共包含了20個樣本砂體的共400個屬性數(shù)據(jù)。

步驟S104，根據(jù)所述樣本屬性集合，確定等效權(quán)重矩陣。

根據(jù)樣本屬性集合，獲得一個等效權(quán)重矩陣，可以為獲取一個等效權(quán)重矩陣，使得該矩陣在識別是否為油砂體時，可以起到替代井間差值的作用，這樣在其他沒有測井資料，只有地震資料的區(qū)域，就可以利用地震資料和該等效權(quán)重矩陣來預(yù)測油砂體。

在本申請的一個實施例中，S104具體實施時，可以通過以下步驟確定等效權(quán)重矩陣。

(1)根據(jù)所述樣本屬性集合中的各個屬性，獲得滿足第三預(yù)設(shè)條件的敏感屬性，并根據(jù)所述敏感屬性構(gòu)建敏感屬性集合，所述敏感屬性集合中包括至少一個地震敏感屬性和至少一個測井敏感屬性。

在本實施例中，由于不同屬性對識別油砂體的敏感性和敏感程度有差別，因此可以將樣本屬性集合中的屬性看成是樣本砂體的特征，并通過一定特征篩選算法將樣本屬性集合中的屬性進(jìn)行篩選，篩選的目的可以為留下樣本屬性集合中對識別油砂體敏感的一些屬性。所述屬性之所以對識別油砂體敏感主要是由于該屬性的數(shù)值，對區(qū)分砂體是否為油砂體是有益的或者敏感的。所述第三預(yù)設(shè)條件可以為對識別油砂體敏感。篩選特征的方法可以包括：ReliefF算法、fisher score特征篩選算法、子集搜索算法等。在本實施例中，特征篩選時還需保證篩選得到的敏感屬性集合中包括至少一個地震敏感屬性和至少一個測井敏感屬性。具體的，在本申請的一個實施例中，樣本屬性集合中包含10個地震屬性，10個測井屬性，經(jīng)過特征篩選之后，得到的敏感屬性集合中包含5個地震敏感屬性和5個測井敏感屬性。

(2)根據(jù)所述敏感屬性集合，確定所述等效權(quán)重矩陣。

在本實施例中，敏感屬性集合中保留的屬性，都是對區(qū)別砂體是否為油砂體敏感的一些特征?？梢愿鶕?jù)敏感屬性集合中各個測井敏感屬性和地震敏感屬性的數(shù)值，確定一個識別油砂體的識別標(biāo)準(zhǔn)，再根據(jù)該識別標(biāo)準(zhǔn)反推出當(dāng)只有地震數(shù)據(jù)時，可以替代井間差值作用的等效權(quán)重矩陣。

在本申請的一個實施例中，可以通過兩個步驟，根據(jù)所述敏感屬性集合確定所述等效權(quán)重矩陣。

1、根據(jù)每一個所述樣本砂體的地震敏感屬性和測井敏感屬性，確定與每一個所述樣本砂體相對應(yīng)的評價值，所有所述樣本砂體的評價值構(gòu)成評價值集合。

每一個樣本砂體的地震敏感屬性和測井敏感屬性可以都不止一個，此時，想要通過這個敏感屬性去識別油砂體，需要根據(jù)這個敏感屬性確定與樣本砂體對應(yīng)的一個評價值，通過該評價值來識別油砂體。具體的，在某一個實施例中，樣本砂體共有10個，每一個樣本砂體有5個地震敏感屬性和5個測井敏感屬性，可以將5個地震敏感屬性和5個測井敏感屬性的加權(quán)和定義為評價值，根據(jù)評價值數(shù)值高低判斷識別為油砂體。

2、根據(jù)所述評價值集合和由每一個所述樣本砂體的地震敏感屬性構(gòu)成的集合，確定所述等效權(quán)重矩陣。

上述步驟2的實質(zhì)為確定一個特征權(quán)重矩陣，使得根據(jù)該特征權(quán)重矩陣和每一個所述樣本砂體的地震敏感屬性，獲得與根據(jù)地震敏感屬性和測井敏感屬性得到的評價值矩陣數(shù)值相同矩陣。該過程可以用如圖4所示的示意圖表示。圖4中，a1…表示地震敏感屬性，b1…表示測井敏感屬性。

步驟S105，獲取所述地震資料中與所述特定地震屬性相對應(yīng)的地震屬性，并結(jié)合所述等效權(quán)重矩陣，預(yù)測所述研究區(qū)內(nèi)的油砂體分布。

在本實施例中，若等效權(quán)重矩陣是根據(jù)樣本屬性集合中所有的屬性獲得的，此時在預(yù)測油砂體分布時，獲取所述地震屬性則應(yīng)該與樣本屬性集合中的特定地震屬性相對應(yīng)。具體的，例如樣本屬性集合中包括的地震屬性有a1，a2，a3，a4，此時需要從地震資料中獲取的地震屬性也為a1，a2，a3，a4對應(yīng)的數(shù)值。本實施例中，由于需要預(yù)測研究區(qū)內(nèi)的油砂體分布，所以根據(jù)地震資料獲得的地震屬性是面屬性數(shù)值，是研究區(qū)內(nèi)的連續(xù)變量。在其他實施例中，還可以通過獲取待識別砂體的特定地震屬性，根據(jù)上述方法識別待識別砂體。在另一些實施例中，從地震資料中獲得地震屬性可以為研究區(qū)內(nèi)某一砂組的地震屬性，此時可以根據(jù)上述方法預(yù)測研究區(qū)內(nèi)某一砂組內(nèi)的油砂體分布。

在本申請的一個實施例中，在獲得樣本屬性集合后，對樣本屬性集合中的屬性進(jìn)行了特征篩選，得到了敏感屬性集合，此時S105具體實施時，可以通過以下步驟實現(xiàn)。

(1)獲取所述地震資料中的第一地震敏感屬性，并根據(jù)所述第一地震敏感屬性構(gòu)建第一地震敏感屬性矩陣，所述第一地震敏感屬性與所述敏感屬性集合中的地震敏感屬性相對應(yīng)。

獲取所述地震資料中的第一地震敏感屬性可以為獲取地震資料中的部分地震屬性，該部分地震屬性與敏感屬性集合中包含的地震敏感屬性相對應(yīng)。獲得的第一地震敏感屬性可以指根據(jù)地震資料獲得到的整個研究區(qū)內(nèi)的屬性連續(xù)數(shù)值，用來預(yù)測整個研究區(qū)內(nèi)的油砂體分布。

(2)將所述第一地震敏感屬性矩陣與所述等效權(quán)重矩陣相乘，獲得研究區(qū)內(nèi)評價值分布。

等效權(quán)重矩陣在獲得評價值時，可以替代井間差值的作用，因此將該等效權(quán)重矩陣與第一地震敏感屬性矩陣相乘，就可以得到整個研究區(qū)內(nèi)的評價值分布。

(3)根據(jù)所述樣本砂體的評價值集合，確定評價值識別閾值。

樣本砂體除了可以獲得等效權(quán)重矩陣之外，還可以根據(jù)各個樣本砂體是否為油砂體與評價值數(shù)值關(guān)系，確定油砂體的評價值識別閾值。

(4)根據(jù)所述研究區(qū)內(nèi)評價值分布和所述評價值識別閾值，預(yù)測所述研究區(qū)的油砂體分布。

由圖1的流程圖可知，本申請實施例通過研究區(qū)內(nèi)已知的樣本砂體數(shù)據(jù)，確定了可以在識別油砂體過程中代替測井?dāng)?shù)據(jù)作用的等效權(quán)重矩陣，再將該等效權(quán)重矩陣與該研究區(qū)內(nèi)地震數(shù)據(jù)相結(jié)合，就可以預(yù)測研究區(qū)內(nèi)的油砂體分布。本申請實施例所提供技術(shù)方案充分利用了地震數(shù)據(jù)橫向連續(xù)性，可以提高預(yù)測預(yù)測油砂體分布的精度。

在本申請的一個實施例中，某研究區(qū)的地震數(shù)據(jù)、測井?dāng)?shù)據(jù)和開發(fā)數(shù)據(jù)已知。可以根據(jù)以下步驟來預(yù)測研究區(qū)內(nèi)油砂體。

1、根據(jù)所述地震數(shù)據(jù)和所述測井?dāng)?shù)據(jù)，進(jìn)行井震標(biāo)定。

2、根據(jù)所述地震數(shù)據(jù)、所述測井?dāng)?shù)據(jù)和所述井震標(biāo)定的分層結(jié)果從研究區(qū)目的層中確定目標(biāo)層。

所述目標(biāo)層可以為目的層中沉積環(huán)境穩(wěn)定的地層范圍。在本實施例中，目標(biāo)層中包含了井震標(biāo)定后的砂組1～砂組5(如圖3所示)，由圖3可知，該五組砂組中沉積環(huán)境均穩(wěn)定，因此，該五組砂組所在的地層為目標(biāo)層。

3、根據(jù)所述測井資料、所述目標(biāo)層和所述開發(fā)資料，確定樣本井。

確定樣本井時要滿足幾個條件：樣本井在研究區(qū)內(nèi)要盡可能分布均勻、樣本井要包含前面確定的目標(biāo)層中的全部或者部分范圍，并且樣本井要具有較完整的測井?dāng)?shù)據(jù)和開發(fā)數(shù)據(jù)。在本實施例中，根據(jù)以上幾個要求最終在研究區(qū)內(nèi)確定了9口樣本井，如圖5中W1～W9所示。

4、從所述樣本井中獲取樣本砂體，所述樣本砂體構(gòu)成了樣本砂體集合。

在本實施例中，樣本井W1～W9中的每一口井中均包含了如圖3所示的5組砂組中的至少一個砂組，實際最終獲得了共42個目標(biāo)研究砂組。由于目前提取特征的極限厚度只能達(dá)到砂組的厚度，因此本實施例中砂組等同于砂體，即只要砂組中含有油砂體，該砂組就可以被稱為油砂體。在本實施例中，將油層、含水油層和油水同層的砂體稱為油砂體，根據(jù)這一標(biāo)準(zhǔn)，結(jié)合開發(fā)資料，對本實施例中的42個樣本砂體進(jìn)行分類，其中包括23個油砂體，19個非油砂體。

5、提取所述樣本砂體的特定地震屬性和特定測井屬性，所述特定地震屬性和特定所述測井屬性構(gòu)成了樣本屬性集合。

每一個砂體對應(yīng)的地震屬性和測井屬性都很多，在本實施例中，根據(jù)需要確定了10個地震屬性和8個測井屬性，如下圖6所示。其中，地震屬性分別為：2ms和4ms時窗的均方根振幅、瞬時頻率、瞬時相位、主頻率屬性(域頻率)，以及層間均方根振幅和層間平均振幅。測井屬性分別為：泥質(zhì)含量、滲透率、含水飽和度的算術(shù)平均值和幾何平均值，孔隙度的算法平均值以及砂體厚度值。在本實施例中，所述樣本屬性集合中共有42個樣本砂體(樣本個數(shù))，每一個樣本砂體對應(yīng)了18個屬性數(shù)值。因此，可以將所述樣本屬性集合表示成一個42行18列的矩陣。

6、使用ReliefF特征篩選算法篩選樣本屬性集合中的屬性，獲得對識別油砂體敏感的敏感屬性，所述敏感屬性構(gòu)成了敏感屬性集合，所述敏感屬性集合中包括至少一個地震敏感屬性和至少一個測井敏感屬性。

在本實施例中，運用ReliefF特征篩選算法提取樣本屬性集合中對識別油砂體敏感的敏感屬性。所述的ReliefF算法是一種基于特征對近距離樣本的區(qū)分能力的特征選擇算法，是一種有效的降維方式，即從訓(xùn)練集中隨機選取一個樣本R，其核心要求是特征集包含的特征與類標(biāo)簽高度相關(guān)。利用ReliefF特征選擇算法計算樣本屬性集合中各個屬性的特征權(quán)重，并重復(fù)計算多次(次數(shù)可自由選定，通常10次以上)。在本實施例中，共重復(fù)計算了20次，每一次計算得到一組特征權(quán)重數(shù)值，將20次計算結(jié)果繪制在一個坐標(biāo)中如圖7所示，圖中特征編碼對應(yīng)圖6中各個特征的編碼。將每一個屬性20次計算的特征權(quán)重值平均，可以得到如圖8所示的結(jié)果(圖中特征編碼對應(yīng)圖6中各個特征的編碼)。設(shè)置特征門檻值，篩選特征，將屬性的特征權(quán)重值平均值大于門檻值的屬性留下，作為識別油砂體的敏感屬性。通常門檻值的選取應(yīng)使篩選的敏感屬性集合滿足三個特征：可區(qū)分性好，各屬性應(yīng)具有明顯區(qū)分是否為油砂體；可靠性好，對同類的樣本數(shù)據(jù)來說，特征權(quán)重值應(yīng)該比較接近；數(shù)量少，要含有盡可能少的屬性項個數(shù)。在本實施例中，設(shè)置的門檻值為0.038，圖7中特征權(quán)重值平均大于0.038的屬性被留下，即泥質(zhì)含量算數(shù)平均值、孔隙度算數(shù)平均值、泥質(zhì)含量幾何平均值、2ms的均方根振幅、2ms主頻率屬性(域頻率)、4ms的均方根振幅和4ms主頻率屬性(域頻率)，共計7種敏感屬性構(gòu)成，如圖9所示，圖9中還列出了各個敏感屬性(即篩選特征項)的特征權(quán)重值平均值。本實施例中，敏感屬性集合中包含42個樣本砂體的共計42×7個敏感屬性，該敏感屬性集合可以表示成42行7列的敏感屬性矩陣。

7、使用灰度關(guān)聯(lián)分析方法，求取敏感屬性集合中各個屬性的權(quán)系數(shù)。

所述的灰度關(guān)聯(lián)分析是對動態(tài)過程做發(fā)展趨勢量化分析，實質(zhì)是通過確定參考序列和若干比較序列的相似程度，判斷其是否緊密，反映了數(shù)據(jù)間的關(guān)聯(lián)程度，關(guān)聯(lián)度越大，說明數(shù)據(jù)間越相似。使用灰度關(guān)聯(lián)分析方法的具體過程中，將敏感屬性集合中特征權(quán)重值最大的屬性項作為參考序列，并作為敏感屬性集合中的第一個元素，其余屬性項依次后移，組成新的敏感屬性集合，如圖10所示，圖中編碼對應(yīng)圖6中各個特征的編碼。按照灰度關(guān)聯(lián)分析方法，對新的敏感屬性集合進(jìn)行權(quán)系數(shù)計算，第k個敏感屬性的權(quán)系數(shù)記為γ_k(k＝0，1，……，7)，權(quán)系數(shù)計算結(jié)果如圖11所示，圖中編碼對應(yīng)圖6中各個特征的編碼。

8、將敏感屬性集合中各個敏感屬性的數(shù)值進(jìn)行最大值標(biāo)準(zhǔn)化，并計算各個樣本砂體的評價值。

本實施例中，進(jìn)行最大值標(biāo)準(zhǔn)化時，要區(qū)分各個敏感屬性與油砂體的正負(fù)相關(guān)性。所述的正負(fù)相關(guān)性，正相關(guān)是指某屬性的值越大，表示油砂體的可能性越大；負(fù)相關(guān)是指某屬性的值越小，表示油砂體的可能性越大。根據(jù)行業(yè)內(nèi)對本實施例所涉及的7個敏感屬性的一般認(rèn)識，將這7個屬性進(jìn)行正負(fù)相關(guān)性劃分結(jié)果如圖12所示(圖中編碼對應(yīng)圖6中各個特征的編碼)，圖中“+”表示正相關(guān)，“-”表示負(fù)相關(guān)，因此可知，孔隙度、均方根振幅、主頻率為正相關(guān)，泥質(zhì)含量為負(fù)相關(guān)。對某一敏感屬性數(shù)值最大值標(biāo)準(zhǔn)化時，若該屬性與油砂體正相關(guān)性，則利用各個樣本砂體本屬性的值除以所有樣本砂體本屬性值的最大值；若該屬性與油砂體負(fù)相關(guān)性，則利用所有樣本砂體本屬性值的最大值減去樣本砂體本屬性的值，然后再除以所有樣本砂體本屬性值的最大值，得到最大值標(biāo)準(zhǔn)化后的樣本數(shù)據(jù)集，例如，敏感屬性集合可以表示成20行10列的矩陣，每一行對應(yīng)一個樣本砂體，每一列對應(yīng)一種敏感屬性，最大值標(biāo)準(zhǔn)化時，若某一列對應(yīng)的敏感屬性與油砂體的正相關(guān)，用該列中每一個數(shù)值除以該列的最大值，得到該列最大值標(biāo)準(zhǔn)化的結(jié)果；若某一列對應(yīng)敏感屬性與油砂體負(fù)相關(guān)，將該列中某一個數(shù)值a進(jìn)行最大值標(biāo)準(zhǔn)化計算時，先用該列的最大值減去待標(biāo)準(zhǔn)化數(shù)值a，然后再除以該列的最大值。本實施例中，將第k個敏感屬性的第i個元素標(biāo)最大值準(zhǔn)化結(jié)果記為D_i(k)(i＝0,1，……，42；k＝0，1，……，7)。

得到各個樣本砂體敏感屬性最大值標(biāo)準(zhǔn)化數(shù)值D_i(k)之后，再按照以下公式，計算各個樣本(即各個樣本砂體)的評價值，記E_i為第i個樣本的評價值(i＝0,1，……，42)。

上述公式計算評價值的公式中，將權(quán)系數(shù)作為加權(quán)值，針對每一個樣本砂體都得到了一個綜合敏感屬性集合中各個敏感屬性的綜合評價值，具體的，例如，在本實施例中，某一個樣本砂體A，其敏感屬性集合中最大值標(biāo)準(zhǔn)化后的泥質(zhì)含量算數(shù)平均值為m₁、最大值標(biāo)準(zhǔn)化后的孔隙度算數(shù)平均值為m₂、最大值標(biāo)準(zhǔn)化后的泥質(zhì)含量幾何平均值為m₃、最大值標(biāo)準(zhǔn)化后的2ms的均方根振幅為m₄、最大值標(biāo)準(zhǔn)化后的2ms主頻率屬性(域頻率)為m₅、最大值標(biāo)準(zhǔn)化后的4ms的均方根振幅為m₆以及最大值標(biāo)準(zhǔn)化后的4ms主頻率屬性(域頻率)為m₇，且泥質(zhì)含量算數(shù)平均值、孔隙度算數(shù)平均值、泥質(zhì)含量幾何平均值、2ms的均方根振幅、2ms主頻率屬性(域頻率)、4ms的均方根振幅和4ms主頻率屬性(域頻率)最大值標(biāo)準(zhǔn)化后的權(quán)系數(shù)分別為D(1)、D(2)…、D(7)。則該樣本砂體A的評價值為D(1)·m₁+D(2)·m₂+…D(7)·m₇。

將本實施例中42個樣本砂體樣本的評價值結(jié)果，進(jìn)行直方圖分析，其中，23個油砂體的直方圖如圖13所示，19個非油砂體的評價值直方圖如圖14所示。根據(jù)油砂體和非油砂體評價值的區(qū)別，確定一個油砂體的最小置信區(qū)間。所述最小置信區(qū)間是區(qū)分是否為油砂體的門檻值，在不小于該門檻值的情況下，認(rèn)為預(yù)測的結(jié)果為油砂體，否則為非油砂體。在本實施例中，將油砂體樣本評價值的最小值和非油砂體樣本評價值的最大值取平均值，將該平均值作為油砂體的最小置信區(qū)間，即0.35。

前面8個步驟根據(jù)已知是否為油砂體的樣本砂體樣本，確定了可以用來識別油砂體的評價值，并確定了評價值識別油砂體的最小置信區(qū)間，后續(xù)識別油砂體時，只要將根據(jù)研究區(qū)內(nèi)各個敏感屬性，加權(quán)計算得到評價值，就可以根據(jù)最小置信區(qū)間判斷是否為油砂體，但是研究區(qū)內(nèi)測井?dāng)?shù)據(jù)并不是所有區(qū)域都有的，因此針對研究區(qū)內(nèi)那些沒有測井?dāng)?shù)據(jù)的區(qū)域，若想要利用評價值的方法識別油砂體，則需要借助其他一些參數(shù)或者方法來識別。

9、根據(jù)各個樣本砂體的評價值，獲得等效權(quán)重矩陣。

等效權(quán)重矩陣可以使得在沒有測井?dāng)?shù)據(jù)的數(shù)據(jù)，同樣可以獲得砂體的評價值。在本實施例中，獲得等效權(quán)重矩陣的方法可以如下步驟所示。

(1)將各個樣本砂體敏感屬性集合中的地震敏感屬性，記為矩陣A，具體的，在本實施例中，A為一個42行3列的矩陣，記為A_42×3，其中，矩陣A的每一行對應(yīng)一個樣本砂體，每一列對應(yīng)一種地震敏感屬性。

(2)在矩陣A的右側(cè)加上一個列向量，該列向量的各個元素均為1，得到一個新的矩陣，記為B_42×4。

(3)將步驟8得到的各個樣本砂體的評價值，按照與矩陣B(或A)中各行相對應(yīng)的順序，組成評價值列向量E_42×1，并求解一個等效權(quán)重矩陣C_4×1，滿足B_42×4·C_4×1＝E_42×1。

矩陣C_4×1在后續(xù)該研究區(qū)內(nèi)其他沒有測井?dāng)?shù)據(jù)的區(qū)域中，就可以取代測井?dāng)?shù)據(jù)的作用。

10、根據(jù)步驟9中得到的等效權(quán)重矩陣，預(yù)測研究區(qū)內(nèi)的油砂體。

具體的，在本實施例中，將研究區(qū)地震資料中的敏感地震屬性數(shù)值按照在矩陣B(或A)中各個列的排列順序，組成矩陣B_n×4。將該矩陣B_n×4與等效權(quán)重矩陣C_4×1相乘，得到研究區(qū)內(nèi)評價值。再以0.35為界限，判斷研究區(qū)內(nèi)油砂體分布。具體的，若待識別砂體的評價值達(dá)到0.35，則為油砂體，否則為非油砂體。

本實施例中，本申請實施例通過研究區(qū)內(nèi)已知的樣本砂體數(shù)據(jù)，確定了可以在識別油砂體過程中代替測井?dāng)?shù)據(jù)作用的等效權(quán)重矩陣，再將該等效權(quán)重矩陣與該研究區(qū)內(nèi)地震數(shù)據(jù)相結(jié)合，就可以預(yù)測研究區(qū)內(nèi)的油砂體分布。本申請實施例所提供技術(shù)方案充分利用了地震數(shù)據(jù)橫向連續(xù)性，可以提高預(yù)測預(yù)測油砂體分布的精度。

為了進(jìn)一步驗證本實施例所識別油砂體的準(zhǔn)確性，在本申請的一個實施例中，還在上一實施例所在研究區(qū)進(jìn)行了驗證試驗。在該研究區(qū)內(nèi)結(jié)合現(xiàn)有的收斂差值算法進(jìn)行油砂體預(yù)測，并對比按照上述實施例步驟1至10所預(yù)測的結(jié)果，以此來驗證本申請實施例所提方法的識別準(zhǔn)確性，以如圖5中的Y1和Y2作為驗證井，并以如圖3中所示2砂組為例。根據(jù)驗證井Y1和Y2的測井和開發(fā)數(shù)據(jù)(如圖15和圖16所示)，可知這兩口井在2砂組段為油砂體層段。圖17為利用圖5中所示樣本井中2砂組評價值計算結(jié)果，按照收斂插值算法，所得到的油砂體平面預(yù)測結(jié)果，圖18是從圖17中提取等值線大于或等于0.35的平面區(qū)域，作為油砂體的平面分布范圍示意圖。圖19是利用圖5中所示樣本井和普通井中2砂組評價值計算結(jié)果，按照收斂插值算法，所得到的油砂體平面預(yù)測結(jié)果，同樣圖20是圖19中提取等值線大于或等于0.35的平面區(qū)域，作為油砂體的平面分布范圍示意圖。圖21是利用圖5中所示樣本井中2砂組的評價值計算結(jié)果，按照上一實施例中步驟1至步驟10的流程，通過求取等效權(quán)重矩陣，將敏感地震屬性項與等效權(quán)重矩陣進(jìn)行乘積求和運算，得到油砂體的平面預(yù)測結(jié)果；圖22是圖21中提取等值線大于或等于0.35的平面區(qū)域，作為油砂體的平面分布范圍示意圖。

在圖18和圖20中，驗證井Y1和Y2中的2砂組為非油砂體，與實際不符，但在圖22中，驗證井Y1和Y2中的2砂組為油砂體，與實際一致，因此可以認(rèn)為本申請實施例所提供的方法識別準(zhǔn)確率更高。

另一方面，對比圖20和圖18可知，在井?dāng)?shù)較少的情況，由于空間約束井?dāng)?shù)較少，在插值得到的預(yù)測平面圖中，井間油砂體的分布范圍不易確定，尤其是油砂體邊界更難以準(zhǔn)確預(yù)測；在井?dāng)?shù)較多的情況，隨著空間約束井?dāng)?shù)的增加，在插值得到的預(yù)測平面圖中，油砂體的分布范圍準(zhǔn)確性有所提升，但是井間油砂體的邊界依然容易受到插值算法的影響。而利用本申請實施例所提供的方法預(yù)測油砂體預(yù)測，通過求取等效權(quán)重矩陣，建立了地震平面屬性與井點油砂體評價值的數(shù)學(xué)關(guān)系，由于不受插值算法的影響，即使在井?dāng)?shù)較少的情況下，依然可以很好地預(yù)測井間油砂體的平面分布情況。

本申請實施例中還提供了一種預(yù)測油砂體分布的裝置，如下面的實施例所述。由于該裝置解決問題的原理與一種預(yù)測油砂體分布的方法相似，因此該裝置的實施可以參見一種預(yù)測油砂體分布的方法實施，重復(fù)之處不再贅述。

如圖23所示，本申請實施例所提供的一種預(yù)測油砂體分布的裝置，可以包括以下幾個模塊。

獲取模塊2301，用于獲取研究區(qū)內(nèi)的地震數(shù)據(jù)、測井?dāng)?shù)據(jù)和開發(fā)數(shù)據(jù)。

樣本砂體確定模塊2302，用于根據(jù)所述地震數(shù)據(jù)、測井?dāng)?shù)據(jù)和開發(fā)數(shù)據(jù)，獲取所述研究區(qū)內(nèi)滿足預(yù)設(shè)條件的樣本砂體，構(gòu)成樣本砂體集合，所述樣本砂體集合中包括至少一個油砂體和至少一個非油砂體。

屬性提取模塊2303，用于提取所述樣本砂體集合中各個樣本砂體的特定地震屬性和特定測井屬性，所述特定地震屬性和特定測井屬性構(gòu)成樣本屬性集合。

矩陣確定模塊2304，用于根據(jù)所述樣本屬性集合，確定等效權(quán)重矩陣。

預(yù)測模塊2305，用于獲取所述地震資料中與所述特定地震屬性相對應(yīng)的地震屬性，并結(jié)合所述等效權(quán)重矩陣，預(yù)測所述研究區(qū)內(nèi)的油砂體分布。

由上述裝置的實施例可知，本申請實施例通過研究區(qū)內(nèi)既有地震數(shù)據(jù)又有測井?dāng)?shù)據(jù)，且已知是否為油砂體的樣本砂體數(shù)據(jù)，確定了可以在識別油砂體過程中代替測井?dāng)?shù)據(jù)作用的等效權(quán)重矩陣，再將該等效權(quán)重矩陣與該研究區(qū)內(nèi)地震數(shù)據(jù)相結(jié)合，就可以預(yù)測研究區(qū)內(nèi)的油砂體分布。本申請實施例所提供技術(shù)方案充分利用了地震數(shù)據(jù)橫向連續(xù)性，且設(shè)計地震屬性范圍也更廣，可以提高預(yù)測預(yù)測油砂體分布的精度。

本申請實施例中所描述的方法的步驟可以直接嵌入硬件、處理器執(zhí)行的軟件模塊、或者這兩者的結(jié)合。軟件模塊可以存儲于RAM存儲器、閃存、ROM存儲器、EPROM存儲器、EEPROM存儲器、寄存器、硬盤、可移動磁盤、CD-ROM或本領(lǐng)域中其它任意形式的存儲媒介中。示例性地，存儲媒介可以與處理器連接，以使得處理器可以從存儲媒介中讀取信息，并可以向存儲媒介存寫信息?？蛇x地，存儲媒介還可以集成到處理器中。處理器和存儲媒介可以設(shè)置于ASIC中，ASIC可以設(shè)置于用戶終端中?？蛇x地，處理器和存儲媒介也可以設(shè)置于用戶終端中的不同的部件中。

在一個或多個示例性的設(shè)計中，本申請實施例所描述的上述功能可以在硬件、軟件、固件或這三者的任意組合來實現(xiàn)。如果在軟件中實現(xiàn)，這些功能可以存儲與電腦可讀的媒介上，或以一個或多個指令或代碼形式傳輸于電腦可讀的媒介上。電腦可讀媒介包括電腦存儲媒介和便于使得讓電腦程序從一個地方轉(zhuǎn)移到其它地方的通信媒介。存儲媒介可以是任何通用或特殊電腦可以接入訪問的可用媒體。例如，這樣的電腦可讀媒體可以包括但不限于RAM、ROM、EEPROM、CD-ROM或其它光盤存儲、磁盤存儲或其它磁性存儲裝置，或其它任何可以用于承載或存儲以指令或數(shù)據(jù)結(jié)構(gòu)和其它可被通用或特殊電腦、或通用或特殊處理器讀取形式的程序代碼的媒介。此外，任何連接都可以被適當(dāng)?shù)囟x為電腦可讀媒介，例如，如果軟件是從一個網(wǎng)站站點、服務(wù)器或其它遠(yuǎn)程資源通過一個同軸電纜、光纖電纜、雙絞線、數(shù)字用戶線(DSL)或以例如紅外、無線和微波等無線方式傳輸?shù)囊脖话谒x的電腦可讀媒介中。所述的碟片(disk)和磁盤(disc)包括壓縮磁盤、鐳射盤、光盤、DVD、軟盤和藍(lán)光光盤，磁盤通常以磁性復(fù)制數(shù)據(jù)，而碟片通常以激光進(jìn)行光學(xué)復(fù)制數(shù)據(jù)。上述的組合也可以包含在電腦可讀媒介中。

以上所述的具體實施例，對本申請的目的、技術(shù)方案和有益效果進(jìn)行了進(jìn)一步詳細(xì)說明，所應(yīng)理解的是，以上所述僅為本申請實施例的具體實施例而已，并不用于限定本申請的保護范圍，凡在本申請的精神和原則之內(nèi)，所做的任何修改、等同替換、改進(jìn)等，均應(yīng)包含在本申請的保護范圍之內(nèi)。