：
1.本发明属于海底矿物质分布模拟方法，特别涉及一种模拟海底多金属结核赋存分布方法。
背景技术：2.占地球表面70％的海洋蕴藏着丰富的矿产资源。多金属结核是国际海底区域一种重要的海底固体矿产资源，富含锰、铜、钴、镍锂等金属元素。多金属结核呈大小不等的球状、不规则状，广泛分布于大洋海盆的沉积物表层。
3.我国在东太平洋cc区和西北太平洋已有3个多金属结核矿区，总面积超过21万平方公里，目前逐步从勘探阶段走向试采和开发阶段，详查和评价工作任务巨大。受限于当前调查程度较低和资源评价工作不足等背景，从检索的公开资料分析，目前尚未有真实、形象直观地反映海底多金属结核矿产赋存分布的技术方法。因此，有必要发明一种专门的反映结核赋存分布的技术方法，用于海底多金属结核矿产资源的评价和勘探开发工作。
技术实现要素：4.本发明的目的就在于提供一种模拟海底多金属结核赋存分布方法，该方法在海底多金属结核勘探开发过程中，能综合利用多种来源数据，有效提高海底多金属结核丰度的空间精度，更真实形象直观地反映海底结核矿产资源的赋存分布。
5.如上构思，本发明的技术方案是：一种模拟海底多金属结核赋存分布方法，包括如下步骤：
6.步骤一：收集研究区域的数据资料形成数据集且设置数学基础，所述研究区域的数据资料包括多金属结核地质采样站位数据、多波束调查数据、海底视像剖面和海底表层沉积物；所述数学基础是确定研究区域的地理范围、分辨率(网格大小r)、规定投影和坐标系，形成研究区网格area
study
；
7.步骤二：将多波束调查数据进行处理，即将海底水深数据网格化形成网格半径为r的网格数据g
dem
、将回波强度数据网格化形成网格半径为r的回波强度网格g
sidescan
，采用网格数据g
dem
生成地形起伏度网格g
bpi
和海底地形粗糙度网格g
rough
；
8.步骤三：将多金属结核地质采样站位数据进行处理，包括多金属结核地质采样站位数据网格化和多金属结核地质采样站位数据整理，所述地质采样站位数据整理即根据地质采样站位数据网格化网格数据，形成用于线性回归的站位数据集d；
9.步骤四：预测多金属结核丰度值g
iab
，即首先获得研究区域多金属结核预测丰度线性回归模型公式f
nodule
，，再调用此公式，利用步骤三得到的网格数据预测多金属结核丰度值g
iab
；
10.步骤五：对地质站位数据进行泰森多边形分析；
11.步骤六：以大小和数量不同的球体模型来模拟研究区所有格点的多金属结核赋存分布。
12.优选地，所述地质采样站位数据包括站位经纬度、多金属结核丰度、覆盖率、结核类型及大中小重量百分比。
13.优选地，所述步骤二的具体步骤包括：
14.①
利用多波束方法测定海底水深数据，采用距离加权法或高斯样条插值法将所测得的海底水深数据网格化，形成网格半径为r的网格数据g
dem
，计算公式如下：
15.其中，式中wi和disi分别为参与网格化点的水深值、计算权重值以及与交点的距离大小；
16.②
利用多波束方法测定回波强度数据，采用距离加权法或高斯样条插值法将所测得的回波强度数据网格化，形成网格半径为r的回波强度网格g
sidescan
；
17.③
采用网格数据g
dem
生成地形起伏度bpi网格g
bpi
和海底地形粗糙度网格g
rough
。
18.优选地，所述步骤三的具体步骤包括：
19.①
将多金属结核地质采样站位数据网格化，即采用克里金网格化方法将原始多金属结核地质采样站位数据的覆盖率、大型结核的重量百分比p
l
和小型结核的重量百分比ps进行网格化，生成覆盖率gcovg、大型结核重量百分比g
pl
和小型结核重量百分比g
ps
网格数据，其分辨率与g
dem
相同；
20.②
将多金属结核地质采样站位数据进行整理，即从多金属结核地质采样站位已有的站号、经纬度、水深、丰度、覆盖率、大型结核重量百分比p
l
和小型结核重量百分比ps生成的网格数据中查询提取站点的回波强度、bpi和粗糙度，形成用于线性回归的站位数据集d。
21.优选地，所述步骤四的具体步骤包括：
22.①
利用步骤三形成的站位数据集d，以站位丰度值为因变量，自变量从站位数据集d的覆盖率、回波强度、bpi、p
l
、ps、水深和粗糙度中选取，进行线性回归分析，求得研究区域多金属结核预测丰度线性回归模型公式f
nodule
；
23.②
调用研究区域多金属结核预测丰度线性回归模型公式f
nodule
，利用上述网格数据gcovg、g
sidescan
、g
dem
、g
bpi
、g
rough
、g
pl
和g
ps
计算研究区域所有格点的丰度预测值g
iab
。
24.优选地，所述步骤五的具体步骤包括：
25.①
采用步骤三形成的站位数据集d进行泰森多边形分析，形成泰森多边形poly
nodule
，每个多边形都包含结核大中小重量百分比属性数据；
26.②
对area
study
的每个网格单元与泰森多边形poly
nodule
叠加分析，使area
study
的每个网格单元都添加结核大中小重量百分比和结核丰度属性数据。
27.优选地，所述步骤六的具体步骤包括：
28.①
根据研究区域的丰度最大值ab
max
和研究区域的数学基础，选取球体模型的合适数量num
max
；
29.②
对area
study
的每个网格单元依据丰度预测值g
iab
、大型结核重量百分比g
pli
、小型结核重量百分比g
psi
、中型结核重量百分比(1-g
psi-g
pli
)，计算网格单元大型结核的数量num
pli
＝num
max
×
(g
iab
÷
ab
max
)
×gpli
、小型结核的数量num
psi
＝num
max
×
(g
iab
÷
ab
max
)
×gpsi
和中型结核的数量num
pmi
＝num
max
×
(g
iab
÷
ab
max
)
×
(1-g
psi-g
pli
)；
30.③
在area
study
的每个网格单元面积范围内，分别生成num
pli
、num
pmi
和num
psi
个随机坐标位置，形成包含坐标和大中小属性的数据文件file
data
。
31.优选地，所述步骤六中的步骤
③
的算法如下：
32.①
读取网格坐标、丰度值g
iab
、结核大中小重量百分比g
pli
、g
psi
和(1-g
psi-g
pli
)；
33.②
计算网格单元内的坐标点总数量n＝num
max
×
(g
iab
÷
ab
max
)，大型结核的数量num
pli
＝n
×gpli
，小型结核的数量num
psi
＝n
×gpsi
，中型结核的数量num
pmi
＝n
×
(1-g
psi-g
pli
)；
34.③
进行投影转换，生成网格区域rgn，其长宽为width和height；
35.④
以下是大型结核的坐标生成的算法伪代码：
[0036][0037][0038]
⑤
小型结核的坐标和中型结核的坐标生成算法同上；
[0039]
⑥
重复以上
①
至
⑤
，直至area
study
的每个网格单元都生成num
pli
、num
pmi
和num
psi
个随机坐标位置并形成包含坐标和大中小属性的数据文件file
data
为止。
[0040]
④
采用三维软件，导入g
dem
生成海底三维地形图，加载数据文件file
data
按照大中小确定球体半径和坐标确定位置进行叠加显示，生成研究区所有格点的多金属结核赋存分布模拟图。
[0041]
优选地，所述线性回归模型的线性回归相关系数r大于0.9。
[0042]
本发明实现了模拟海底结核矿产资源的赋存分布，并给出了详细的技术流程与方法，在海底多金属结核勘探开发过程中，能综合利用多种来源数据，有效提高海底多金属结核丰度的空间精度，更真实形象直观地反映海底结核矿产资源的赋存分布，具有较高的实际应用价值。
附图说明：
[0043]
图1是本发明的工作流程图；
[0044]
图2是本发明步骤六的工作流程图；
[0045]
图3是研究区所有格点的多金属结核赋存分布模拟图。
具体实施方式：
[0046]
参考附图，本发明提供一种模拟海底多金属结核赋存分布方法，是按照下述步骤建立的，如图1所示：
[0047]
步骤一：收集研究区域的数据资料，设置数学基础。
[0048]
1、收集研究区域的数据资料形成数据集：
[0049]
地质采样站位数据(包括站位经纬度、多金属结核丰度、覆盖率、结核类型及大中小重量比)、多波束调查数据、海底视像剖面、海底表层沉积物等相关数据；
[0050]
2、设置数学基础，确定研究区域的地理范围，分辨率(网格大小r)、规定投影和坐标系，形成研究区网格area
study
；
[0051]
步骤二：多波束调查数据处理：
[0052]
1、对于全覆盖的多波束方法测定海底水深数据，因其数据密集、数据量巨大，采用距离加权法或高斯样条插值法将所测得的海底水深数据网格化，形成网格半径为r的网格数据g
dem
。计算公式如下：
[0053]
其中，
[0054]
式中wi和disi分别为参与网格化点的水深值、计算权重值以及与交点的距离大小。
[0055]
2、利用多波束方法测定回波强度数据，采用距离加权法或高斯样条插值法将所测得的回波强度数据网格化，形成网格半径为r的回波强度网格g
sidescan
；
[0056]
3、采用g
dem
生成地形起伏度(又称地形位置指数(bpi))网格g
bpi
和海底地形粗糙度网格g
rough
；
[0057]
步骤三：多金属结核地质采样站位数据处理：
[0058]
1、地质采样站位数据网格化
[0059]
对原始地质采样站位数据采用克里金(kriging)网格化方法将地质采样站位数据的覆盖率(包括站位覆盖率、海底视像提取的结核覆盖率)、大型结核(结核直径大于6cm)重量比p
l
和小型结核(直径小于3cm)重量比ps进行网格化，生成覆盖率gcovg、大型结核重量比g
pl
和小型结核重量比g
ps
网格数据，其分辨率与g
dem
相同；
[0060]
2、地质采样站位数据整理
[0061]
地质采样站位已有站号、经纬度、水深、丰度、覆盖率、大型结核重量比p
l
、小型结核重量比ps，再从上述生成的网格数据中查询提取站点的回波强度、bpi、粗糙度，形成用于线性回归的站位数据集d；
[0062]
步骤四：预测多金属结核丰度值：
[0063]
1、利用站位数据集d，以站位丰度值为因变量，自变量从站位数据集d的覆盖率、回波强度、bpi、p
l
、ps、水深和粗糙度中选取，进行线性回归分析，一般要求线性回归相关系数r大于0.9，求得研究区域多金属结核预测丰度线性回归模型公式f
nodule
；
[0064]
2、调用研究区域多金属结核预测丰度线性回归模型公式f
nodule
，利用上述网格数据gcovg、g
sidescan
、g
dem
、g
bpi
、g
rough
、g
pl
和g
ps
计算研究区所有格点的丰度预测值；
[0065]
步骤五：地质站位数据进行泰森多边形分析：
[0066]
1、采用多金属结核站位数据集d进行泰森(thiessen)多边形分析，形成泰森多边形poly
nodule
，每个多边形都包含结核大中小重量百分比属性数据；
[0067]
2、对area
study
的每个网格单元与泰森多边形poly
nodule
叠加分析，使area
study
的每个网格单元都添加结核大中小重量百分比和结核丰度属性数据；
[0068]
步骤六：以大小和数量不同的球体模型来模拟研究区所有格点的多金属结核赋存分布：
[0069]
1、根据研究区的丰度最大值ab
max
和研究区比例尺等数学基础选取球体模型的合适数量num
max
；
[0070]
2、对area
study
的每个网格单元依据丰度预测值g
iab
、大型结核重量百分比g
pli
、小型结核重量百分比g
psi
、中型结核重量百分比(1-g
psi-g
pli
)，计算网格单元大型结核的数量num
pli
＝num
max
×
(g
iab
÷
ab
max
)
×gpli
，小型结核的数量num
psi
＝num
max
×
(g
iab
÷
ab
max
)
×gpsi
，中型结核的数量num
pmi
＝num
max
×
(g
iab
÷
ab
max
)
×
(1-g
psi-g
pli
)；
[0071]
3、在area
study
的每个网格单元面积范围内，分别生成num
pli
、num
pmi
和num
psi
个随机坐标位置，形成包含坐标和大中小属性的数据文件file
data
；
[0072]
其流程如图2，其算法如下：
[0073]
(1)读取网格坐标、丰度g
iab
、结核大中小重量百分比g
pli
、g
psi
和(1-g
psi-g
pli
)；
[0074]
(2)计算网格单元内的坐标点总数量n＝num
max
×
(g
iab
÷
ab
max
)，大型结核的数量num
pli
＝n
×gpli
，小型结核的数量num
psi
＝n
×gpsi
，中型结核的数量num
pmi
＝n
×
(1-g
psi-g
pli
)；
[0075]
(3)进行投影转换，生成网格区域rgn，其长宽为width和height；
[0076]
(4)以下是大型结核的坐标生成的伪代码：
[0077][0078]
(5)小型结核的坐标和中型结核的坐标生成算法同上；
[0079]
(6)重复以上(1)至(5)，直至area
study
的每个网格单元都生成num
pli
、num
pmi
和num
psi
个随机坐标位置并形成包含坐标和大中小属性的数据文件file
data
为止；
[0080]
4、采用三维软件，导入g
dem
生成海底三维地形图，加载数据文件file
data
按照结核大中小确定球状半径、坐标确定位置生成若干球体进行叠加显示，生成研究区所有格点的多金属结核赋存分布模拟图，如图3所示；
[0081]
5也可以研制三维软件系统，以地形g
dem
为基底，根据场景距离远近和比例尺动态选取球体模型的数量num
max
，采用步骤六的第2条和第3条中的方法动态计算相应每个网格单元面积范围内的大中小球体的数量和坐标位置生成球体并渲染，从而实现海底多金属结核赋存分布模拟技术。
[0082]
以上所述实施例仅表达了本技术的具体实施方式，其描述较为具体和详细，但并不能因此而理解为对本技术保护范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本技术技术方案构思的前提下，还可以做出若干变形和改进，这些都属于本技术的保护范围。
技术特征：1.一种模拟海底多金属结核赋存分布方法，其特征在于：包括如下步骤：步骤一：收集研究区域的数据资料形成数据集且设置数学基础，所述研究区域的数据资料包括多金属结核地质采样站位数据、多波束调查数据、海底视像剖面和海底表层沉积物；所述数学基础是确定研究区域的地理范围、分辨率(网格大小r)、规定投影和坐标系，形成研究区网格area
study
；步骤二：将多波束调查数据进行处理，即将海底水深数据网格化形成网格半径为r的网格数据g
dem
、将回波强度数据网格化形成网格半径为r的回波强度网格g
sidescan
，采用网格数据g
dem
生成地形起伏度网格g
bpi
和海底地形粗糙度网格g
rough
；步骤三：将多金属结核地质采样站位数据进行处理，包括多金属结核地质采样站位数据网格化和多金属结核地质采样站位数据整理，所述地质采样站位数据整理即根据地质采样站位数据网格化网格数据，形成用于线性回归的站位数据集d；步骤四：预测多金属结核丰度值g
iab
，即首先获得研究区域多金属结核预测丰度线性回归模型公式f
nodule，
，再调用此公式，利用步骤三得到的网格数据预测多金属结核丰度值g
iab
；步骤五：对地质站位数据进行泰森多边形分析；步骤六：以大小和数量不同的球体模型来模拟研究区所有格点的多金属结核赋存分布。2.根据权利要求1所述的一种模拟海底多金属结核赋存分布方法，其特征在于：所述地质采样站位数据包括站位经纬度、多金属结核丰度、覆盖率、结核类型及大中小重量百分比。3.根据权利要求1所述的一种模拟海底多金属结核赋存分布方法，其特征在于：所述步骤二的具体步骤包括：
①
利用多波束方法测定海底水深数据，采用距离加权法或高斯样条插值法将所测得的海底水深数据网格化，形成网格半径为r的网格数据g
dem
，计算公式如下：其中，式中w
i
和dis
i
分别为参与网格化点的水深值、计算权重值以及与交点的距离大小；
②
利用多波束方法测定回波强度数据，采用距离加权法或高斯样条插值法将所测得的回波强度数据网格化，形成网格半径为r的回波强度网格g
sidescan
；
③
采用网格数据g
dem
生成地形起伏度bpi网格g
bpi
和海底地形粗糙度网格g
rough
。4.根据权利要求1所述的一种模拟海底多金属结核赋存分布方法，其特征在于：所述步骤三的具体步骤包括：
①
将多金属结核地质采样站位数据网格化，即采用克里金网格化方法将原始多金属结核地质采样站位数据的覆盖率、大型结核的重量百分比p
l
和小型结核的重量百分比p
s
进行网格化，生成覆盖率gcovg、大型结核重量百分比g
pl
和小型结核重量百分比g
ps
网格数据，其分辨率与g
dem
相同；
②
将多金属结核地质采样站位数据进行整理，即从多金属结核地质采样站位已有的站号、经纬度、水深、丰度、覆盖率、大型结核重量百分比p
l
和小型结核重量百分比p
s
生成的网格数据中查询提取站点的回波强度、bpi和粗糙度，形成用于线性回归的站位数据集d。
5.根据权利要求1所述的一种模拟海底多金属结核赋存分布方法，其特征在于：所述步骤四的具体步骤包括：
①
利用步骤三形成的站位数据集d，以站位丰度值为因变量，自变量从站位数据集d的覆盖率、回波强度、bpi、p
l
、p
s
、水深和粗糙度中选取，进行线性回归分析，求得研究区域多金属结核预测丰度线性回归模型公式f
nodule
；
②
调用研究区域多金属结核预测丰度线性回归模型公式f
nodule
，利用上述网格数据gcovg、g
sidescan
、g
dem
、g
bpi
、g
rough
、g
pl
和g
ps
计算研究区域所有格点的丰度预测值g
iab
。6.根据权利要求1所述的一种模拟海底多金属结核赋存分布方法，其特征在于：所述步骤五的具体步骤包括：
①
采用步骤三形成的站位数据集d进行泰森多边形分析，形成泰森多边形poly
nodule
，每个多边形都包含结核大中小重量百分比属性数据；
②
对area
study
的每个网格单元与泰森多边形poly
nodule
叠加分析，使area
study
的每个网格单元都添加结核大中小重量百分比和结核丰度属性数据。7.根据权利要求1所述的一种模拟海底多金属结核赋存分布方法，其特征在于：所述步骤六的具体步骤包括：
①
根据研究区域的丰度最大值ab
max
和研究区域的数学基础，选取球体模型的合适数量num
max
；
②
对area
study
的每个网格单元依据丰度预测值g
iab
、大型结核重量百分比g
pli
、小型结核重量百分比g
psi
、中型结核重量百分比(1-g
psi-g
pli
)，计算网格单元大型结核的数量num
pli
＝num
max
×
(g
iab
÷
ab
max
)
×
g
pli
、小型结核的数量num
psi
＝num
max
×
(g
iab
÷
ab
max
)
×
g
psi
和中型结核的数量num
pmi
＝num
max
×
(g
iab
÷
ab
max
)
×
(1-g
psi-g
pli
)；
③
在area
study
的每个网格单元面积范围内，分别生成num
pli
、num
pmi
和num
psi
个随机坐标位置，形成包含坐标和大中小属性的数据文件file
data
。8.根据权利要求8所述的一种模拟海底多金属结核赋存分布方法，其特征在于：所述步骤六中的步骤
③
的算法如下：
①
读取网格坐标、丰度值g
iab
、结核大中小重量百分比g
pli
、g
psi
和(1-g
psi-g
pli
)；
②
计算网格单元内的坐标点总数量n＝num
max
×
(g
iab
÷
ab
max
)，大型结核的数量num
pli
＝n
×
g
pli
，小型结核的数量num
psi
＝n
×
g
psi
，中型结核的数量num
pmi
＝n
×
(1-g
psi-g
pli
)；
③
进行投影转换，生成网格区域rgn，其长宽为width和height；
④
以下是大型结核的坐标生成的算法伪代码：
⑤
小型结核的坐标和中型结核的坐标生成算法同上；
⑥
重复以上
①
至
⑤
，直至area
study
的每个网格单元都生成num
pli
、num
pmi
和num
psi
个随机坐标位置并形成包含坐标和大中小属性的数据文件file
data
为止。
④
采用三维软件，导入g
dem
生成海底三维地形图，加载数据文件file
data
按照大中小确定球体半径和坐标确定位置进行叠加显示，生成研究区所有格点的多金属结核赋存分布模拟图。9.根据权利要求5所述的一种模拟海底多金属结核赋存分布方法，其特征在于：所述线性回归模型的线性回归相关系数r大于0.9。
技术总结一种模拟海底多金属结核赋存分布方法，包括如下步骤：1、收集研究区域的数据资料形成数据集且设置数学基础；2、将多波束调查数据进行处理；3、将多金属结核地质采样站位数据进行处理；4、预测多金属结核丰度值G
技术研发人员：程永寿 相文玺 章伟艳 程淦 李正刚 李小虎 翟洪昌
受保护的技术使用者：国家海洋信息中心
技术研发日：2022.03.18
技术公布日：2022/7/5
