CN111159869B - 一种砂岩铀矿床与古河道间定量参数的自动计算方法 - Google Patents
一种砂岩铀矿床与古河道间定量参数的自动计算方法 Download PDFInfo
- Publication number
- CN111159869B CN111159869B CN201911327726.5A CN201911327726A CN111159869B CN 111159869 B CN111159869 B CN 111159869B CN 201911327726 A CN201911327726 A CN 201911327726A CN 111159869 B CN111159869 B CN 111159869B
- Authority
- CN
- China
- Prior art keywords
- vector
- ancient river
- point
- sliding window
- vector point
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 229910052770 Uranium Inorganic materials 0.000 title claims abstract description 102
- JFALSRSLKYAFGM-UHFFFAOYSA-N uranium(0) Chemical compound [U] JFALSRSLKYAFGM-UHFFFAOYSA-N 0.000 title claims abstract description 101
- 238000004364 calculation method Methods 0.000 title claims abstract description 35
- 230000033558 biomineral tissue development Effects 0.000 claims abstract description 78
- 238000005452 bending Methods 0.000 claims abstract description 60
- 238000000034 method Methods 0.000 claims abstract description 17
- 230000000694 effects Effects 0.000 claims description 6
- 238000011166 aliquoting Methods 0.000 claims description 3
- 238000003491 array Methods 0.000 claims description 3
- 238000005553 drilling Methods 0.000 abstract description 10
- 238000010586 diagram Methods 0.000 description 10
- 238000007418 data mining Methods 0.000 description 3
- 238000000605 extraction Methods 0.000 description 3
- 230000011218 segmentation Effects 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000009189 diving Effects 0.000 description 1
- 239000011229 interlayer Substances 0.000 description 1
- 239000010410 layer Substances 0.000 description 1
- 230000004807 localization Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 239000005416 organic matter Substances 0.000 description 1
- 230000003647 oxidation Effects 0.000 description 1
- 238000007254 oxidation reaction Methods 0.000 description 1
- 229910052704 radon Inorganic materials 0.000 description 1
- SYUHGPGVQRZVTB-UHFFFAOYSA-N radon atom Chemical compound [Rn] SYUHGPGVQRZVTB-UHFFFAOYSA-N 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 239000013049 sediment Substances 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Landscapes
- Complex Calculations (AREA)
Abstract
本发明属于参数计算领域，具体涉及一种砂岩铀矿床与古河道间定量参数的自动计算方法,该方法包括按照一定距离等分古河道矢量线段线状要素，记录等分古河道的矢量点要素集合，通过滑窗方法自动计算古河道的弯曲度系数集合，通过求解弯曲度系数集合的局部极大值，得到古河道局部极大弯曲位置。本发明通过读入铀矿化矢量点数据，可以自动计算铀矿化矢量点要素位置距离最近古河道弯曲位置的最小距离参数，自动计算计算钻孔矢量点要素位置距离最近古河道矢量线要素的最小距离参数,自动且定量地计算古河道与铀矿之间的特征参数，所求得的特征参数更加精确，速度也更快。
Description
技术领域
本发明属于参数计算领域，具体涉及一种砂岩铀矿床与古河道间定量参数的自动计算方法。
背景技术
古河道砂岩型铀矿是指产于古河道沉积砂岩层中的一种后生铀矿床，铀主要经过潜水/潜水-层间氧化作用富集成矿。该类铀矿单个矿床一般规模较小，品位较低，但矿体具成群成带产出特征，分布面积达若干平方公里，总储量较大。因此，古河道砂岩型铀矿具有较好的找矿潜力。由于其可采用成本低，古河道砂岩型铀矿是铀工业中非常具有经济价值和竞争力的铀资源之一。
铀资源勘查数据挖掘应用最终目的是分析和挖掘多源异构的铀资源数据的隐含信息，得到规律性认识，为铀成矿预测服务。需要从这些数据中辨析和抽取成矿信息并与矿化信息集成后才可以应用，即对数据预处理中抽取、辨析方法进行研究。如提取铀矿化或钻孔与古河道位置关系参数、砂体有机质含量、岩性结构、土壤氡气、高精度磁测、水中铀、地化异常U、Th、K等数值、航磁、航放数值等等，与矿化级别的矿化品位、矿化厚度等数值进行集成，完成试验区铀成矿信息大数据集合建设，为下一步的数据挖掘和预测奠定基础。其中，定量提取铀矿钻孔与古河道之间关系参数非常重要。
古河道型砂岩铀矿化在空间上受古河道控制、并严格定位于某范围内的砂岩铀成床，铀矿化多产于古河道的特殊部位。现在分析铀矿钻孔位置与古河道分布关系，往往采用定性的方式描述，如已知钻孔位于古河道的拐弯、交汇或变缓部位。这种定性的描述很难运用数据挖掘方法进行已知钻孔与古河道关系的分析。
因此需要提出一种砂岩铀矿床与古河道间定量参数的自动计算方法，来自动求解古河道局部弯曲部位，并计算钻孔位置与古河道局部极大弯曲部位距离、钻孔位置与古河道之间最小距离等重要特征参数。
发明内容
本发明的目的在于提供一种定量计算铀矿钻孔与古河道之间特征参数的方法，通过自动求解古河道局部弯曲部位，自动计算铀矿化位置与古河道最大弯曲部位距离、铀矿化位置与古河道之间最小距离等重要参数，用于定量描述古河道与铀矿之间的关系，为后期定量分析铀矿化与古河道等找矿要素之间关系提供数据支撑。
本发明的技术方案是：
一种砂岩铀矿床与古河道间定量参数的自动计算方法，该方法包括以下步骤：
步骤(1)按照一定距离等分古河道矢量线段线状要素；
步骤(2)记录步骤(1)中获取的等分的位置，得到等分古河道的矢量点要素集合；
步骤(3)针对步骤(2)获取的等分古河道的矢量点要素集合，通过滑窗方法自动计算古河道的弯曲度系数集合；
步骤(4)利用步骤(3)获取的古河道的弯曲度系数集合，求解弯曲度系数集合的局部极大值，得到古河道局部极大弯曲位置；
步骤(5)读入铀矿化矢量点数据；
步骤(6)计算步骤(5)读入的铀矿化矢量点要素位置距离步骤(4)获取的最近古河道弯曲位置的最小距离参数；
步骤(7)计算步骤(5)读入的铀矿化矢量点要素位置距离步骤(2)获取的最近古河道等分位置矢量点的最小距离参数。
所述的步骤(1)中将古河道矢量线段线状要素按照一定距离等分，古河道矢量线状要素记为L，将线状要素L等分(s-1)段，s为等分点的总点数，s值可根据实际情况调整。
所述的步骤(2)中记录等分古河道线状要素位置的矢量点要素集合，古河道矢量线状要素L等分为(s-1)条线段，其等分位置矢量点集合记为Ai(xi，yi)，其中i为古河道等分位置矢量点编号，s为等分矢量点总点数，i∈[1，s]，xi为编号为i的矢量点要素的经度值，yi为编号为i的矢量点要素的纬度值。
所述的步骤(3)具体包含以下子步骤：
步骤(3.1)确定滑窗大小m，即滑窗内点要素的数量值，设为m，m∈[3，s-1]，m取奇数，具体滑窗大小数值可根据效果调整；
步骤(3.2)按照步长1滑动滑窗，遍历矢量点集合Ai(xi，yi)，依次计算滑窗范围内点集合中第b点到第(m+b)点的欧式距离之和，使用一维数组N(k)进行记录，其中，k为各滑动窗口的集合序号，b为滑动窗口中第一个矢量点的序号，当滑动窗口步长为1时，k值与b值相等，s为滑窗的大小，即滑窗内点要素的数量值，设该一维数组元素总数为K，K＝s-m+1，k∈[1，s-m+1]，N(k)的计算公式为：
步骤(3.3)按照步长1滑动滑窗，遍历矢量点集合Ai(xi，yi)，依次计算滑窗内首尾矢量点间的欧式距离，即计算Ai(xi，yi)中第b点Ab(xb，yb)与第(m+b)点Am+b(xm+b，ym+b)两点之间的欧式距离n(k)，即使用一维数组n(k)进行记录该距离数值，其中k为各滑动窗口的集合序号，b为各滑动窗口中第1个矢量点的编号，计算公式为：
步骤(3.4)依次计算一维数组N(k)与n(k)之间各元素的比值，得到弯曲度系数集合数组一维数组w(k)，该数组中的元素可以反映滑窗范围内古河道的弯曲程度，其中k＝s-m+1，计算公式为：
w(k)＝N(k)/n(k)
所述的步骤(4)具体包含以下子步骤：
步骤(4.1)设置一维数组WQ(e)用来记录弯曲度系数集合w(k)中极大值的数值序号，其中，e为求解出的弯曲度极大值序号，e∈[1，n]，n为求解出的弯曲度极大值的总量，遍历弯曲度系数集合w(k)，k为各滑动窗口的集合序号，循环计算w(k+1)-w(k)和w(k+2)-w(k+1)，k∈[1，s-m+1]，当同时满足w(k+1)-w(k)>0且w(k+2)-w(k+1)>0且w(k+1)>a时，记录该一维数组w(k)的极大值数值序号，即WQ(e)＝k+1，a为满足古河道弯曲率的最小数值，设置a值可以去除一些弯曲率不显著的部位，a值可根据实际效果自行设置，a∈(1，+∞]；
步骤(4.2)通过弯曲度系数集合中极大值的数值序号数组WQ(e)，求解处于弯曲极值处矢量点序号，设弯曲极值点矢量点序号数组为d(e)，d(e)＝WQ(e)+INT(m/2)，其中，s为滑动窗口大小，INT(m/2)值为窗口大小数值除以2，并进行取整计算，e为求解出的弯曲度极大值序号，e∈[1，n]，n为求解出的弯曲度极大值的总量；
步骤(4.3)根据求解处于弯曲极值处矢量点序号数组d(e)，得到古河道局部极大弯曲位置矢量点集合Bf(xf，yf)＝Ad(e)(xd(e)，yd(e))其中，f为古河道局部极大弯曲位置矢量点集合编号，f∈[1，n]，n为求解出的弯曲度极大值的总量。
所述的步骤(5)中读入铀矿化矢量点数据；矢量点集合记为Zj(xxj，yyj)，其中j为铀矿化矢量点编号，铀矿化矢量点总数量为p，则j∈[1，p]，xxj为编号为j的铀矿化矢量点的经度值，yyj为编号为j的的铀矿化矢量点的纬度值。
所述的步骤(6)中计算铀矿化矢量点要素集合距离最近古河道弯曲位置矢量点集合的最小距离参数；即计算铀矿化矢量点要素集合Zj(xxj，yyj)到古河道弯曲位置矢量点集合Bf(xf，yf)的最小欧式距离集合，该最小欧式距离集合用一维数组G(j)表示，j为铀矿化矢量点编号，铀矿化矢量点总数量为p，j∈[1，p]，一维数组G(j)中的各元素反映了已知矿化点集合与其最邻近古河道弯曲拐点的距离参数，遍历铀矿化矢量点要素集合Zj(xxj，yyj)，j∈[1，p]，依次计算与古河道弯曲位置矢量点集合Bf(xf，yf)中各点的欧式距离，f∈[1，n]，并将求解的最小值记录在一维数组G(j)中。
所述的步骤(7)中计算铀矿化矢量点要素集合距离最近古河道古河道等分位置矢量点的最小距离参数；即计算铀矿化矢量点要素集合Zj(xxj，yyj)到古河道等分位置矢量点集合Ai(xi，yi)的最小欧式距离集合，该最小欧式距离集合用一维数组H(j)表示，j为铀矿化矢量点编号，铀矿化矢量点总数量为p，j∈[1，p]，一维数组H(j)中的各元素反映了已知矿化点集合与其最邻近古河道的距离参数，遍历铀矿化矢量点要素集合Zj(xxj，yyj)，j∈[1，p]，依次计算与古河道等分位置矢量点集合Ai(xi，yi)中各点的欧式距离，i∈[1，s]，并将求解的最小值记录在一维数组H(j)中。
本发明的有益效果：
本发明属于一种参数提取方法，具体公开一种砂岩铀矿床与古河道间定量参数的自动计算方法。该方法包括按照一定距离等分古河道矢量线段线状要素，记录等分古河道的矢量点要素集合，通过滑窗方法自动计算古河道的弯曲度系数集合。通过求解弯曲度系数集合的局部极大值，得到古河道局部极大弯曲位置；本发明读入铀矿化矢量点数据，自动计算铀矿化矢量点要素位置距离最近古河道弯曲位置的最小距离参数，自动计算计算钻孔矢量点要素位置距离最近古河道矢量线要素的最小距离参数。
通过本方法，自动且定量地计算古河道与铀矿之间的特征参数，所求得的特征参数更加精确，速度也更快。
附图说明
图1为本发明所提供的古河道的等分矢量点示意图；
图2为本发明所提供的第一个滑动窗口内等分矢量点求和计算示意图；
图3为本发明所提供的第二个滑动窗口内等分矢量点求和计算示意图；
图4为本发明所提供的第三个滑动窗口内等分矢量点求和计算示意图；
图5为本发明所提供的第一个滑动窗口内首尾矢量点欧式距离计算示意图；
图6为本发明所提供的第五个滑动窗口内等分矢量点求和计算示意图；
图7为本发明所提供的第五个滑动窗口内首尾矢量点欧式距离计算示意图；
图8为本发明所提供的求解处于古河道弯曲极值处的矢量点序号示意图；
图9为本发明所提供的矿化点到最近弯曲点距离及最近等分矢量点的距离参数示意图；
具体实施方式
下面结合附图与实施例对本发明进行进一步的介绍：
本发明提供了一种砂岩铀矿化与古河道间定量参数计算方法，包含以下步骤：
步骤(1)按照一定距离等分古河道矢量线段线状要素；
将古河道矢量线段线状要素按照一定距离等分。古河道矢量线状要素记为L，将线状要素L等分(s-1)段，s为等分点的总点数，s值可根据实际情况调整。
读入古河道矢量线段文件“古河道.shp”,该文件的投影坐标系统为“Beijing_1954_GK_Zone_15N”，即北京54大地坐标系，6度分带，分带号为15。该文件包含一条线状要素，即古河道的分布位置，记为L。使用GIS功能中的分割功能，将L平均分割为(s-1)＝25条线段。在实际应用中可以根据效果调整分割线段数值。该线段分割功能为常规GIS功能。
步骤(2)记录步骤(1)中获取的等分的位置，得到等分古河道的矢量点要素集合；
古河道矢量线状要素L等分为(s-1)条线段，其等分位置矢量点集合记为Ai(xi,yi)，其中i为古河道等分位置矢量点编号，s为等分矢量点总点数，i∈[1,s]。xi为编号为i的矢量点要素的经度值，yi为编号为i的矢量点要素的纬度值，如图1所示。
使用GIS功能自动记录古河道线状要素L分割后的端点位置，如图1。其等分位置矢量点集合记为Ai(xi,yi)，其中i为古河道等分位置矢量点编号，i∈[1,26]，即L被分割为25条线段，分割后的端点数s＝26。这26个矢量点要素集合具体坐标值表如表1。
表1等分位置矢量点集合坐标值
步骤(3)针对步骤(2)获取的等分古河道的矢量点要素集合，通过滑窗方法自动计算古河道的弯曲度系数集合；该步骤包含以下子步骤：
步骤(3.1)确定滑窗大小m，即滑窗内点要素的数量值，设为m，m∈[3,s-1]，m取奇数，具体滑窗大小数值可根据效果调整。
步骤(3.2)按照步长1滑动滑窗，遍历矢量点集合Ai(xi,yi)，依次计算滑窗范围内点集合中第b点到第(m+b)点的欧式距离之和,使用一维数组N(k)进行记录。其中，k为各滑动窗口的集合序号，b为滑动窗口中第一个矢量点的序号，当滑动窗口步长为1时，k值与b值相等。m为滑窗的大小，即滑窗内点要素的数量值。设该一维数组元素总数为K，K＝s-m+1，k∈[1,s-m+1]。则在此实例中，K＝s-m+1＝26-5+1＝22,k∈[1,22]，在本实例中，有26个矢量点要素集合，滑动窗口宽度为5，则滑动窗口总数为22；b为各滑动窗口中第1个矢量点的编号，b∈[1,22]。N(k)的计算公式为：
其中，i为古河道等分位置矢量点编号。
由于各矢量点是等间距的，滑窗大小为5，因此滑窗5个矢量点的欧式距离之和是一样的，在本实例中，一维数组N(k)各元素为：
N(k)＝[8071.4,8071.4,8071.4,8071.4,8071.4,8071.4,8071.4,8071.4,8071.4,8071.4,8071.4,8071.4,8071.4,8071.4,8071.4,8071.4,8071.4,8071.4,8071.4,8071.4,8071.4,8071.4]。
图2、图3、图4和图6分别第一、第二、第三和第五滑动窗口内等分矢量点欧式距离求和示意图，可观察滑动窗口的运动规律。
步骤(3.3)按照步长1滑动滑窗，遍历矢量点集合Ai(xi,yi)，依次计算滑窗内首尾矢量点间的欧式距离，即计算Ai(xi,yi)中第b点Ab(xb,yb)与第(m+b)点Am+b(xm+b,ym+b)两点之间的欧式距离n(k)，即使用一维数组n(k)进行记录该距离数值，其中k为各滑动窗口的集合序号，b为各滑动窗口中第1个矢量点的编号，计算公式为：
通过计算，在本实例中，一维数组n(k)各元素为：
n(k)＝[6378.7，6435.6，6129.1，4930.2，3915.9，6432.8，6073.9，6265.5，6292.3，6302.9，5854.7，5157.0，5036.6，5765.3，6287.6，5829.3，4333.7，3733.5，5183.2，6082.1，5642.9，5637.2]。
图5和图7分别为第一和第五个滑动窗口内首尾矢量点欧式距离计算示意图。
步骤(3.4)依次计算一维数组N(k)与n(k)之间各元素的比值，得到弯曲度系数集合数组一维数组w(k)，该数组中的元素可以反映滑窗范围内古河道的弯曲程度，其中k＝s-m+1，计算公式为：
w(k)＝N(k)/n(k)
通过计算，在本实例中，一维数组w(k)各元素为：
w(k)＝[1.265，1.254，1.317，1.637，2.061，1.255，1.329，1.288，1.283，1.281，1.379，1.565，1.603，1.400，1.284，1.385，1.862，2.162，1.557，1.327，1.430，1.432]
步骤(4)利用步骤(3)获取的古河道的弯曲度系数集合，求解弯曲度系数集合的局部极大值，得到古河道局部极大弯曲位置；
步骤(4.1)设置一维数组WQ(e)用来记录弯曲度系数集合w(k)中极大值的数值序号，其中，e为求解出的弯曲度极大值序号，e∈[1,n]，n为求解出的弯曲度极大值的总量。遍历弯曲度系数集合w(k)，k为各滑动窗口的集合序号，循环计算w(k+1)-w(k)和w(k+2)-w(k+1)，在本实例中，k∈[1,22]，当同时满足w(k+1)-w(k)>0且w(k+2)-w(k+1)>0且w(k+1)>1.5时，记录该一维数组w(k)的极大值数值序号，即WQ(e)＝k+1。
在数组w(k)中，符号上述条件的元素分别为：w(5)＝2.016；w(13)＝1.603；w(18)＝2.163。
因此，弯曲度系数集合中极大值的数值序号数组WQ(e)＝[5,13,18]。
步骤(4.2)通过弯曲度系数集合中极大值的数值序号数组WQ(e)，求解处于弯曲极值处矢量点序号。设弯曲极值点矢量点序号数组为d(e)，
d(e)＝WQ(e)+INT(m/2)，其中，m为滑动窗口大小，INT(m/2)值为窗口大小数值除以2，并进行取整。e为求解出的弯曲度极大值序号，e∈[1,n]，n为求解出的弯曲度极大值的总量。在本实例中，m＝5,则d(e)＝WQ(e)+INT(m/2)＝d(e)＝WQ(e)+INT(5/2)＝WQ(e)+2；其中，WQ(e)＝[5,13,18]；则d(e)＝WQ(e)+2＝[7,15,20]。即在本实例中，处于弯曲极值处矢量点序号分别为A7,A15和A20，如图8所示。
步骤(4.3)根据求解处于弯曲极值处矢量点序号数组d(e)＝[7,15,20]，得到古河道局部极大弯曲位置矢量点集合Bf(xf,yf)＝Ad(e)(xd(e),yd(e)),其中，f为古河道局部极大弯曲位置矢量点集合编号，f∈[1,n]。n为求解出的弯曲度极大值的总量，在本实例中，n＝3。处于古河道弯曲极值处矢量点集合Bf(xf,yf)中各元素如表2。
表2古河道弯曲极值处矢量点集合坐标值
步骤(5)读入铀矿化矢量点数据；
本实例中，读入铀矿化矢量点文件“矿化.shp”，如图9。矿化点集合记为Zj(xxj,yyj)，其中j为铀矿化矢量点编号，铀矿化矢量点总数量为3个，则j∈[1,3]。xxj为编号为j的铀矿化矢量点的经度值，yyj为编号为j的的铀矿化矢量点的纬度值，矿化矢量点的具体坐标值如表3。
表3读入的矿化矢量点集合坐标值
步骤(6)计算步骤(5)读入的铀矿化矢量点要素位置距离步骤(4)获取的最近古河道弯曲位置的最小距离参数；
在本实例中，即计算铀矿化矢量点要素集合Zj(xxj,yyj)到古河道弯曲位置矢量点集合Bf(xf,yf)的最小欧式距离集合。该最小欧式距离集合用一维数组G(j)表示，j为铀矿化矢量点编号，铀矿化矢量点总数量为3，j∈[1,3]。
遍历铀矿化矢量点要素集合Zj(xxj,yyj)，j∈[1,3]，依次计算与古河道弯曲位置矢量点集合Bf(xf,yf)中各点的欧式距离，f∈[1,n]，并将求解的最小值记录在一维数组G(j)中。
在本实例中，铀矿化矢量点集合中编号为1的矢量点Z1(xx1,yy1)到古河道弯曲位置矢量点集合Bf(xf,yf)的各点欧式距离分别为7315米、10942米和15769米，取其中的最小值7315米，则G(1)＝7315。
同理，铀矿化矢量点集合中编号为2的矢量点Z2(xx2,yy2)到古河道弯曲位置矢量点集合Bf(xf,yf)的各点欧式距离分别为2958米、9292米和11911米，取其中的最小值2958米，则G(2)＝2958。
同理，铀矿化矢量点集合中编号为3的矢量点Z3(xx3,yy3)到古河道弯曲位置矢量点集合Bf(xf,yf)的各点欧式距离分别为8695米、3385米和5827米，取其中的最小值2958米，则G(3)＝3385。
由此，在本实例中，铀矿化矢量点要素集合距离古河道弯曲位置矢量点集合的最小距离参数集合为G(j)＝[7315,2958,3385]，如图9所示。
步骤(7)计算步骤(5)读入的铀矿化矢量点要素位置距离步骤(2)获取的最近古河道等分位置矢量点的最小距离参数。
在本实例中，即计算铀矿化矢量点要素集合Zj(xxj,yyj)到古河道等分位置矢量点集合Ai(xi,yi)的最小欧式距离集合。该最小欧式距离集合用一维数组H(j)表示，j为铀矿化矢量点编号，铀矿化矢量点总数量为3，j∈[1,3]。一维数组H(j)中的各元素反映了已知矿化点集合与其最邻近古河道的距离参数。
遍历铀矿化矢量点要素集合Zj(xxj,yyj)，j∈[1,p]，依次计算与古河道等分位置矢量点集合Ai(xi,yi)中各点的欧式距离，i∈[1,26]，并将求解的最小值记录在一维数组H(j)中。
求解方法类似于步骤六，在此不再赘述，铀矿化矢量点要素集合距离古河道等分位置矢量点的最小距离参数集合为H(j)＝[1724，1275，2546]，如图9所示。
上面结合实施例对本发明作了详细说明，但是本发明并不限于上述实施例，在本领域普通技术人员所具备的知识范围内，还可以在不脱离本发明宗旨的前提下作出各种变化。本发明中未作详细描述的内容均可以采用现有技术。
Claims (6)
1.一种砂岩铀矿床与古河道间定量参数的自动计算方法，其特征在于，该方法包括以下步骤：
步骤(1)按照一定距离等分古河道矢量线段线状要素；
步骤(2)记录步骤(1)中获取的等分的位置，得到等分古河道的矢量点要素集合；
步骤(3)针对步骤(2)获取的等分古河道的矢量点要素集合，通过滑窗方法自动计算古河道的弯曲度系数集合；
所述的步骤(3)具体包含以下子步骤：
步骤(3.1)确定滑窗大小m，即滑窗内点要素的数量值，设为m，m∈[3，s-1]，m取奇数，具体滑窗大小数值可根据效果调整；
步骤(3.2)按照步长1滑动滑窗，遍历矢量点集合Ai(xi，yi)，依次计算滑窗范围内点集合中第b点到第(m+b)点的欧式距离之和，使用一维数组N(k)进行记录，其中，k为各滑动窗口的集合序号，b为滑动窗口中第一个矢量点的序号，当滑动窗口步长为1时，k值与b值相等，s为滑窗的大小，即滑窗内点要素的数量值，设该一维数组元素总数为K，K＝s-m+1，k∈[1，s-m+1]，N(k)的计算公式为：
步骤(3.3)按照步长1滑动滑窗，遍历矢量点集合Ai(xi，
yi)，依次计算滑窗内首尾矢量点间的欧式距离，即计算Ai(xi，yi)中第b点Ab(xb，yb)与第(m+b)点Am+b(xm+b，ym+b)两点之间的欧式距离n(k)，即使用一维数组n(k)进行记录该距离数值，其中k为各滑动窗口的集合序号，b为各滑动窗口中第1个矢量点的编号，计算公式为：
步骤(3.4)依次计算一维数组N(k)与n(k)之间各元素的比值，得到弯曲度系数集合数组一维数组w(k)，该数组中的元素可以反映滑窗范围内古河道的弯曲程度，其中k＝s-m+1，计算公式为：
w(k)＝N(k)/n(k)
步骤(4)利用步骤(3)获取的古河道的弯曲度系数集合，求解弯曲度系数集合的局部极大值，得到古河道局部极大弯曲位置；
所述的步骤(4)具体包含以下子步骤：
步骤(4.1)设置一维数组WQ(e)用来记录弯曲度系数集合w(k)中极大值的数值序号，其中，e为求解出的弯曲度极大值序号，e∈[1，n]，n为求解出的弯曲度极大值的总量，遍历弯曲度系数集合w(k)，k为各滑动窗口的集合序号，循环计算w(k+1)-w(k)和w(k+2)-w(k+1)，k∈[1，s-m+1]，当同时满足w(k+1)-w(k)>0且w(k+2)-w(k+1)>0且w(k+1)>a时，记录该一维数组w(k)的极大值数值序号，即WQ(e)＝k+1，a为满足古河道弯曲率的最小数值，设置a值可以去除一些弯曲率不显著的部位，a值可根据实际效果自行设置，a∈(1，+∞]；
步骤(4.2)通过弯曲度系数集合中极大值的数值序号数组WQ(e)，求解处于弯曲极值处矢量点序号，设弯曲极值点矢量点序号数组为d(e)，d(e)＝WQ(e)+INT(m/2)，其中，s为滑动窗口大小，INT(m/2)值为窗口大小数值除以2，并进行取整计算，e为求解出的弯曲度极大值序号，e∈[1，n]，n为求解出的弯曲度极大值的总量；
步骤(4.3)根据求解处于弯曲极值处矢量点序号数组d(e)，得到古河道局部极大弯曲位置矢量点集合Bf(xf，yf)＝Ad(e)(x d(e)，yd(e))其中，f为古河道局部极大弯曲位置矢量点集合编号，f∈[1，
n]，n为求解出的弯曲度极大值的总量；
步骤(5)读入铀矿化矢量点数据；
步骤(6)计算步骤(5)读入的铀矿化矢量点要素位置距离步骤(4)获取的最近古河道弯曲位置的最小距离参数；
步骤(7)计算步骤(5)读入的铀矿化矢量点要素位置距离步骤(2)获取的最近古河道等分位置矢量点的最小距离参数。
2.如权利要求1所述的一种砂岩铀矿床与古河道间定量参数的自动计算方法，其特征在于：所述的步骤(1)中将古河道矢量线段线状要素按照一定距离等分，古河道矢量线状要素记为L，将线状要素L等分(s-1)段，s为等分点的总点数，s值可根据实际情况调整。
3.如权利要求2所述的一种砂岩铀矿床与古河道间定量参数的自动计算方法，其特征在于：所述的步骤(2)中记录等分古河道线状要素位置的矢量点要素集合，古河道矢量线状要素L等分为(s-1)条线段，其等分位置矢量点集合记为Ai(xi，yi)，其中i为古河道等分位置矢量点编号，s为等分矢量点总点数，i∈[1，s]，xi为编号为i的矢量点要素的经度值，yi为编号为i的矢量点要素的纬度值。
4.如权利要求3所述的一种砂岩铀矿床与古河道间定量参数的自动计算方法，其特征在于：所述的步骤(5)中读入铀矿化矢量点数据；矢量点集合记为Zj(xxj，yyj)，其中j为铀矿化矢量点编号，铀矿化矢量点总数量为p，则j∈[1，p]，xxj为编号为j的铀矿化矢量点的经度值，yyj为编号为j的铀矿化矢量点的纬度值。
5.如权利要求4所述的一种砂岩铀矿床与古河道间定量参数的自动计算方法，其特征在于：所述的步骤(6)中计算铀矿化矢量点要素集合距离最近古河道弯曲位置矢量点集合的最小距离参数；即计算铀矿化矢量点要素集合Zj(xxj，yyj)到古河道弯曲位置矢量点集合Bf(xf，yf)的最小欧式距离集合，该最小欧式距离集合用一维数组G(j)表示，j为铀矿化矢量点编号，铀矿化矢量点总数量为p，j∈[1，p]，一维数组G(j)中的各元素反映了已知矿化点集合与其最邻近古河道弯曲拐点的距离参数，遍历铀矿化矢量点要素集合Zj(xxj，yyj)，j∈[1，p]，依次计算与古河道弯曲位置矢量点集合Bf(xf，yf)中各点的欧式距离，f∈[1，n]，并将求解的最小值记录在一维数组G(j)中。
6.如权利要求5所述的一种砂岩铀矿床与古河道间定量参数的自动计算方法，其特征在于：所述的步骤(7)中计算铀矿化矢量点要素集合距离最近古河道等分位置矢量点的最小距离参数；即计算铀矿化矢量点要素集合Zj(xxj，yyj)到古河道等分位置矢量点集合Ai(xi，yi)的最小欧式距离集合，该最小欧式距离集合用一维数组H(j)表示，j为铀矿化矢量点编号，铀矿化矢量点总数量为p，j∈[1，p]，一维数组H(j)中的各元素反映了已知矿化点集合与其最邻近古河道的距离参数，遍历铀矿化矢量点要素集合Zj(xxj，yyj)，j∈[1，p]，依次计算与古河道等分位置矢量点集合Ai(xi，yi)中各点的欧式距离，i∈[1，s]，并将求解的最小值记录在一维数组H(j)中。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911327726.5A CN111159869B (zh) | 2019-12-20 | 2019-12-20 | 一种砂岩铀矿床与古河道间定量参数的自动计算方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911327726.5A CN111159869B (zh) | 2019-12-20 | 2019-12-20 | 一种砂岩铀矿床与古河道间定量参数的自动计算方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111159869A CN111159869A (zh) | 2020-05-15 |
| CN111159869B true CN111159869B (zh) | 2023-11-14 |
Family
ID=70557539
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201911327726.5A Active CN111159869B (zh) | 2019-12-20 | 2019-12-20 | 一种砂岩铀矿床与古河道间定量参数的自动计算方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111159869B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113627657A (zh) * | 2021-07-23 | 2021-11-09 | 核工业北京地质研究院 | 一种使用机器学习模型的砂岩型铀成矿有利区预测方法 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104237964A (zh) * | 2014-08-14 | 2014-12-24 | 核工业北京地质研究院 | 古河道型砂岩铀矿定位方法 |
| CN108241176A (zh) * | 2016-12-26 | 2018-07-03 | 核工业北京地质研究院 | 一种放射性综合信息异常远景区圈定方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| FR2868848B1 (fr) * | 2004-04-07 | 2006-06-23 | Earth Decision Sciences Sa | Procede, dispositif et produit-programme de simulation stochastique de paleo chenaux ou de paleo structures geologiques |
-
2019
- 2019-12-20 CN CN201911327726.5A patent/CN111159869B/zh active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104237964A (zh) * | 2014-08-14 | 2014-12-24 | 核工业北京地质研究院 | 古河道型砂岩铀矿定位方法 |
| CN108241176A (zh) * | 2016-12-26 | 2018-07-03 | 核工业北京地质研究院 | 一种放射性综合信息异常远景区圈定方法 |
Non-Patent Citations (2)
| Title |
|---|
| 对航空放射性测量数据进行无分布性分析来识别富铀或贫铀区;Fredric L.Prikle 等;国外铀矿地质;全文 * |
| 砂岩型铀矿弱信息提取及综合评价方法技术研究;韩绍阳 等;世界核地质科学;第21卷(第04期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111159869A (zh) | 2020-05-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109441422B (zh) | 一种页岩气井间距优化开采方法 | |
| CN108562950B (zh) | 一种基于测井信息对地层层位进行智能化划分的方法 | |
| CN105089662B (zh) | 碳酸盐岩地层录井综合柱状图的校正方法及系统 | |
| CN104237964A (zh) | 古河道型砂岩铀矿定位方法 | |
| CN104018828A (zh) | 一种基于演化进程的曲流河砂体储层建筑结构分析方法 | |
| CN106570287A (zh) | 一种基于三维离散裂隙网络预测隧道突涌水量的方法 | |
| CN111159869B (zh) | 一种砂岩铀矿床与古河道间定量参数的自动计算方法 | |
| CN106777707B (zh) | 一种利用改进的蜘蛛网图进行测井岩性定量识别的方法 | |
| CN110850057A (zh) | 一种基于自相似性理论的储层裂缝建模方法及系统 | |
| CN111598340B (zh) | 基于分数阶希尔伯特变换的薄砂体平面展布预测方法 | |
| CN111199107A (zh) | 三角洲相砂岩圈闭评价新方法 | |
| CN116482757A (zh) | 基于浅剖井震联合的地层结构识别分析方法 | |
| CN113534283B (zh) | 一种砂岩型铀矿成矿要素特征的定量评价方法 | |
| CN111999768B (zh) | 一种钻井分层设计中基于构造模型的构造海拔确定方法 | |
| CN111768104B (zh) | 单斜砂泥岩互层型坝基岩性变化程度定量评价方法 | |
| CN110989034B (zh) | 一种回归-分形插值法反演测井横波时差方法 | |
| CN111485868B (zh) | 一种基于开发方案的煤层气田动用储量估算方法 | |
| CN106646608B (zh) | 一种用于观测系统参数论证的层位模型建立方法 | |
| CN112950084A (zh) | 一种反渗透海水淡化厂选址方法 | |
| CN104199107A (zh) | 基于垂直地震多波数据的钻前深度预测方法以及系统 | |
| CN111274736A (zh) | 一种基于有导师学习神经网络算法的导水裂隙带预测方法 | |
| CN117251739B (zh) | 用于砂岩型铀资源定量预测的样本集的构建方法 | |
| CN111175825A (zh) | 深度域速度建模方法 | |
| CN112147713A (zh) | 泥页岩总有机碳含量分段预测方法 | |
| CN110320572B (zh) | 一种识别沉积相的方法及系统 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |