[發明專利]一種地震正演模擬中的多級斷點保護方法

有效

• 申請號: 201710486742.3
• 申請日: 2017-06-23
• 公開（公告）號: CN109116417B
• 公開（公告）日: 2020-05-19
• 發明（設計）人: 王小青;楊文廣;曹永生;楊尚琴;陳金煥;陳海洋
• 申請（專利權）人: 中國石油化工股份有限公司;中國石油化工股份有限公司石油物探技術研究院
• 主分類號: G01V1/30
• 分類號: G01V1/30
• 代理公司: 北京思創畢升專利事務所 11218
• 代理人: 孫向民;廉莉莉

1.一種地震正演模擬中的多級斷點保護方法，其特征在于，所述方法包括：

1)用戶輸入參數DCHECKPT的值；

2)對所輸入的參數DCHECKPT的值進行判斷：

如果DCHECKPT＝0，則進入炮級別的斷點保護，即主進程在每炮循環的末尾處向文件CHECKPT_INF中覆蓋輸出當前炮的激發編號SHOTID；

如果DCHECKPT＝N，且N為正整數，則進入時間級別的斷點保護，即當it％DCHECKPT＝0時，主進程在時間循環的末尾處向文件CHECKPT_INF中覆蓋輸出斷點時間nt_checkpt，當it％DCHECKPT≠0時不輸出；所有進程在各自的節點上生成獨立的CHECKPT_FILE文件，并且向各自的CHECKPT_FILE文件中覆蓋輸出正演模擬所需的所有中間變量及數組信息；以及主進程在炮循環的末尾處向文件CHECKPT_INF中覆蓋輸出當前炮的激發編號SHOTID。

2.根據權利要求1所述的地震正演模擬中的多級斷點保護方法，其特征在于，重啟正演模擬運算時，執行以下步驟：

1’)用戶輸入參數CHECKPTREAD的值；

2’)對所輸入的參數CHECKPTREAD的值進行判斷：

如果CHECKPTREAD＝0，從第一個炮點開始進行炮點循環；

如果CHECKPTREAD＝1，則執行下一步驟；

3’)對所輸入的參數DCHECKPT的值進行判斷：

如果DCHECKPT＝0，則所有進程在炮循環開始前讀入文件CHECKPT_INF中所保存的SHOTID值，并令重啟的正演模擬作業的激發循環從SHOTID+1開始；

如果DCHECKPT＝N，則所有進程在炮循環開始前讀入文件CHECKPT_INF中所保存的SHOTID值，并令重啟的正演模擬作業的激發循環從SHOTID+1開始；所有進程在時間循環開始前讀入文件CHECKPT_INF中保存的nt_checkpt值，并令時間循環從nt_checkpt+1開始；所有進程讀入各自的CHECKPT_FILE文件中的中間變量及數組的值賦給對應的變量和數組；

4’)重復執行步驟1)和步驟2)，其中，nt_checkpt為斷點時間。

3.根據權利要求2所述的地震正演模擬中的多級斷點保護方法，其特征在于，在步驟4’)之前執行以下步驟：

31’)用戶輸入參數CHECKPTWRITE的值；

32’)對所輸入的參數CHECKPTWRITE的值進行判斷：

如果CHECKPTWRITE＝0，則不輸出斷點保護信息；

如果CHECKPTWRITE＝1，繼續執行步驟4’)。

4.根據權利要求1所述的地震正演模擬中的多級斷點保護方法，其特征在于，在DCHECKPT＝N的情況下，正演模擬所需的所有中間變量及數組信息N個時間間隔更新一次。

5.根據權利要求2所述的地震正演模擬中的多級斷點保護方法，其特征在于，對地震記錄進行實時計算并實時輸出到地震記錄文件中：

當ishot＝1且it＝1時，從地震記錄文件的起始位置依次輸出3600字節的卷頭、240字節的道頭和地震記錄；

當ishot≠1且it＝1時輸出240字節的道頭和地震記錄；以及

當it≠1時輸出地震記錄，其中，ishot為當前震源的序號，it為當前的時間采樣點的序號。