本發(fā)明涉及地震資料處理領(lǐng)域，更具體地，涉及一種基于貝葉斯反演框架的反射系數(shù)估計(jì)方法和一種基于貝葉斯反演框架的反射系數(shù)估計(jì)裝置。

背景技術(shù)：

眾所周知，地震記錄通常是帶限的，由于接收和大地非彈性參數(shù)傳播效應(yīng)的影響，地震資料一般都缺少低頻和高頻分量。低頻的缺失會(huì)導(dǎo)致反演的多解性；而高頻的缺失使得資料的分辨率比較低。反褶積效果的好壞對(duì)整個(gè)處理的影響至關(guān)重要，然而不同地區(qū)不同條件下提高分辨率的難度是不同的。大地吸收、噪聲干擾、幾何擴(kuò)散以及地震資料處理中的一系列問(wèn)題都會(huì)給分辨率造成不同的影響。提高分辨率處理的依據(jù)主要有：不同頻率成分有不同的信噪比；不同信噪比的頻率成分對(duì)分辨率的貢獻(xiàn)是不同的。對(duì)于地震資料而言，提高分辨率處理的重點(diǎn)在于高頻段信號(hào)的可靠恢復(fù)，如反q濾波、反褶積、譜白化等方法。多道預(yù)測(cè)反褶積的主要目的是壓制面波等干擾；地表一致性反褶積的主要目的是調(diào)整相位，使有效波一致化，頻率適當(dāng)提高。

傳統(tǒng)的褶積模型把地震記錄看作是地震子波與地下反射系數(shù)序列相褶積的結(jié)果；而反褶積則是通過(guò)壓縮地震子波從地震記錄中反演地下的反射系數(shù)序列。然而，地震反褶積實(shí)際上是一個(gè)“盲”過(guò)程，因?yàn)橥ǔ５乇肀ぐl(fā)的地震子波是未知的，大地濾波效應(yīng)(反射系數(shù)序列)也是未知的。而在地震子波和反射系數(shù)都未知時(shí)，常常要做統(tǒng)計(jì)性假設(shè)，如：地震子波是最小相位的或地層反射系數(shù)是高斯白噪聲等，因此傳統(tǒng)的反褶積也是一種統(tǒng)計(jì)性反褶積。

傳統(tǒng)的反褶積方法通常根據(jù)這些假設(shè)條件，用地震記錄的自相關(guān)代替子波的自相關(guān)，同時(shí)使用基于二階統(tǒng)計(jì)量的線性濾波技術(shù)(如：預(yù)測(cè)反褶積，脈沖反褶積等)來(lái)實(shí)現(xiàn)子波估計(jì)與反褶積。這些假設(shè)和相應(yīng)的方法在實(shí)際應(yīng)用中取得了一定的效果，但是它們往往不符合實(shí)際情況，并且基于二階統(tǒng)計(jì)量的反褶積方法不包含子波的相位信息。實(shí)際地下介質(zhì)中傳播的地震子波往往是混合相位的；地層反射系數(shù)也不是高斯白噪聲。這些問(wèn)題阻礙了反褶積估計(jì)地震反射系數(shù)的精度。

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

本發(fā)明提出了能夠獲得更為準(zhǔn)確的反射系數(shù)估計(jì)的方法和裝置。針對(duì)上述反褶積中沿用的一系列不確切的假設(shè)，從而獲得更為準(zhǔn)確的反射系數(shù)估計(jì)。本發(fā)明還提出了相應(yīng)的裝置。

根據(jù)本發(fā)明的一方面，提出了一種基于貝葉斯反演框架的反射系數(shù)估計(jì)方法，包括：

步驟101，基于下式求解子波序列w＝[w0,…,wp]^t：

其中，r表示反射系數(shù)序列r的褶積矩陣，r＝[r1,r2,...,rm]^t,r0為初始反射系數(shù)序列，m是反射系數(shù)序列r中的元素個(gè)數(shù)，qw是由[μjr(r)+αji(r)]^t構(gòu)成的對(duì)角陣，ρr(ri)表示反射系數(shù)序列r中元素ri的先驗(yàn)概率，c為積分算子矩陣，表示為ξ＝cr，μ和α均為給定的約束因子，d＝[d1,d2,...,dn]^t是觀測(cè)到的一道地震數(shù)據(jù),d＝gr+n,g是n×m維子波褶積矩陣,n＝[n1,n2,...,nn]^t表示觀測(cè)噪聲,進(jìn)入步驟102；

步驟102，基于當(dāng)前得到的w重新計(jì)算反射系數(shù)序列r，進(jìn)入步驟103；

步驟103，基于當(dāng)前得到的w和r判斷盲反褶積目標(biāo)函數(shù)j(w,r)是否收斂，如果收斂則確定當(dāng)前得到的r即為估計(jì)的反射系數(shù)序列；如果不收斂，則將r0更新為當(dāng)前得到的r，返回步驟101，進(jìn)行下一輪迭代，直至盲反褶積目標(biāo)函數(shù)j(w,r)收斂，其中，盲反褶積目標(biāo)函數(shù)j(w,r)為：

j(w,r)＝j(luò)0(w,r)+μjr(r)+αji(r)，

其中，

根據(jù)本發(fā)明的另一方面，提出了一種基于貝葉斯反演框架的反射系數(shù)估計(jì)裝置，包括：子波序列求解單元，用于基于下式求解子波序列w＝[w0,…,wp]^t：

其中，r表示反射系數(shù)序列r的褶積矩陣，r＝[r1,r2,...,rm]^t,r0為初始反射系數(shù)序列，m是反射系數(shù)序列r中的元素個(gè)數(shù)，qw是由[μjr(r)+αji(r)]^t構(gòu)成的對(duì)角陣，ρr(ri)表示反射系數(shù)序列r中元素ri的先驗(yàn)概率，c為積分算子矩陣，表示為ξ＝cr，μ和α均為給定的約束因子，d＝[d1,d2,...,dn]^t是觀測(cè)到的一道地震數(shù)據(jù),d＝gr+n,g是n×m維子波褶積矩陣,n＝[n1,n2,...,nn]^t表示觀測(cè)噪聲,反射系數(shù)求解單元，用于基于子波序列求解單元當(dāng)前得到的w重新計(jì)算反射系數(shù)序列r；控制單元，用于基于子波序列求解單元當(dāng)前得到的w和反射系數(shù)求解單元當(dāng)前得到r判斷盲反褶積目標(biāo)函數(shù)j(w,r)是否收斂，如果收斂則確定當(dāng)前得到的r即為估計(jì)的反射系數(shù)序列；如果不收斂，則將r0更新為當(dāng)前得到的r，返回步驟101，進(jìn)行下一輪迭代，直至盲反褶積目標(biāo)函數(shù)j(w,r)收斂，其中，盲反褶積目標(biāo)函數(shù)j(w,r)為：

j(w,r)＝j(luò)0(w,r)+μjr(r)+αji(r)，

其中，

本發(fā)明的各方面，基于貝葉斯反演框架的盲反褶積，通過(guò)統(tǒng)計(jì)學(xué)原理，減少反褶積算法中的各項(xiàng)假設(shè)，相比于現(xiàn)有技術(shù)中反褶積沿用的一系列不確切的假設(shè)，能顯著提高得到的反射系數(shù)的準(zhǔn)確度。

附圖說(shuō)明

通過(guò)結(jié)合附圖對(duì)本發(fā)明示例性實(shí)施方式進(jìn)行更詳細(xì)的描述，本發(fā)明的上述以及其它目的、特征和優(yōu)勢(shì)將變得更加明顯，其中，在本發(fā)明示例性實(shí)施方式中，相同的參考標(biāo)號(hào)通常代表相同部件。

圖1示出了根據(jù)本發(fā)明的一個(gè)實(shí)施例的基于貝葉斯反演框架的反射系數(shù)估計(jì)方法的流程圖。

圖2示出了中國(guó)西部某工區(qū)的地震成像剖面。

圖3示出了采用傳統(tǒng)反褶積對(duì)圖2所示地震成像剖面進(jìn)行處理得到的結(jié)果。

圖4示出了采用本發(fā)明對(duì)圖2所示地震成像剖面進(jìn)行處理得到的結(jié)果。

圖5示出了圖2、圖3、圖4的頻譜分析對(duì)比示意圖。

具體實(shí)施方式

下面將參照附圖更詳細(xì)地描述本發(fā)明的優(yōu)選實(shí)施方式。雖然附圖中顯示了本發(fā)明的優(yōu)選實(shí)施方式，然而應(yīng)該理解，可以以各種形式實(shí)現(xiàn)本發(fā)明而不應(yīng)被這里闡述的實(shí)施方式所限制。相反，提供這些實(shí)施方式是為了使本發(fā)明更加透徹和完整，并且能夠?qū)⒈景l(fā)明的范圍完整地傳達(dá)給本領(lǐng)域的技術(shù)人員。

此處，先對(duì)本發(fā)明中涉及的原理進(jìn)行簡(jiǎn)單說(shuō)明。

在貝葉斯理論中，后驗(yàn)概率分布函數(shù)ppdf以p(m|d)表示，在給定數(shù)據(jù)向量d下，模型向量m的隨機(jī)分布為：

其中p(d|m)為似然函數(shù)，p(m)為先驗(yàn)隨機(jī)分布函數(shù)。如果僅關(guān)心ppdf的形狀，分母p(d)是一個(gè)可被忽略的常數(shù)。則有：

p(m|d)∝p(m)p(d|m)，公式2

考慮地震褶積模型：

d＝gr+n，公式3

其中d＝[d1,d2,...,dn]^t是觀測(cè)到的一道地震數(shù)據(jù)，n表示地震數(shù)據(jù)的樣點(diǎn)個(gè)數(shù)，r＝[r1,r2,...,rm]^t是反射系數(shù)序列，m是反射系數(shù)序列r中的元素個(gè)數(shù),g是n×m維子波褶積矩陣，n＝[n1,n2,...,nn]^t表示觀測(cè)噪聲。設(shè)子波序列為w＝[w0,…,wp]^t，則有n＝m+p-1。由貝葉斯公式2可得以下近似：

p(r|d)∝p(r)p(d|r)，公式4

式中：p(r|d)表示反射系數(shù)序列r的后驗(yàn)概率，p(r)表示反射系數(shù)序列r的先驗(yàn)概率，p(d|r)表示似然函數(shù)。

假定噪聲向量n服從均值為(通常設(shè))，方差為σn²的高斯分布即：且獨(dú)立，則似然函數(shù)可等價(jià)于噪聲的概率密度分布，表示為：

其中

同樣不妨假設(shè)先驗(yàn)概率依然服從高斯分布即：

其中

把式公式5、公式6代入式公式4，則r的后驗(yàn)概率可表示為：

使上式取最大值的解即為反射系數(shù)序列r的最優(yōu)解，也就是所謂的最大后驗(yàn)概率估計(jì)(maximumaposterioriestimator，map)。

對(duì)公式7兩邊同時(shí)取對(duì)數(shù)：

并略去常數(shù)項(xiàng)-lnk1k2，定義新的目標(biāo)函數(shù)為：

j＝(gr-d)^t(gr-d)+μr^tr公式9

其中顯然，對(duì)公式9的最小化即等價(jià)于對(duì)公式10的最大化。因此可定義目標(biāo)函數(shù)的一般形式為：

其中，ρr(ri)是表示反射系數(shù)序列r中元素ri的先驗(yàn)概率。公式10即為針對(duì)非最小相位子波和非白噪反射系列情況下的盲反褶積數(shù)學(xué)框架，其中j0(w,r)表示數(shù)據(jù)誤差最小二乘擬合，其目的是為了尋求最佳子波和反射系數(shù)使噪聲能量最小；jr(r)表示反射系數(shù)的先驗(yàn)約束項(xiàng)，它體現(xiàn)了反射系數(shù)的隨機(jī)性和統(tǒng)計(jì)學(xué)特征。

可以采用修正的柯西準(zhǔn)則，減小強(qiáng)弱對(duì)比效應(yīng)，使弱反射也能較好的體現(xiàn)出來(lái)，其表達(dá)式為：

其中：

在有井資料、或其他地質(zhì)資料可利用的情況下，通?？蛇\(yùn)用一些簡(jiǎn)單的數(shù)學(xué)變換引入先驗(yàn)信息，來(lái)得到地球物理反演問(wèn)題更適定、更有意義的解。由于波阻抗可看作是反射系數(shù)對(duì)時(shí)間的積分，因此我們把它引入到目標(biāo)函數(shù)中作為反射系數(shù)的先驗(yàn)約束。

由反射系數(shù)的定義可知：

其中，i(ti)表示第i層波阻抗。當(dāng)反射系數(shù)足夠小的情況下(|r(i)|＜0.3)，有以下近似：

上式對(duì)時(shí)間的積分，可以得到相對(duì)波阻抗的表達(dá)式：

i(t0)表示初始波阻抗值，公式14的離散形式表示為：

把公式15寫(xiě)成矩陣形式：

ξ＝cr公式16

其中c為積分算子矩陣，其形式可表示為：

則由最小二乘法可定義波阻抗約束項(xiàng)如下：

最終建立盲反褶積目標(biāo)函數(shù)為：

j(w,r)＝j(luò)0(w,r)+μjr(r)+αji(r)公式19

其中第一項(xiàng)、第二項(xiàng)如公式10中定義；第三項(xiàng)波阻抗約束，用來(lái)控制反演結(jié)果的準(zhǔn)確性和穩(wěn)定性。μ、α分別為稀疏約束和波阻抗約束因子，μ越大，反射系數(shù)越稀疏，α越大，反演結(jié)果越準(zhǔn)確。

根據(jù)上述約束方程，可采用松弛交替迭代法同時(shí)求取反射系數(shù)和地震子波。具體步驟如下：

1)對(duì)任意反射系數(shù)r，目標(biāo)函數(shù)可寫(xiě)為：

將上式寫(xiě)成矩陣形式：

其中r表示反射系數(shù)r的褶積矩陣，qw是由[μjr(r)+αji(r)]^t構(gòu)成的對(duì)角陣。為了進(jìn)一步加快收斂速度，可假設(shè)一個(gè)預(yù)條件方程：

則原矩陣方程式21可改為：

設(shè)則有：

2)、把第1)步中得到的地震子波序列w代入公式19，并求r的偏導(dǎo)：

當(dāng)時(shí)，

即：

把上式寫(xiě)成矩陣形式：

其中矩陣q的對(duì)角元素為：

最終可以得到反射系數(shù)序列r的表達(dá)式為：

r＝(g^tg+μq+αc^tc)^-1(g^td+αc^tξ)公式31

實(shí)施例1

圖1示出了根據(jù)本發(fā)明的一個(gè)實(shí)施例的基于貝葉斯反演框架的反射系數(shù)估計(jì)方法的流程圖。在本實(shí)施例中，該方法包括：

步驟101，基于下式求解子波序列w＝[w0,…,wp]^t：

其中，r表示反射系數(shù)序列r的褶積矩陣，r＝[r1,r2,...,rm]^t,r0為初始反射系數(shù)序列，m是反射系數(shù)序列r中的元素個(gè)數(shù)，qw是由[μjr(r)+αji(r)]^t構(gòu)成的對(duì)角陣，ρr(ri)表示反射系數(shù)序列r中元素ri的先驗(yàn)概率，c為積分算子矩陣，表示為ξ＝cr，μ和α均為給定的約束因子，d＝[d1,d2,...,dn]^t是觀測(cè)到的一道地震數(shù)據(jù),d＝gr+n,g是n×m維子波褶積矩陣,n＝[n1,n2,...,nn]^t表示觀測(cè)噪聲,進(jìn)入步驟102；

步驟102，基于求解得到的w重新計(jì)算反射系數(shù)序列r，進(jìn)入步驟103；

步驟103，基于當(dāng)前的w和r判斷盲反褶積目標(biāo)函數(shù)j(w,r)是否收斂，如果收斂則確定當(dāng)前的r即為估計(jì)得到的反射系數(shù)序列；如果不收斂，則將r0更新為當(dāng)前的r，返回步驟101，進(jìn)行下一輪迭代，直至盲反褶積目標(biāo)函數(shù)j(w,r)收斂，其中，盲反褶積目標(biāo)函數(shù)j(w,r)為：

j(w,r)＝j(luò)0(w,r)+μjr(r)+αji(r)，

其中，

本實(shí)施例基于貝葉斯反演框架的盲反褶積，通過(guò)統(tǒng)計(jì)學(xué)原理，減少反褶積算法中的各項(xiàng)假設(shè)，相比于現(xiàn)有技術(shù)中反褶積沿用的一系列不確切的假設(shè)，能顯著提高得到的反射系數(shù)的準(zhǔn)確度。

在一個(gè)示例中，在上述步驟101中，可以采用共軛梯度迭代求解公式24以得到，迭代終止條件為：

其中，j^k(w,r)表示該共軛梯度迭代中第k次迭代后得到的j(w,r)，ε1為預(yù)設(shè)參數(shù)?？上扔泄曹椞荻鹊?jì)算再運(yùn)用得到地震子波序列w＝[w0,…,wp]^t。

其中，共軛梯度迭代的具體計(jì)算格式可以如下：①定義初始值：r←fx-d,β←0。②迭代：{δx←f′r，s←δx+βs，δr←fs，α←-γ/(δr·δr)，x←x+αs，r←r+αδr}。式中f是過(guò)渡矩陣，x表示模型空間，d表示數(shù)據(jù)空間，r表示誤差，β是調(diào)節(jié)系數(shù)。

在一個(gè)示例中，在步驟102中，可以通過(guò)直接求逆的方法重新計(jì)算反射系數(shù)序列r。具體地，可以基于下式重新計(jì)算反射系數(shù)序列r：

r＝(g^tg+μq+αc^tc)^-1(g^td+αc^tξ)，公式31

其中，q為對(duì)角矩陣，其第i行第i列的元素q(m_cauchy)ii為：

σr為給定值，表示在反射系數(shù)序列r的先驗(yàn)概率分布符合高斯分布的情況下反射系數(shù)序列r的方差。

在另一個(gè)示例中，在步驟102中，參考公式25，可以基于下式重新計(jì)算反射系數(shù)序列r：

在這種情況下，可以采用共軛梯度迭代求解公式32以得到r，迭代終止條件為：

其中，j^k(w,r)表示該共軛梯度迭代中第k次迭代后得到的j(w,r)，ε2為預(yù)設(shè)參數(shù)。

實(shí)施例2

在本發(fā)明的另一實(shí)施例中，公開(kāi)了一種基于貝葉斯反演框架的反射系數(shù)估計(jì)裝置。在本實(shí)施例中，該裝置包括子波序列求解單元、反射系數(shù)求解單元和控制單元。

子波序列求解單元，用于基于下式求解子波序列w＝[w0,…,wp]^t：

其中，r表示反射系數(shù)序列r的褶積矩陣，r＝[r1,r2,...,rm]^t,r0為初始反射系數(shù)序列，m是反射系數(shù)序列r中的元素個(gè)數(shù)，qw是由[μjr(r)+αji(r)]^t構(gòu)成的對(duì)角陣，ρr(ri)表示反射系數(shù)序列r中元素ri的先驗(yàn)概率，c為積分算子矩陣，表示為ξ＝cr，μ和α均為給定的約束因子，d＝[d1,d2,...,dn]^t是觀測(cè)到的一道地震數(shù)據(jù),d＝gr+n,g是n×m維子波褶積矩陣,n＝[n1,n2,...,nn]^t表示觀測(cè)噪聲,

反射系數(shù)求解單元，用于基于子波序列求解單元當(dāng)前得到的w重新計(jì)算反射系數(shù)序列r。

控制單元，用于基于子波序列求解單元當(dāng)前得到的w和反射系數(shù)求解單元當(dāng)前得到r判斷盲反褶積目標(biāo)函數(shù)j(w,r)是否收斂，如果收斂則確定當(dāng)前得到的r即為估計(jì)的反射系數(shù)序列；如果不收斂，則將r0更新為當(dāng)前得到的r，返回步驟101，進(jìn)行下一輪迭代，直至盲反褶積目標(biāo)函數(shù)j(w,r)收斂，其中，盲反褶積目標(biāo)函數(shù)j(w,r)為：

j(w,r)＝j(luò)0(w,r)+μjr(r)+αji(r)，

其中，

在一個(gè)示例中，在子波序列求解單元中，可以采用共軛梯度迭代求解公式24以得到，迭代終止條件為：

其中，j^k(w,r)表示該共軛梯度迭代中第k次迭代后得到的j(w,r)，ε1為預(yù)設(shè)參數(shù)。

在一個(gè)示例中，在反射系數(shù)求解單元中，可以基于下式重新計(jì)算反射系數(shù)序列r：

r＝(g^tg+μq+αc^tc)^-1(g^td+αc^tξ)，公式31其中，q為對(duì)角矩陣，其第i行第i列的元素q(m_cauchy)ii為：

σr為給定值，表示在反射系數(shù)序列r的先驗(yàn)概率分布符合高斯分布的情況下反射系數(shù)序列r的方差。

在另一示例中，在反射系數(shù)求解單元中，可以基于下式重新計(jì)算反射系數(shù)序列r：

在這種情況下，可以采用共軛梯度迭代求解公式32以得到r，迭代終止條件為：

其中，j^k(w,r)表示該共軛梯度迭代中第k次迭代后得到的j(w,r)，ε2為預(yù)設(shè)參數(shù)。

應(yīng)用示例

為便于理解本發(fā)明實(shí)施例的方案及其效果，以下給出一個(gè)具體應(yīng)用示例。本領(lǐng)域技術(shù)人員應(yīng)理解，該示例僅為了便于理解本發(fā)明，其任何具體細(xì)節(jié)并非意在以任何方式限制本發(fā)明。

圖2示出了中國(guó)西部某工區(qū)的地震成像剖面。

圖3示出了采用傳統(tǒng)預(yù)測(cè)反褶積對(duì)圖2所示地震成像剖面進(jìn)行處理得到的結(jié)果。這種反褶積方法中，假設(shè)地震子波為最小相位，地震反射系數(shù)為白噪音。對(duì)比圖3和圖2可以看出，采用傳統(tǒng)預(yù)測(cè)反褶積可以使同相軸變細(xì)，使得整體剖面的分辨率得以提升，但其結(jié)果依然不能令人滿意。

圖4示出了采用本發(fā)明對(duì)圖2所示地震成像剖面進(jìn)行處理得到的結(jié)果。對(duì)比圖4和圖1，采用本發(fā)明可以使同相軸變細(xì)，提高剖面的分辨率。在對(duì)比圖4和圖3，可以看出，圖4中示出的同相軸更細(xì)，意味著分辨率更高，同時(shí)整體剖面更為干凈，意味著噪音更小，相比于圖3，圖4得到的反射系數(shù)估計(jì)更為準(zhǔn)確。

圖5示出了圖2、圖3、圖4的頻譜分析對(duì)比示意圖。通過(guò)對(duì)比頻譜可以看出，通過(guò)這兩種反褶積，均可展寬頻譜，因此可提升分辨率。相比而言，傳統(tǒng)預(yù)測(cè)反褶積的展寬效果較差，而采用本發(fā)明的展寬效果更為顯著，特別在低頻，端，可有效展寬了低頻部分，這對(duì)于當(dāng)前的地震反演成像等應(yīng)用都非常重要。

本發(fā)明可以是系統(tǒng)、方法和/或計(jì)算機(jī)程序產(chǎn)品。計(jì)算機(jī)程序產(chǎn)品可以包括計(jì)算機(jī)可讀存儲(chǔ)介質(zhì)，其上載有用于使處理器實(shí)現(xiàn)本發(fā)明的各個(gè)方面的計(jì)算機(jī)可讀程序指令。

計(jì)算機(jī)可讀存儲(chǔ)介質(zhì)可以是可以保持和存儲(chǔ)由指令執(zhí)行設(shè)備使用的指令的有形設(shè)備。計(jì)算機(jī)可讀存儲(chǔ)介質(zhì)例如可以是――但不限于――電存儲(chǔ)設(shè)備、磁存儲(chǔ)設(shè)備、光存儲(chǔ)設(shè)備、電磁存儲(chǔ)設(shè)備、半導(dǎo)體存儲(chǔ)設(shè)備或者上述的任意合適的組合。計(jì)算機(jī)可讀存儲(chǔ)介質(zhì)的更具體的例子(非窮舉的列表)包括：便攜式計(jì)算機(jī)盤、硬盤、隨機(jī)存取存儲(chǔ)器(ram)、只讀存儲(chǔ)器(rom)、可擦式可編程只讀存儲(chǔ)器(eprom或閃存)、靜態(tài)隨機(jī)存取存儲(chǔ)器(sram)、便攜式壓縮盤只讀存儲(chǔ)器(cd-rom)、數(shù)字多功能盤(dvd)、記憶棒、軟盤、機(jī)械編碼設(shè)備、例如其上存儲(chǔ)有指令的打孔卡或凹槽內(nèi)凸起結(jié)構(gòu)、以及上述的任意合適的組合。這里所使用的計(jì)算機(jī)可讀存儲(chǔ)介質(zhì)不被解釋為瞬時(shí)信號(hào)本身，諸如無(wú)線電波或者其他自由傳播的電磁波、通過(guò)波導(dǎo)或其他傳輸媒介傳播的電磁波(例如，通過(guò)光纖電纜的光脈沖)、或者通過(guò)電線傳輸?shù)碾娦盘?hào)。

這里所描述的計(jì)算機(jī)可讀程序指令可以從計(jì)算機(jī)可讀存儲(chǔ)介質(zhì)下載到各個(gè)計(jì)算/處理設(shè)備，或者通過(guò)網(wǎng)絡(luò)、例如因特網(wǎng)、局域網(wǎng)、廣域網(wǎng)和/或無(wú)線網(wǎng)下載到外部計(jì)算機(jī)或外部存儲(chǔ)設(shè)備。網(wǎng)絡(luò)可以包括銅傳輸電纜、光纖傳輸、無(wú)線傳輸、路由器、防火墻、交換機(jī)、網(wǎng)關(guān)計(jì)算機(jī)和/或邊緣服務(wù)器。每個(gè)計(jì)算/處理設(shè)備中的網(wǎng)絡(luò)適配卡或者網(wǎng)絡(luò)接口從網(wǎng)絡(luò)接收計(jì)算機(jī)可讀程序指令，并轉(zhuǎn)發(fā)該計(jì)算機(jī)可讀程序指令，以供存儲(chǔ)在各個(gè)計(jì)算/處理設(shè)備中的計(jì)算機(jī)可讀存儲(chǔ)介質(zhì)中。

用于執(zhí)行本發(fā)明操作的計(jì)算機(jī)程序指令可以是匯編指令、指令集架構(gòu)(isa)指令、機(jī)器指令、機(jī)器相關(guān)指令、微代碼、固件指令、狀態(tài)設(shè)置數(shù)據(jù)、或者以一種或多種編程語(yǔ)言的任意組合編寫(xiě)的源代碼或目標(biāo)代碼，所述編程語(yǔ)言包括面向?qū)ο蟮木幊陶Z(yǔ)言—諸如smalltalk、c++等，以及常規(guī)的過(guò)程式編程語(yǔ)言—諸如“c”語(yǔ)言或類似的編程語(yǔ)言。計(jì)算機(jī)可讀程序指令可以完全地在用戶計(jì)算機(jī)上執(zhí)行、部分地在用戶計(jì)算機(jī)上執(zhí)行、作為一個(gè)獨(dú)立的軟件包執(zhí)行、部分在用戶計(jì)算機(jī)上部分在遠(yuǎn)程計(jì)算機(jī)上執(zhí)行、或者完全在遠(yuǎn)程計(jì)算機(jī)或服務(wù)器上執(zhí)行。在涉及遠(yuǎn)程計(jì)算機(jī)的情形中，遠(yuǎn)程計(jì)算機(jī)可以通過(guò)任意種類的網(wǎng)絡(luò)—包括局域網(wǎng)(lan)或廣域網(wǎng)(wan)—連接到用戶計(jì)算機(jī)，或者，可以連接到外部計(jì)算機(jī)(例如利用因特網(wǎng)服務(wù)提供商來(lái)通過(guò)因特網(wǎng)連接)。在一些實(shí)施例中，通過(guò)利用計(jì)算機(jī)可讀程序指令的狀態(tài)信息來(lái)個(gè)性化定制電子電路，例如可編程邏輯電路、現(xiàn)場(chǎng)可編程門陣列(fpga)或可編程邏輯陣列(pla)，該電子電路可以執(zhí)行計(jì)算機(jī)可讀程序指令，從而實(shí)現(xiàn)本發(fā)明的各個(gè)方面。

這里參照根據(jù)本發(fā)明實(shí)施例的方法、裝置(系統(tǒng))和計(jì)算機(jī)程序產(chǎn)品的流程圖和/或框圖描述了本發(fā)明的各個(gè)方面。應(yīng)當(dāng)理解，流程圖和/或框圖的每個(gè)方框以及流程圖和/或框圖中各方框的組合，都可以由計(jì)算機(jī)可讀程序指令實(shí)現(xiàn)。

這些計(jì)算機(jī)可讀程序指令可以提供給通用計(jì)算機(jī)、專用計(jì)算機(jī)或其它可編程數(shù)據(jù)處理裝置的處理器，從而生產(chǎn)出一種機(jī)器，使得這些指令在通過(guò)計(jì)算機(jī)或其它可編程數(shù)據(jù)處理裝置的處理器執(zhí)行時(shí)，產(chǎn)生了實(shí)現(xiàn)流程圖和/或框圖中的一個(gè)或多個(gè)方框中規(guī)定的功能/動(dòng)作的裝置。也可以把這些計(jì)算機(jī)可讀程序指令存儲(chǔ)在計(jì)算機(jī)可讀存儲(chǔ)介質(zhì)中，這些指令使得計(jì)算機(jī)、可編程數(shù)據(jù)處理裝置和/或其他設(shè)備以特定方式工作，從而，存儲(chǔ)有指令的計(jì)算機(jī)可讀介質(zhì)則包括一個(gè)制造品，其包括實(shí)現(xiàn)流程圖和/或框圖中的一個(gè)或多個(gè)方框中規(guī)定的功能/動(dòng)作的各個(gè)方面的指令。

也可以把計(jì)算機(jī)可讀程序指令加載到計(jì)算機(jī)、其它可編程數(shù)據(jù)處理裝置、或其它設(shè)備上，使得在計(jì)算機(jī)、其它可編程數(shù)據(jù)處理裝置或其它設(shè)備上執(zhí)行一系列操作步驟，以產(chǎn)生計(jì)算機(jī)實(shí)現(xiàn)的過(guò)程，從而使得在計(jì)算機(jī)、其它可編程數(shù)據(jù)處理裝置、或其它設(shè)備上執(zhí)行的指令實(shí)現(xiàn)流程圖和/或框圖中的一個(gè)或多個(gè)方框中規(guī)定的功能/動(dòng)作。

附圖中的流程圖和框圖顯示了根據(jù)本發(fā)明的多個(gè)實(shí)施例的系統(tǒng)、方法和計(jì)算機(jī)程序產(chǎn)品的可能實(shí)現(xiàn)的體系架構(gòu)、功能和操作。在這點(diǎn)上，流程圖或框圖中的每個(gè)方框可以代表一個(gè)模塊、程序段或指令的一部分，所述模塊、程序段或指令的一部分包含一個(gè)或多個(gè)用于實(shí)現(xiàn)規(guī)定的邏輯功能的可執(zhí)行指令。在有些作為替換的實(shí)現(xiàn)中，方框中所標(biāo)注的功能也可以以不同于附圖中所標(biāo)注的順序發(fā)生。例如，兩個(gè)連續(xù)的方框?qū)嶋H上可以基本并行地執(zhí)行，它們有時(shí)也可以按相反的順序執(zhí)行，這依所涉及的功能而定。也要注意的是，框圖和/或流程圖中的每個(gè)方框、以及框圖和/或流程圖中的方框的組合，可以用執(zhí)行規(guī)定的功能或動(dòng)作的專用的基于硬件的系統(tǒng)來(lái)實(shí)現(xiàn)，或者可以用專用硬件與計(jì)算機(jī)指令的組合來(lái)實(shí)現(xiàn)。

以上已經(jīng)描述了本發(fā)明的各實(shí)施例，上述說(shuō)明是示例性的，并非窮盡性的，并且也不限于所披露的各實(shí)施例。在不偏離所說(shuō)明的各實(shí)施例的范圍和精神的情況下，對(duì)于本技術(shù)領(lǐng)域的普通技術(shù)人員來(lái)說(shuō)許多修改和變更都是顯而易見(jiàn)的。本文中所用術(shù)語(yǔ)的選擇，旨在最好地解釋各實(shí)施例的原理、實(shí)際應(yīng)用或?qū)κ袌?chǎng)中的技術(shù)的改進(jìn)，或者使本技術(shù)領(lǐng)域的其它普通技術(shù)人員能理解本文披露的各實(shí)施例。