一种基于地震波波动方程的照明分析计算方法及计算装置

摘要

本发明提供一种基于地震波波动方程的照明分析计算方法包括：建立地质模型，读入观测系统参数；将地质模型离散化成网格模型；读取第一层的数据准备进行计算；设定一个初始波场；计算下一层相应位置处的波场；根据计算得到的下一层各点的波场值，计算下一层各点的照明能量；计算装置输出能量分布图。本发明提供的一种基于地震波波动方程的照明分析计算方法，炮点和检波点位置相同的点只计算一次，减少了计算量；计算时，将波场值存储在计算机内存中，降低了对计算机存储空间的要求并提高了计算的速度。本发明还提供了一种基于地震波波动方程的照明分析计算装置实施例，该实施例可以实现本发明基于地震波波动方程的照明分析计算。

主权项

一种基于地震波波动方程的照明分析计算方法，其特征在于，包括：S1：计算装置根据工区地质构造建立地质模型，并读入地震观测系统的参数；S2：计算装置根据地震观测系统，将建立的地质模型离散化成多层的网格模型；S3：读取上述网格模型的第一层的数据至内存，所述第一层作为当前层，准备进行计算；S4：计算装置以炮点或检波点位置为中心，在计算区域内设定一个初始波场；S5：计算装置根据初始波场和传播算子，计算当前层的下一层相应位置处的波场；S6：计算装置判断当前作为中心的炮点或检波点是否是观测系统中当前层的最后一个炮点或检波点，并根据判断的结果执行相应的操作；S7：计算装置根据计算得到的所述当前层的下一层各点的波场值，计算所述当前层的下一层各点的照明能量；S8：计算装置判断计算出照明能量的层是否是网格模型的最后一层，并根据判断的结果执行相应的操作；S9：计算装置输出每一层的能量分布图；其中，所述S1中观测系统的参数，具体包括：地震波频率、炮点位置和检波点位置；所述S6，具体包括：判断当前作为中心的炮点或检波点是否是观测系统中当前层的最后一个炮点或检波点，若是最后一个炮点或检波点，则丢弃内存中原来存储的当前层的波场值，将计算所得的当前层的下一层的波场值存储在内存中；若当前作为中心的作用点不是观测系统中当前层的最后一个炮点或检波点，则将上述S4中的中心向右偏移至下一个炮点或检波点，并重复S4‑S6进行计算和判断，直至该层所有观测系统中的炮点和检波点都计算结束，丢弃内存中原来存储的当前层的波场值，将计算所得的当前层的下一层的波场值存放在内存中；所述S8，具体包括：先判断S7计算出的照明能量是否为网格模型的最后一层的照明能量，如是，则计算结束；若计算出照明能量的层不是网格模型的最后一层，则将该计算出照明能量的层作为当前层，并重复S4‑S8进行计算和判断，直至计算出网格模型中最后一层的照明能量，则计算结束；所述S7中计算装置根据计算得到的下一层各点的波场值，计算下一层各点的照明能量，包括：先将S5计算得到的下一层的波数域波场值作傅立叶变换至空间域；然后在空间域根据观测系统中的炮检关系，对于每一个作用点，先计算出该点作为中心点时在下一层的相同位置的波场值，然后计算出与其相关的接收点作为中心点时，该点在各区域计算得到的下一层相同位置的波场值，将该点作为中心点时得到的下一层的波场值进行平方后分别和与其相关的接收点作为中心点时该点的下一层波场值的平方相乘并求和，计算所得的结果记作该计算区域中心点位置的照明能量DI，并存储在计算机的硬盘上。