| 主权项 |
一种地震三维观测系统的共中心点CMP面元属性的分析方法，其特征在于，包括：根据计算资源所允许容纳的数据量，将地震三维观测系统的待勘探区域的面元格网划分为多个覆盖次数计算分区；采用多个计算资源并行分别计算多个覆盖次数计算分区中面元的覆盖次数；根据所述多个覆盖次数计算分区中面元的覆盖次数和计算资源所允许容纳的中点数量，将所述待勘探区域的所述面元格网划分为多个中点属性计算分区；采用多个计算资源并行分别计算多个中点属性计算分区的中点属性值，其中，所述中点属性值包括：炮检距、炮检方位角；根据计算资源所允许容纳的数据量，将所述地震三维观测系统的待勘探区域的面元格网划分为多个覆盖次数计算分区，包括：在大地坐标系下，对所述待勘探区域的面元格网中的面元进行编号，其中，以所述大地坐标系原点处的面元为起点，先按照以RMin所在的方向从左到右的顺序，再按照RMax所在的方向从下到上的顺序对所述面元格网中的面元进行递增编号，所述RMin是所述待勘探区域横向长度的最小者，所述RMax是所述待勘探区域纵向长度的最大者；根据以下公式确定覆盖次数计算分区的数量：<img file="FDA0000854496880000011.GIF" wi="838" he="151" />其中，所述M是覆盖次数计算分区的数量，<maths num="0001"><math><![CDATA[<mrow><mi>λ</mi><mo>=</mo><mfrac><mrow><mo>(</mo><mi>M</mi><mi>a</mi><mi>x</mi><mi>B</mi><mi>i</mi><mi>n</mi><mi>R</mi><mi>o</mi><mi>w</mi><mi>I</mi><mi>n</mi><mi>l</mi><mo>+</mo><mn>1</mn><mo>)</mo><mo>×</mo><mo>(</mo><mi>M</mi><mi>a</mi><mi>x</mi><mi>B</mi><mi>i</mi><mi>n</mi><mi>C</mi><mi>o</mi><mi>l</mi><mi>C</mi><mi>r</mi><mi>l</mi><mo>+</mo><mn>1</mn><mo>)</mo><mo>×</mo><mi>F</mi><mi>o</mi><mi>l</mi><mi>d</mi><mi>D</mi><mi>a</mi><mi>t</mi><mi>a</mi><mi>P</mi><mi>e</mi><mi>r</mi><mi>B</mi><mi>i</mi><mi>n</mi></mrow><mrow><mi>D</mi><mi>a</mi><mi>t</mi><mi>a</mi><mi>P</mi><mi>e</mi><mi>r</mi><mi>K</mi><mi>e</mi><mi>r</mi><mi>n</mi><mi>e</mi><mi>l</mi></mrow></mfrac><mo>,</mo></mrow>]]></math><img file="FDA0000854496880000012.GIF" wi="1325" he="128" /></maths><img file="FDA0000854496880000013.GIF" wi="70" he="78" />为向上取整符号，NumOfKernel参与计算的计算资源的数量，DataPerKernel为每个计算资源所允许容纳的数据量，FoldDataPreBin为每个面元的覆盖次数属性项数据，MaxBinRowInl为所述待勘探区域的纵向最大面元行号，MaxBinColCrl为所述待勘探区域的横向最大面元列号；根据计算资源所允许容纳的数据量，将所述地震三维观测系统的待勘探区域的面元格网划分为多个覆盖次数计算分区，包括：通过以下公式确定每个覆盖次数计算分区包括的面元行数：BinLinePerZone＝min(L<sub>max</sub>,η)，其中，<maths num="0002"><math><![CDATA[<mrow><mi>η</mi><mo>=</mo><mfrac><mrow><mi>D</mi><mi>a</mi><mi>t</mi><mi>a</mi><mi>P</mi><mi>e</mi><mi>r</mi><mi>K</mi><mi>e</mi><mi>r</mi><mi>n</mi><mi>e</mi><mi>l</mi></mrow><mrow><mi>F</mi><mi>o</mi><mi>l</mi><mi>d</mi><mi>D</mi><mi>a</mi><mi>t</mi><mi>a</mi><mi>P</mi><mi>e</mi><mi>r</mi><mi>B</mi><mi>i</mi><mi>n</mi><mo>×</mo><mrow><mo>(</mo><mi>M</mi><mi>i</mi><mi>n</mi><mi>B</mi><mi>i</mi><mi>n</mi><mi>I</mi><mi>D</mi><mi>O</mi><mi>f</mi><mi>L</mi><mi>i</mi><mi>n</mi><mi>e</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow></mrow></mfrac><mo>,</mo></mrow>]]></math><img file="FDA0000854496880000021.GIF" wi="925" he="135" /></maths><img file="FDA0000854496880000022.GIF" wi="629" he="147" /><img file="FDA0000854496880000023.GIF" wi="68" he="71" />为向上取整符号，BinLinePerZone为每个覆盖次数计算分区包括的面元行数，min为取二个数值中的较小者，DataPerKernel为每个计算资源所允许容纳的数据量，FoldDataPreBin为每个面元的覆盖次数属性项数据，MaxBinIDOfLine是每个覆盖次数计算分区的面元行数和列数中的最大者，MinBinIDOfLine是每个覆盖次数计算分区的面元行数和列数中的最小者，NumOfKernel是参与计算的计算资源的数量；通过以下公式确定每个覆盖次数计算分区的起始面元的编号：BinIDBegin<sub>i</sub>＝BinLinePerZone×(MinBinIDOfLine+1)×i，其中，所述i是覆盖次数计算分区的编号，i∈[0,M‑1]；通过以下公式确定每个覆盖次数计算分区的终止面元的编号：BinIDEnd<sub>i</sub>＝BinLinePerZone×(MinBinIDOfLine+1)×(i+1)‑1；根据所述多个覆盖次数计算分区中面元的覆盖次数和所述计算资源所允许容纳的中点数量，将所述待勘探区域的所述面元格网划分为所述多个中点属性计算分区，包括：根据以下公式确定每个中点属性计算分区的起始面元的编号：BinIDBegin<sub>k+1</sub>＝BinIDEnd<sub>k</sub>+1，其中，所述K是所述中点属性计算分区的编号，k∈[0,N‑1]，BinIDBegin<sub>k+1</sub>为第k+1个中点属性计算分区的起始面元的编号，BinIDEnd<sub>k</sub>为第k个中点属性计算分区的终止面元的编号；根据以下公式确定每个中点属性计算分区的终止面元的编号：BinIDEnd<sub>k+1</sub>＝BinIDBegin<sub>k+1</sub>+R‑1，其中，BinIDEnd<sub>k+1</sub>为第k+1个中点属性计算分区的终止面元的编号，所述R为所述中点属性计算分区包括的面元数。 |
