[發明專利]基于遺傳算法的剩余靜校正求取方法

在審

• 申請號: 201510945879.1
• 申請日: 2015-12-16
• 公開（公告）號: CN105608493A
• 公開（公告）日: 2016-05-25
• 發明（設計）人: 周強;張曉斌;何光明;陳愛萍;劉鴻;曹琳昱
• 申請（專利權）人: 中國石油集團川慶鉆探工程有限公司地球物理勘探公司
• 主分類號: G06N3/12
• 分類號: G06N3/12;G06N3/00
• 代理公司: 成都中璽知識產權代理有限公司 51233
• 代理人: 安宇宏;譚昌馳

技術領域

本發明涉及一種基于遺傳算法的剩余靜校正求取方法，該方法可自動、 科學確定遺傳算法計算參數，應用于石油地震勘探，屬于地震勘探資料處理 領域，是一種獲取精確的高頻靜校正量的有效方法。

背景技術

我國年石油天然氣對外依存度極高，而我國剩余油氣資源超過一半聚集 在山地復雜構造中，及時發現探明山地復雜構造油氣資源迫在眉睫。這些地 區地震地表條件復雜，靜校正問題嚴重，導致多次覆蓋資料同一反射點因為 地震道之間存在延遲時差，疊加成像困難，其精度難以滿足后期的構造、儲 層解釋，由這類品質的資料提供的建議井位可靠性不高，將會增加油氣勘探 的成本和風險。

為降低這類因素帶來的成本和風險，物探人員嘗試采用遺傳算法求剩余 靜校正量，減少共反射點道集的延遲時差，其基本思想是從模型群體開始搜 索，把模型參數編碼，將模型空間的點映射到染色體空間的染色體，然后通 過選擇、交換和變異等遺傳操作對模型群體進行繁殖，逐次迭代，最后求取 非線性反演問題整體極值所對應的最優解或近似最優解。

然而，常規的遺傳算法求取剩余靜校正的實際性能受控制參數(如種群數 目、種群大小、變異概率等)的影響，如何確定這些控制參數成為一直困擾物 探人員的瓶頸問題，也直接制約了該技術在生產上的推廣應用。

發明內容

為解決常規遺傳算法參數難以確定的問題，申請人提出一種基于遺傳算 法的剩余靜校正求取方法，可以根據允許的最大時移量計算出種群大小、種 群代數和個體變異率，自動計算出剩余靜校正量。

本發明提供了一種基于遺傳算法的剩余靜校正求取方法。所述方法用于 處理二維地震資料，其包括以下步驟：

A、將不同的炮點或檢波點的真實值分別作為遺傳算法的不同個體中的 單個基因，將多個個體構成初始種群，種群大小popSize通過式1確定，式1 為：popSize＝(maxSta/sampInv)×f_popSize，式1中maxSta為最大時移量， 在10～60ms之間選擇；sampInv為地震數據的采樣間隔，單位為ms；f_popSize為 種群大小的調控因子，在0.5～1.5之間取值；

B、計算種群中每個個體的適應度值，隨后對種群中的所有個體按適應度 值由大到小排序，選擇保留適應度值最高的個體作為交換對象；

C、在種群的第k個個體和所述交換對象上各隨機產生兩個交換位置，將 所述第k個個體的兩個交換位置之間的基因段與所述交換對象的兩個交換位 置之間的基因段進行交換，產生子代新個體，其中，k＝popSize×0.4，k須 取整；

D、對步驟C中所述進行交換后獲得的種群，保留適應度值最高的個體作為精英進入子代，對于其余個體，若隨機數r大于該個體的變異概率，則該個體隨機產生兩個變異位置，并將該個體的兩個變異位置處的基因隨機變為允許移動的最大時移量maxSta范圍內的隨機值，從而產生子代新個體，若隨機數r不大于該個體的變異概率，則不變異，同時該個體直接進入子代作為子代個體，其中，隨機數r在[0，1]間隨機產生，個體的變異概率通過以下方式獲得，先對該個體所在種群中的所有個體按適應度值由大到小排序，排序后，第i個個體的變異概率根據式2得到，式2為：