c3演算法

① 斷層信息提取研究進展

斷裂對油氣聚集帶或油氣藏控製作用的實質則是斷裂封閉性問題，斷層在油氣運移和聚集過程中既是運移通道又是聚集成藏的封堵屏障。斷層活動期油氣可通過斷層進入新的儲層中，斷層成為油氣運移的通道；斷層停止活動時則阻擋了油氣的進一步運移，在斷層靠近油源一側的構造高部位聚集成藏。因此，研究斷裂的封閉性非常重要，其中，斷層檢測與信息提取則是斷裂封閉性評價研究的基礎。

斷層信息的識別與解釋是含油氣盆地地震構造解釋的關鍵。常規的三維地震資料一體化解釋方法是通過相干體技術、可視化技術、層位自動追蹤技術、面塊切片技術、地震屬性處理功能等手段，採用點、線、面結合的空間解釋，在建立斷層骨架的基礎上，再對小斷層進行解釋，最終使斷層和層位達到線道閉合。這種方法大大提高了斷層解釋的精度和效率，為復雜斷塊的小斷層、小斷塊圈閉識別的可靠性提供了依據。

但是，隨著我國油氣田勘探和開發強度的增加，再發現大的構造圈閉幾率已愈來愈小，而隱蔽圈閉（subtle trap）與油氣藏的勘探越來越占據重要地位。其中，岩性圈閉和小斷層（裂縫）的識別是尋找隱蔽圈閉和油氣藏的關鍵。因此，綜合利用地質、鑽井、測井、測試、高解析度三維地震資料等，運用一些信息處理新方法新技術對岩性圈閉和小斷層識別預測。在此基礎上，進一步研究岩性圈閉和小斷層的控油氣作用，這不僅可以為我國復雜隱蔽油氣藏的勘探提供技術上和理論分析方面的支持，而且對於我國石油工業增產上儲具有重要意義。

由此可見，在地震構造解釋中，斷層的檢測和識別非常重要，尤其是對那些隱蔽小斷層（裂縫）的識別。在這方面，作者通過參與國家「十五」科技攻關計劃項目「中國油氣資源發展關鍵技術研究」 的「隱蔽油氣藏描述新方法新技術研究」課題研究工作提出了以下三種提取斷層信息的新方法。

（1）基於高階統計量的相干體估計演算法

相干就是地震資料的連續性分析。通過計算相鄰道與道之間的相似性，把計算出的值以最大值或最小值的方式計算形成一個新的數據體，突出同相軸的突變點或不連續點。當地下存在斷層和其他地層不連續變化時，在這些不連續點的兩側，地震道會表現出不同的反射波特徵，從而導致局部道與道之間的相乾性突變。通過計算道與道之間的這種不相似性來描述地層、岩性的非均一性。這種技術人為因素相對小，其水平時間切片顯示了任意方向的斷層，可克服平行同相軸的斷層難以解釋的困難，提供另一個視角更直觀的觀測結果。

1995年Bahorich等提出的相干體演算法是地震資料解釋發展史上的一個重大突破，相干體能夠很好地描述地層、岩性的橫向不均勻性，使斷層的自動解釋成為可能。第一代相干體演算法是基於二階統計量的C1演算法，這一演算法簡單、高效、易於實現，但它在處理帶雜訊數據時的穩健性不好。

Marfurt等（1998）又提出了利用多個地震道的C2演算法，這一演算法可以利用任意數目的地震道來估計相乾性。與C1演算法相比，這一演算法對雜訊的抑制能力明顯提高。但是，在用C2演算法計算相乾性時，為了抑制雜訊，必須用更多的地震道，所以得到的相干體的解析度會降低，同時，計算量也會顯著增加。

Gersztenkorn等（1999）提出了一種基於特徵結構的相干估計演算法（C3）。這一方法首先構建一包含若干條地震道的三維數據塊並用這一數據塊構成一協方差矩陣，然後用協方差矩陣的最大特徵值和矩陣跡的比值作為相干估計值。實踐證明，C3演算法比C1、C2演算法更能抑制數據中的雜訊。然而，由於C3演算法要構建規模較大的協方差矩陣並計算其主特徵值，所以計算量很大。C3演算法的另一缺點是沒有考慮傾角對相干估計值的影響，因此，在有很強傾角的區域，C3演算法不能給出好的相干估計值。Marfurt（1999）提出改進的C3演算法，這一方法先用C2演算法估計出傾角並對傾角估計值進行平滑，然後用平滑後的傾角值計算相干值，但這一改進演算法的計算量很大。

Cohen（2002）提出一種估計局部結構熵的演算法（LSE），這一方法先構建一三維數據塊，然後把這一三維數據塊分成四個子塊，並用這四個子塊構建一個（4×4）的相關矩陣。最後，用歸一化的矩陣的跡作為結構熵的估計值。與C3演算法相比，LSE演算法不須計算規模較大的協方差矩陣的主特徵值，因而計算量較小。但是，LSE演算法同樣未考慮傾角對局部結構熵的影響。

王西文等（2002）提出了一種在小波變換域中進行相干分析的演算法。這一方法在不同的頻帶進行相干分析，對不同的頻帶得到的相干體進行加強或減弱，然後利用所有頻帶的相干體進行重構。基於小波變換的方法能有效地提高斷層解釋的解析度。

由於高階統計量能保持信號相位並且能夠自動抑制信號中方差未知的加性高斯雜訊，因此，它被成功地應用於信號處理領域。相干分析的一個核心問題是時延估計，高階統計量在相關高斯雜訊情況下可以比二階統計量給出更好的時延估計結果。Srinivasan（2001）利用高階統計量進行地震同相軸的自動拾取，相對二階統計量，高階統計量在處理地震信號的時延估計問題時可以給出更高解析度的時延估計。傳統的基於互相關的C1演算法先分步估計不同空間方向地震信號的時延並得到其對應的相關度，然後加以融合。與C1演算法不同的是，基於高階統計量的相干體估計演算法同時利用三條地震道計算它們的四階矩的一個零時延切片，然後在該二維切片上尋找最大相關點，並將該點的四階矩作為相干估計值。為了在相干分析中同時利用更多的地震道，提出了一種用多個地震道逐點交織組合構成超道的方法，和原始地震道一樣，新形成的超道同樣可以很方便地進行傾角掃描.利用超道技術和HOSC、C1進行組合，可以得到兩種高效的相干分析演算法，分別稱之為ST -HOSC和ST-C1。利用實際資料處理結果表明，該演算法獲得的相干體包含更多的斷層細節信息。

（2）偏移後三維地震數據體空間去模糊法

偏移後三維地震數據體為地震解釋提供了大量的信息，但由於偏移運算元的孔徑影響，偏移得到的三維地震數據體在空間上被模糊了，不利於進一步的斷層檢測。

為了獲得高解析度的斷層信息，一方面可以對相干分析演算法本身進行研究；另一方面，可以通過提高輸入數據體的空間解析度來達到同樣的目的。為了提高地震數據體的空間解析度，我們採用了圖像去模糊技術。所謂的圖像去模糊技術，就是假設觀測圖像可以用一個多維褶積模型表示，即由真實圖像和一個點擴散函數（Point Spread Function，PSF）褶積得到，而去模糊就是一個多維反褶積過程，以消除PSF影響，恢復真實圖像。我們假設三維地震數據體的每一個時間切片，都是由地下結構圖像和一個二維的PSF褶積得到的，且對於所有的時間切片，PSF都是不變的。另外，為了保持原數據體的相位信息，設定所利用的PSF是零相位。利用所有時間切片估計出PSF後，基於維納濾波的反褶積方法被採用來去掉PSF的影響，得到提高空間解析度的數據體。實際資料處理結果表明，運用該演算法先對地震數據體進行預處理，然後再進行相干分析，可以得到高解析度的相干體圖像，達到提取更多斷層信息的目的。

基於空間反褶積的偏移後三維地震數據體預處理技術，來壓制在空間上對應地層變化的低頻成分，而突出對應斷層的高頻成分。該方法可以有效地增強地震數據體中的對應斷層的信號，可以得到更高解析度、細節更豐富的相干體圖像。

（3）基於相干濾波的相干體圖像增強法

在以往所有的相干分析演算法中，均需要進行精細的參數選擇，兩個最重要的參數是所用的地震道數和分析時窗長度。為了壓制沿地層形成的低相干帶和隨機雜訊帶來的干擾，現有的相干分析技術都需要用多個地震道和較大的分析時窗。但多道地震信號和大分析時窗的採用，則會導致所獲得的相干體的空間和時間解析度降低，同時，採用太多地震道和太大的分析時窗，將導致斷層空間上的位移及擴展，後者將同樣造成假斷層現象。為了獲得高解析度的相干體圖像，一個簡單可行的方法是用少的地震道和小的分析時窗，而且使用少的地震道和小的分析時窗，可以降低計算量，但其缺點是產生了強的沿地層的低相干帶，從而在時間水平切片上形成假的斷層。通過分析，我們發現在縱向相干體剖面上，沿地層形成的低相干帶和由斷層產生的低相干帶存在方向上的差別。一般來說，沿地層形成的低相干帶是接近水平方向的，而由斷層產生的低相干帶是接近垂直方向的。利用上述兩種低相干帶的方向上的差異，提出了一種基於相干濾波的相干體圖像增強技術，通過對相干體進行後濾波處理來壓制不需要的低相干帶和隨機雜訊帶來的干擾，從而使得在相干分析時可以使用較少的地震道和較短的分析時窗來獲取高解析度和高信噪比斷層信息。本書提出的方法假設由斷層產生的低相干帶在局部范圍內是線性相干信號，從而通過局部方向掃描和多項式濾波來估計由斷層產生的低相干帶，達到濾除不需要的低相干帶和隨機雜訊帶來的干擾的目的。實際資料處理結果表明，該演算法可以處理任何相干分析演算法採用相對少的地震道和小分析時窗時獲得的相干體，得到真實的高解析度和高信噪比的相干體圖像。

② 有6級台階,我從下向上走,若每次只能跨一級或兩級,我走上去有幾種不同走法

這是個排列組合問題。
解題方法是：
每次都跨一級：1種
有一次跨兩級：C5 1=5種
有兩次跨兩級：C4 2=6種
有三次跨兩級：C3 3=1種
1+5+6+1=13種
這就是最後答案。
C5 1寫法是在C的右面上標寫1，下標寫5。
C5 1表示的是把同一次跨過的那兩級台階算作一個，這樣就一共有5級台階，那合在一起的台階共有5種放法。其他幾個同理。
C的演算法是（拿C4 2、C3 3舉例）：C4 2=（4*3）/（2*1）=6、C3 3=（3*2*1）/（3*2*1）=1

至於八級台階，照著上面的方法推就行了，只給你個式子:
1+C7 1+C6 2+C5 3+C4 4=34（種）