Quantitative Research for Site-selection of Antarctic Year-round Research Stations Based on Logistic Regression
-
摘要: 南极常年考察站是人类在南极地区科考、探索以及观光等活动的据点,其选址对南极自然环境产生一定影响并制约考察站本身功能。本文突破传统考察站选址方法,基于地理信息系统,采用逻辑回归模型模拟已建站情况,使用定量化方法建立指标体系,评价常年站建站适宜性程度。实验过程摒除了人工指标适宜性分类以及专家打分等过程,降低了人为因素带来的不确定性。实验证明,其对已建站情况的拟合精度大于全指标模型,可以更好地模拟常年站的选址情况。Abstract: Antarctic year-round research stations are the foothold for Antarctic expedition, exploration, and sightseeing activities, whose selection affects the natural environment of Antarctica as well as the function and operational efficiency of the stations. Different from traditional site selection methods, this study utilizes the logistic regression, a quantitative method, to build the criteria system and identify their weights and model the construction environment of Antarctica with the aids from geographical information systems. Due to the elimination of artificial data classification and expert scoring process, the logistic model generated in this study largely reduced the uncertainties brought by human decisions. It turned out that the logistic model was more suitable to the construction environment than the model with all criteria, and more precise and detailed, showing the feasibility of the results.
-
利用地磁场方向的稳定性,通过测量传感器轴线与地磁场水平分量间的夹角实现载体航向的准确估计,是地磁定向导航及航位推算导航的基础,因此,基于磁传感器的航向估计技术得到了较为广泛的关注[1, 2]。现有的电子磁罗盘在无地磁异常干扰的情况下定向精度可达到±0.5°,但随着人类活动的加剧,建筑物、电缆及地下掩埋物等多种人工磁源逐渐增多,导致地球浅场磁异常的加剧,使得经过罗差校正的磁罗盘在磁异常干扰条件下的定向误差[3]达到数10°,严重影响了其在导航中的应用。因此,本文针对存在浅场地磁异常情况下的磁航向估计问题进行研究,设计了新的传感器阵列结构,提出了基于磁异常反演的磁航向误差实时补偿方法,以提高磁航向技术在实际应用中的鲁棒性。
1 磁传感器阵列设计及磁航向扰动探测
1.1 磁传感器阵列设计
浅场磁异常叠加在背景地磁场上,要实现精确的磁航向估计,必须对浅场磁异常进行有效分离。基于磁传感器测量得到的标量及矢量数据均为两种场的叠加,无法实现较好的分离[4],而背景地磁场在小范围内可认为是匀强磁场,即地磁梯度场较小且可忽略不计[5, 6],因此,利用传感器测量得到的磁梯度场可认为仅是浅场磁干扰源产生的异常梯度场,进而可实现背景地磁场和浅场磁异常在梯度测量上的分离。但由于磁航向估计技术是基于磁矢量信息的,必须利用浅场磁异常的梯度场反演得到其产生的磁矢量场,因此,为实现浅场磁异常梯度场的测量和矢量场的准确反演,本文设计了如图 1所示的磁传感器阵列结构。该传感器阵列由5个三轴磁通门传感器构成,相同轴上两个传感器间的基线距离为d,外围的4个传感器可构建平面十字形磁梯度张量系统用于测量多个不同方向的梯度场,进而得到磁梯度张量矩阵用于浅场磁异常源的反演。
由于地磁场及浅场异常场可看作无源的静磁场,因此,磁感应强度的散度和旋度为零,磁梯度张量的9个分量中仅有5个分量是独立的。建立如图 1所示的右手坐标系用于磁场建模,则由外围的4个传感器测量得到的磁梯度张量矩阵为:
(1) 式中,Bij表示第i个传感器在j方向上的磁场矢量值,i=1,2,3,4,j=x,y,z;Bpq表示不同方向的磁场梯度值,p,q=x,y,z。
值得注意的是,在传感器阵列构建并将其安装在载体上的过程中,传感器阵列自身的误差以及载体中铁磁材料引起的磁干扰将会严重影响磁航向的估计,在实际使用前,必须将此两种磁干扰有效消除和补偿。而本文主要针对已经完成上述两种磁干扰补偿的系统在实际应用中受到外部浅场磁异常目标干扰进而影响磁航向估计的问题进行研究,利用磁异常反演的方法补偿浅场磁异常引起的磁航向估计误差,最终实现较为准确的磁航向估计。
1.2 磁航向扰动探测
载体航向估计时,磁传感器与载体捷联在一起,故载体坐标系和传感器坐标系一致,本文建立的右手坐标系如图 2所示。假设无浅场磁异常存在时,三轴磁传感器测得的水平方向上的两个磁场分量分别为Hy和Hx,测量点处地磁场的磁偏角、磁倾角及总场值分别为D、I和F,则载体的磁航向可由式(2)估计得到:
(2) 由图 2可知,浅场磁源的存在将影响地磁场水平方向的磁矢量扰动,进而引起磁航向估计的误差,但并不是所有的浅场磁异常都会改变磁航向,且磁异常较小的情况下也可能引起较大的磁航向误差,因此,必须构建磁航向扰动探测器,准确判断磁航向变化是由于载体自身的运动扰动还是由于浅场磁异常引起的。
浅场磁异常源引起不同方向上梯度测量值的变化,原则上可根据传感器阵列测量得到的梯度值判断是否有磁异常源的存在,但由于载体在运动过程中姿态的变化对测得的梯度值具有一定的影响[7],因此,直接采用梯度值进行磁异常源是否存在的判断具有一定的非鲁棒性。磁梯度张量的张量不变量是由不同方向的梯度分量组合得到的,不受载体姿态变化的影响,因此,本文将其用于磁异常源是否存在的判断。
磁梯度张量存在多个张量不变量可用于磁异常探测,考虑磁航向估计对实时性的要求,选择I1和I2作为判断因子,其数学表达式分别为:
(3) 由式(3)可知,相比于I1,不变量I2随着距离的增大衰减更快,对于同时存在多个磁异常源的复杂情况,其优先描述浅部源的特性。
若磁传感器阵列测得的地磁场不存在浅场磁异常,则可直接利用式(2)进行磁航向的计算,若存在磁异常,需要进一步判断该异常是否引起磁航向的扰动,比如水平方向的两个磁矢量值等比例放大或缩小将不影响磁航向的计算。利用零号传感器测得的地磁场7要素中的总场F、水平分量H、竖直分量Hz和磁倾角I进行磁异常是否影响磁航向的联合判断。若磁异常对磁航向没有干扰或干扰极小,则直接利用零号传感器的测量值进行磁航向计算,若磁异常对磁航向的计算产生干扰,则利用传感器阵列测得的磁梯度张量值进行磁异常的反演,利用反演得到的传感器阵列中心点O位置处的磁矢量值补偿零号传感器测得的磁矢量值,进而在背景地磁场准确测量的情况下实现磁航向的准确估计。磁航向扰动探测流程如图 3所示。
2 磁异常反演及航向误差实时补偿
由于磁航向估计利用磁矢量信息进行载体运动方向的估计,因此,需要利用传感器阵列测量得到的浅层磁异常源的磁梯度张量数据反演其在测量中心点产生的磁矢量值,然后将零号传感器的测量值减去磁异常源产生的磁矢量值,即为背景地磁场的磁矢量值,进而可实现较为精确的磁航向估计。
在磁异常反演中,当探测距离大于2.5倍磁性物体长度时,可将磁性目标简化为磁偶极子模型,并由6个未知量描述,即磁性目标的三维位置和三维磁矩[8]。因此,要实现磁航向估计的实时误差补偿,必须实现磁异常的实时反演,即在单个测点上实时地完成磁偶极子位置及磁矩的反演。
目前,基于单点磁梯度张量数据的磁偶极子定位方法受到了较多关注,其中,Nara[9, 10]提出了联合单点测量的磁梯度张量及磁矢量数据进行磁偶极子定位的方法,但在地磁背景下无法获得磁偶极子的磁矢量信息;Roy[11]为克服背景地磁场的影响,提出了基于张量不变量的单点定位算法,但该算法需要较为复杂的磁梯度张量测量系统。Frahm[12]基于单点磁梯度张量数据的特征分析,实现了磁偶极子单位位置向量及单位磁矩向量的计算,但该方法估计得到4个不同的解且无法实现磁偶极子磁矢量信息的计算,本文在此方法基础上,联合位置、磁矩及磁梯度张量系统中包含的多点磁标量信息构建非线性方程组,在消除三个虚假解的同时,解算磁偶极子的磁矩模及测量点与磁性目标间的距离,进而实现磁矢量信息的正确反演。
假设Frahm方法得到的4个解的单位位置向量及磁矩向量为 ${{\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\frown}$}}{r}}_{k}} {{\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\frown}$}}{m}}_{k}},k=1,2,3,4$ 。若测量点与磁偶极子间的真实距离为R,磁偶极子的真实磁矩模为M,则磁偶极子的磁矩矢量及测点相对于磁偶极子的位置矢量分别为 ${{m}_{k}}=M\centerdot {{\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\frown}$}}{m}}_{k}}=\left( {{m}_{kx}},{{m}_{ky}},{{m}_{kz}} \right)$ 和 $\begin{align} & {{m}_{k}}=M\centerdot {{{\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\frown}$}}{m}}}_{k}}=\left( {{m}_{kx}},{{m}_{ky}},{{m}_{kz}} \right) \\ & {{r}_{k}}=R\centerdot {{{\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\frown}$}}{r}}}_{k}}=\left( {{r}_{kx}},{{r}_{ky}},{{r}_{kz}} \right) \\ \end{align}$ 。
由传感器的阵列结构可知,针对不同的解,外围的4个传感器相对于磁偶极子的位置矢量可分别表示为:
(4) 则传感器与磁偶极子之间的距离可表示为Rkl=‖rkl‖,l=1,2,3,4。
设4个传感器测得的磁场总场值分别为F1ce,F2ce,F3ce和F4ce,则由磁偶极子产生的磁场总场正演公式可得以下非线性方程组:
(5) 式中,仅存在R、M和F 3个未知变量,求解方程即可得到磁偶极子的位置和磁矩矢量 rk、mk。
根据磁偶极子的磁梯度张量正演公式,计算4个不同的磁偶极子在测量点产生的磁梯度张量值,并将其与真实值相比较即可实现虚假解的去除和磁性目标的单点磁梯度张量定位及磁矩反演,进而完成磁异常的三分量磁矢量反演计算。
假设反演得到磁偶极子的磁矩为m0,相对于测量中心点的位置为r0,则其在测量中心点O处产生的三分量磁矢量值为:
(6) 补偿得到的背景地磁场的三分量值为:
(7) 式中,B0为零号传感器测得的三分量磁矢量值。
则在存在影响磁航向估计的磁异常场的情况下,真实磁航向可利用式(7)中反演得到的背景地磁场的三分量值进行计算。
3 实验与分析
为验证本文提出的磁航向误差补偿方法的可行性,采用5个三轴磁通门传感器构建传感器阵列,在地面上放置永磁体模拟浅场磁异常进行磁航向实验,如图 4所示。为详细分析磁航向估计中磁异常反演的正确性,实验采用单点测量的方式进行,实验人员手持传感器阵列且距离地面约50 cm,在测量曲线上每隔10 cm进行一次数据采集。
在永磁体左右两侧各选取一条测线进行磁航向的计算,测量得到的张量不变量及补偿前后的方位角如图 5、图 6所示。
由图 5、 图 6可知,磁梯度张量的张量不变量I1和I2可实现磁异常源是否存在的准确判断。实验所用永磁体的磁矢量场与地磁矢量场间存在较大的方向偏差,因此,在传感器阵列距离永磁体较近时,产生了10~20°左右的磁航向估计误差。在利用磁传感器阵列测得的磁梯度张量数据进行磁异常源产生的磁矢量场反演的基础上,进行了磁航向计算的补偿,从补偿结果可以看出,论文所提算法可实现浅场磁异常源引起的磁航向误差的补偿。补偿后得到的磁航向方位角仍存在一定的波动,这是由于实验人员手持传感器阵列无法一直保持水平姿态及方向不变而导致的。
4 结 语
本文针对由地磁矢量信息进行磁航向估计时易受浅场磁异常源干扰的问题,设计了传感器阵列结构,通过测量磁异常源产生的磁梯度场进而反演其产生的磁异常矢量场,在此基础上计算得到较为准确的背景地磁矢量场,实现磁航向误差的实时补偿。实验结果表明,利用磁梯度张量的单点反演方法可有效消除磁异常目标引起的磁航向误差,具有较高的补偿精度。
-
表 1 逻辑回归的全指标模型及最优模型
Table 1 Model Including all Criteria and the Best Model Based on Logistic Regression
全指标模型 最优模型 b Sig Exp (b) b Sig Exp (b) X1 -0.001 0.01 0.999 / / / X2 -0.583 0.00 0.558 -0.666 0.00 0.514 X3 -0.002 0.00 0.998 -0.002 0.00 0.998 X4 -0.007 0.01 0.993 -0.006 0.02 0.994 X5 0 0.49 1 / / / X6 0.003 0.03 1.003 0.004 0.00 1.004 X7 0 0.93 1 -0.002 0.00 0.998 X8 -0.006 0.02 0.994 / / / X9 -0.006 0.03 0.994 / / / X10 0.002 0.06 1.002 0.001 0.02 1.001 X11 -0.008 0.00 0.992 -0.008 0.00 0.992 X12 0.007 0.00 1.007 / / / X13 -0.003 0.04 0.997 / / / X14 -0.002 0.06 0.998 / / / X15 0 0.85 1 / / / X16 0 0.81 1 / / / X17 0 0.82 1 -0.001 0.00 0.999 X18 0.439 0.03 1.551 / / / X19 0.576 0.00 1.779 0.712 0.00 2.038 X20 -0.076 0.18 0.927 / / / X21 -0.045 0.27 0.956 / / / a -10.309 0.00 0 -7.102 0.00 0.001 表 2 Hosmer和Lemeshow检验结果
Table 2 Result of Hosmer and Lemeshow Test
全指标模型 最优模型 卡方检验 14.568 9.394 df (自由度) 8 8 Sig (P) 0.068 0.310 -
[1] 陈杰, Blume H P.人类活动对南极陆地生态系统的影响[J].极地研究, 2000, 12(1):62-75 http://www.cnki.com.cn/Article/CJFDTOTAL-JDYZ200001006.htm Chen Jie, Blume H P. Anthropic Impacts on the Antarctic Terrestrial Ecosystem[J]. Chinese Journal of Polar Research, 2000, 12(1):62-75 http://www.cnki.com.cn/Article/CJFDTOTAL-JDYZ200001006.htm
[2] Belgian Federal Science Policy Office. The International Polar Foundation. Final Comprehensive Environmental Evaluation: Construction and Operation of the New Belgian Research Station, Dronning Maud Land, Antarctica[C]. The 28th Antarctic Treaty Consultative Meeting, Stockholm, Sweden, 2005
[3] British Antarctic Survey. Final Comprehensive Environmental Evaluation: Proposed Construction and Operation of Halley Ⅵ Research Station, and Demolition and Removal of Halley Ⅴ Research Station, Brunt ice Shelf, Antarctica[C]. The 28th Antarctic Treaty Consultative Meeting, Stockholm, Sweden, 2005
[4] Korea Polar Research Institute, Korea Environment Institute. Final Comprehensive Environmental Evaluation: Construction and Operation of the Jang Bogo Antarctic Research Station, Terra Nova Bay, Antarctica[C]. The 35th Antarctic Treaty Consultative Meeting, Hobart, Australia, 2012
[5] Pang Xiaoping, Liu Haiyan, Zhao Xi.Selecting Suitable Sites for an Antarctic Research Station: A Case for a New Chinese Research Station[J]. Antarctic Science, 2014, 26(5): 479-490 doi: 10.1017/S0954102013000965
[6] 刘海燕, 庞小平.基于GIS和模糊层次分析法的南极考察站选址研究[J].武汉大学学报·信息科学版, 2015, 40(2):249-252 http://ch.whu.edu.cn/CN/abstract/abstract3193.shtml Liu Haiyan, Pang Xiaoping. Antarctic Research Station Site Selection Based on GIS and Fuzzy AHP[J]. Geomatics and Information Science of Wuhan University, 2015, 40(2):249-252 http://ch.whu.edu.cn/CN/abstract/abstract3193.shtml
[7] Ho Zhiyong, Lo C P. Modeling Urban Growth in Atlanta Using Logistic Regression[J]. Computer Environment and Urban Systems, 2007, 31: 667-688 doi: 10.1016/j.compenvurbsys.2006.11.001
[8] Menard S.应用Logistic回归分析[M].上海:格致出版社, 2009 Menard S. Applied Logistic Regression Analysis[M]. Shanghai: True & Wisdom Press, 2009
[9] 王济川, 郭志刚. Logistic回归模型—方法与应用[M].北京:高等教育出版社, 2001 Wang Jichuan, Guo Zhigang. Logistic Regression Model-Method and Application[M]. Beijing: Higher Education Press, 2001
[10] 姜文亮, 张晓通, 李霖, 等.基于GIS和空间逻辑回归模型的城市空间拓展预测[J].测绘科学, 2008, 33(4):172-174 http://www.cnki.com.cn/Article/CJFDTOTAL-CHKD200804064.htm Jiang Wenliang, Zhang Xiaotong, Li Lin, et al. Urban Spatial Expansion Forecast Based on GIS & Spatial Logistic Regression Model[J]. Science of Surveying and Mapping, 2008, 33(4):172-174 http://www.cnki.com.cn/Article/CJFDTOTAL-CHKD200804064.htm
[11] 凌晓良, 陈丹红, 张侠, 等.南极特别保护区的现状与展望[J].极地研究, 2008, 20(1):48-63 http://www.cnki.com.cn/Article/CJFDTOTAL-JDYZ200801009.htm Ling Xiaoliang, Chen Danhong, Zhang Xia, et al. Review of the Status, Feature and Prospect of Antarctic Specially Protected Areas[J]. Chinese Journal of Polar Research, 2008, 20(1):48-63 http://www.cnki.com.cn/Article/CJFDTOTAL-JDYZ200801009.htm
[12] 顾悦婷, 孙波, 陈丹红, 等.南极特别管理区现状分析与未来展望[J].极地研究, 2010, 22(4):431-440 http://www.cnki.com.cn/Article/CJFDTOTAL-JDYZ201004012.htm Gu Yueting, Sun Bo, Chen Danhong, et al. Review of the Current Status, Feature and Future Development of Antarctic Specially Managed Area[J]. Chinese Journal of Polar Research, 2010, 22(4):431-440 http://www.cnki.com.cn/Article/CJFDTOTAL-JDYZ201004012.htm
-
期刊类型引用(3)
1. 汪金花,张博,郭立稳,刘暑明,张恒嘉. 井下巷道地磁匹配特征的CEA卷积增强的分析. 武汉大学学报(信息科学版). 2022(09): 1422-1431 . 百度学术
2. 贾文抖,林春生,林朋飞,翟国君. 基于正弦拟合的三轴磁力仪标定方法. 武汉大学学报(信息科学版). 2019(09): 1328-1333+1406 . 百度学术
3. 申家双,葛忠孝,陈长林. 我国海洋测绘研究进展. 海洋测绘. 2018(04): 1-10+21 . 百度学术
其他类型引用(4)