基于神經(jīng)網(wǎng)絡(luò)方法獲得最優(yōu)化月球內(nèi)部結(jié)構(gòu)模型

廖彬彬,徐建橋，陳曉東*,孫和平,周江存

1 中國(guó)科學(xué)院精密測(cè)量科學(xué)與技術(shù)創(chuàng)新研究院 大地測(cè)量與地球動(dòng)力學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室，武漢 430077 2 中國(guó)科學(xué)院大學(xué)，北京 100049

0 引言

月球是離地球最近的天體，是人類邁向太空的第一站.研究月球內(nèi)部結(jié)構(gòu)對(duì)人類進(jìn)行空間探索、了解地月系統(tǒng)的演化、開(kāi)采并利用月球資源等方面具有重要意義(Smith et al.,2017;Zhao et al.,2018；Sun et al.,2019).20世紀(jì)阿波羅計(jì)劃在月球表面開(kāi)展了一系列月震觀測(cè)實(shí)驗(yàn)，這些月震觀測(cè)數(shù)據(jù)使人類對(duì)月球內(nèi)部結(jié)構(gòu)有了初步的認(rèn)識(shí).早期月震的分析結(jié)果表明月球內(nèi)部主要存在以下幾個(gè)構(gòu)造圈層：一個(gè)厚度約為50～70 km的月殼、一個(gè)厚度較大的月幔且在月幔深處存在一個(gè)低速層(Nakamura et al.,1974;Goins et al.,1981).后來(lái)月球極移和月球磁場(chǎng)的觀測(cè)結(jié)果證明了月球存在一個(gè)液態(tài)的外核(Williams et al.,2001；Shimizu et al.,2013).在阿波羅計(jì)劃之后，高分辨率遙感衛(wèi)星技術(shù)的應(yīng)用以及近年來(lái)對(duì)早期月震數(shù)據(jù)的重新處理使得對(duì)月球結(jié)構(gòu)的認(rèn)識(shí)有了新的突破(Hao et al.,2018).Sugano 和 Heki(2004)利用LP(Lunar Prospector，月球勘探者)衛(wèi)星的重力觀測(cè)數(shù)據(jù)推斷月殼厚度約為20～60 km.Khan等(2000)利用非線性貝葉斯方法對(duì)阿波羅月震數(shù)據(jù)進(jìn)行了重新處理，結(jié)果表明在月殼內(nèi)深度為20 km以及45±5 km處存在速度間斷面，并且月幔內(nèi)深度為500 km處P波速度與S波速度顯著增加.Lognonné等(2003)利用體波走時(shí)數(shù)據(jù)通過(guò)隨機(jī)搜索方法發(fā)現(xiàn)月幔內(nèi)深度500 km以及750 km處均存在速度間斷面.對(duì)于月球更深部分的分層結(jié)構(gòu)，Konopliv等(1998)基于LP衛(wèi)星得到的月球質(zhì)量、平均轉(zhuǎn)動(dòng)慣量以及重力數(shù)據(jù)推測(cè)：如果月核完全由固態(tài)鐵組成，則半徑為220 km到370 km之間；而如果月核為液態(tài)的Fe-FeS共晶體，則半徑為330～590 km之間.Garcia等(2011)基于大地測(cè)量觀測(cè)數(shù)據(jù)(包括月球質(zhì)量、轉(zhuǎn)動(dòng)慣量以及勒夫數(shù))以及月震數(shù)據(jù)得到了月球結(jié)構(gòu)模型VPREMOON，其中的月核半徑為380±40 km，月核平均密度為5200±1000 kg·m-3.Yan等(2015)通過(guò)Monte Carlo方法得到月核半徑為370 km.而Weber等(2011)通過(guò)反演得到月核由固態(tài)內(nèi)核以及液態(tài)外核組成，其中固態(tài)內(nèi)核半徑為240 km，液態(tài)外核半徑范圍為240～330 km．Williams等 (2014)以Weber等(2011)月球內(nèi)部結(jié)構(gòu)模型為初始模型，利用GRAIL(Gravity Recovery and Interior Laboratory,重力重建與內(nèi)部結(jié)構(gòu)實(shí)驗(yàn)室)重力場(chǎng)數(shù)據(jù)計(jì)算了月球2階固體潮勒夫數(shù)、平均密度以及平均轉(zhuǎn)動(dòng)慣量等物理量，通過(guò)調(diào)整低速區(qū)彈性參數(shù)、低速區(qū)厚度以及內(nèi)外核厚度，得到了滿足GRAIL觀測(cè)數(shù)據(jù)的5個(gè)月球內(nèi)部結(jié)構(gòu)模型 (GPM 1—5).然而Matsuyama等(2016)認(rèn)為Williams等(2014)采用的勒夫數(shù)k2=0.02422±0.000012是粘彈性月球的引潮位響應(yīng)，經(jīng)過(guò)改正后給出了彈性月球的勒夫數(shù)為k2=0.0220±0.0017.并且Matsuyama等(2016)利用Markov Chain-Monte Carlo方法進(jìn)行反演，結(jié)果顯示月幔底部低速區(qū)存在的概率約為50%.綜上所述，從上面目前國(guó)內(nèi)外對(duì)月球結(jié)構(gòu)研究結(jié)果可以看出，現(xiàn)階段對(duì)月球結(jié)構(gòu)模型的認(rèn)識(shí)，特別是月球深部結(jié)構(gòu)的認(rèn)識(shí)，仍然存在許多的不確定性，這也是至今沒(méi)有一個(gè)月球內(nèi)部結(jié)構(gòu)模型被廣泛認(rèn)可和采用的根本原因.

目前可以用來(lái)反演月球內(nèi)部結(jié)構(gòu)的數(shù)學(xué)方法有很多，各種方法或多或少都存在一定的局限性.傳統(tǒng)線性反演方法要求模型參數(shù)與觀測(cè)數(shù)據(jù)之間的映射是(擬)線性的，然而實(shí)際中觀測(cè)數(shù)據(jù)與模型的關(guān)系往往是非線性的，甚至可能都不是一一對(duì)應(yīng)的映射關(guān)系(Tarantola,2005).已有的研究結(jié)果表明利用非線性Bayes反演方法可以很好解決傳統(tǒng)線性反演方法的不足(De Wit et al.,2013).Bayes反演方法基于Bayes理論，利用觀測(cè)數(shù)據(jù)d、模型的先驗(yàn)分布σ(m)和似然函數(shù)L(d|m)計(jì)算模型參數(shù)后驗(yàn)概率密度分布σ(m|d)，從而獲得模型的最優(yōu)解，即最大后驗(yàn)概率對(duì)應(yīng)的模型(De Wit et al.,2013).由于采用了更加復(fù)雜的概率模型與參數(shù)估計(jì)方法，Bayes反演方法完全擺脫最小二乘方法要求觀測(cè)數(shù)據(jù)和模型參數(shù)之間存在線性關(guān)系的限制.Bayes反演方法通常利用Monte Carlo方法來(lái)實(shí)現(xiàn)(K?ufl et al.,2013)，然而Monte Carlo方法受限于“維度詛咒”(Bellman,1961)，往往計(jì)算量巨大，需要較長(zhǎng)的計(jì)算時(shí)間.另一種可行的Bayes反演方法是混合密度神經(jīng)網(wǎng)絡(luò)方法(Mixture Density Neural Network，MDN)，其基本原理是建立一個(gè)以觀測(cè)數(shù)據(jù)為輸入，以模型參數(shù)的后驗(yàn)概率密度分布為輸出的神經(jīng)網(wǎng)絡(luò)(Bishop,1995).與Monte Carlo方法相比，MDN方法可以避免“維度詛咒”，從而極大地提高計(jì)算效率(De Wit et al.，2013).此外MDN方法還可以得到觀測(cè)數(shù)據(jù)與模型參數(shù)之間的函數(shù)關(guān)系，為進(jìn)一步分析月球內(nèi)部結(jié)構(gòu)模型的精度提供了可能.

本研究采用Williams等(2014)基于GRAIL觀測(cè)任務(wù)給出的平均密度ρ和平均轉(zhuǎn)動(dòng)慣量MOI以及Matsuyama等(2016)給出的彈性勒夫數(shù)k2作為觀測(cè)約束，通過(guò)MDN方法獲得了月球內(nèi)部結(jié)構(gòu)模型的低速層、外內(nèi)核的半徑、P波和S波速度等8個(gè)模型參數(shù)的后驗(yàn)概率密度分布，計(jì)算了平均月球內(nèi)部結(jié)構(gòu)模型Mean、最大后驗(yàn)概率密度對(duì)應(yīng)的MAP模型以及滿足1-σ準(zhǔn)則的月球內(nèi)部結(jié)構(gòu)模型的參數(shù)范圍.最后利用觀測(cè)值與月球內(nèi)部結(jié)構(gòu)模型的參數(shù)之間的函數(shù)關(guān)系，采用不同觀測(cè)值觀測(cè)誤差范圍內(nèi)月球內(nèi)部結(jié)構(gòu)模型參數(shù)的變化幅度作為度量，定量地分析了不同觀測(cè)值的精度對(duì)月球內(nèi)部結(jié)構(gòu)模型的影響.

1 數(shù)據(jù)和方法

對(duì)月球內(nèi)部結(jié)構(gòu)的反演工作主要由兩部分組成：(1)以Weber等(2011)的月球內(nèi)部結(jié)構(gòu)模型為初始模型，對(duì)內(nèi)外核及低速區(qū)的厚度、P波及S波速度進(jìn)行均勻采樣，得到所有可能的月球內(nèi)部結(jié)構(gòu)模型，然后正演計(jì)算得到對(duì)應(yīng)的固體潮勒夫數(shù)、平均密度及平均轉(zhuǎn)動(dòng)慣量；(2)利用混合密度神經(jīng)網(wǎng)絡(luò)方法和上面計(jì)算的三個(gè)觀測(cè)值，得到對(duì)應(yīng)月球內(nèi)部結(jié)構(gòu)模型的后驗(yàn)概率密度.

1.1 觀測(cè)數(shù)據(jù)的正演計(jì)算

假設(shè)月球是一個(gè)非自轉(zhuǎn)、完全彈性、各向同性以及橫向均勻的分層球體，并且各個(gè)圈層的密度以及彈性參數(shù)均為常數(shù).Matsumoto等(2015)利用Bayes方法反演了月球內(nèi)部結(jié)構(gòu)，他們的結(jié)果表明：月幔底部低速區(qū)以及月核的半徑、彈性參數(shù)和密度都存在很大的不確定性，而低速區(qū)以上的部分模型參數(shù)的不確定性相對(duì)較小.因此本研究假設(shè)在低速區(qū)以上部分月球內(nèi)部結(jié)構(gòu)模型參數(shù)為定值，均采用Weber等(2011)給出月球內(nèi)部結(jié)構(gòu)模型中的數(shù)值，而低速區(qū)以下的模型參數(shù)則為反演的自由度.此外由于觀測(cè)值個(gè)數(shù)較少，出于簡(jiǎn)化反演問(wèn)題的考慮，低速區(qū)與內(nèi)外核的密度同樣也采用Weber等(2011)給出月球內(nèi)部結(jié)構(gòu)模型中的對(duì)應(yīng)數(shù)值.并且為了與嫦娥1號(hào)、Clementine等其他月球任務(wù)給出的月球半徑觀測(cè)數(shù)據(jù)吻合(Ping et al.,2009;Zuber et al.,2013;Smith et al.,2017)，Williams等(2014)建議將Weber等(2011)的月球內(nèi)部結(jié)構(gòu)模型中的半徑調(diào)整為R=1737.151 km.最終在反演問(wèn)題中需要考慮的模型參數(shù)總共8個(gè)，包括內(nèi)核半徑ri、內(nèi)核P波速度VPi、內(nèi)核S波速度VSi、外核半徑ro、外核P波速度VPo、低速區(qū)半徑rl、低速區(qū)P波速度VPl以及低速區(qū)S波速度VSl，所有參數(shù)先驗(yàn)的變化范圍如表2中“先驗(yàn)?zāi)Ｐ头秶币涣兴?假設(shè)8個(gè)模型參數(shù)的先驗(yàn)概率分布滿足均勻分布，在表2給出的先驗(yàn)?zāi)Ｐ头秶鷥?nèi)對(duì)它們進(jìn)行隨機(jī)采樣，總共生成100000個(gè)月球內(nèi)部結(jié)構(gòu)模型樣本.

研究中觀測(cè)值固體潮勒夫數(shù)k2、平均密度ρ以及平均轉(zhuǎn)動(dòng)慣量MOI的具體數(shù)值以及誤差范圍如表1所示.本文中固體潮勒夫數(shù)k2采用的是Matsuyama等(2016)給出的彈性勒夫數(shù)k2=0.0220±0.0017.為了獲得上述100000個(gè)月球內(nèi)部結(jié)構(gòu)模型樣本的模擬觀測(cè)值，首先通過(guò)正演數(shù)值計(jì)算方法獲得每個(gè)月球內(nèi)部結(jié)構(gòu)模型對(duì)應(yīng)固體潮勒夫數(shù)、平均密度以及平均轉(zhuǎn)動(dòng)慣量的模擬觀測(cè)值.其中在求解二階潮汐勒夫數(shù)時(shí)，需要求解一個(gè)由動(dòng)量平衡方程與泊松方程組成的二階偏微分方程組邊值問(wèn)題，本研究采用廖彬彬等(2019)給出的譜元法進(jìn)行數(shù)值計(jì)算.對(duì)于平均密度和平均轉(zhuǎn)動(dòng)慣量，則用Gauss-Legendre數(shù)值積分方法進(jìn)行計(jì)算(Quarteroni，2009).此外由于實(shí)際觀測(cè)值存在誤差，因此需要對(duì)上面正演計(jì)算得到的模擬觀測(cè)值加上觀測(cè)誤差才能得到實(shí)際觀測(cè)值.研究中假設(shè)三個(gè)觀測(cè)值在誤差范圍內(nèi)服從均勻分布，對(duì)其進(jìn)行隨機(jī)采樣得到相應(yīng)的觀測(cè)誤差，然后加上對(duì)應(yīng)的觀測(cè)值便可得到之后訓(xùn)練神經(jīng)網(wǎng)絡(luò)所需要的模擬觀測(cè)值.

表1 反演計(jì)算采用的觀測(cè)值及其精度Table 1 Observations and their accuracies used in the inversion calculation

1.2 用混合密度神經(jīng)網(wǎng)絡(luò)方法反演后驗(yàn)概率密度

神經(jīng)網(wǎng)絡(luò)反演方法主要特點(diǎn)在于它可以解決反演過(guò)程中的“非線性”和“不唯一性”等問(wèn)題.該方法與一般采樣反演方法相比，計(jì)算效率更高，關(guān)鍵是還可得到模型參數(shù)后驗(yàn)概率密度分布與觀測(cè)值間的映射關(guān)系.這種映射關(guān)系可用來(lái)進(jìn)一步研究觀測(cè)值精度對(duì)模型精度的影響，也是本文研究模擬觀測(cè)值精度對(duì)月球內(nèi)部結(jié)構(gòu)模型精度影響采用的主要數(shù)學(xué)方法.神經(jīng)網(wǎng)絡(luò)的數(shù)學(xué)本質(zhì)是實(shí)現(xiàn)從輸入層(Input layer)到輸出層(Output layer)之間任意映射關(guān)系的近似表達(dá)，其映射關(guān)系如圖1所示，可用下面公式(1)和(2)來(lái)描述(De Wit et al.,2013)：

圖1 神經(jīng)網(wǎng)絡(luò)映射關(guān)系示意圖Fig.1 Schematic diagram of the neural network mapping relationship

(1)

(2)

混合密度神經(jīng)網(wǎng)絡(luò)是以觀測(cè)數(shù)據(jù)為輸入層，以其反演的結(jié)構(gòu)模型后驗(yàn)概率密度函數(shù)為輸出層的兩層神經(jīng)網(wǎng)絡(luò)(Bishop,2006).MDN方法輸出層的后驗(yàn)概率密度函數(shù)為高斯混合分布，其定義如下：

(3)

式中w表示公式(1)和(2)中的權(quán)重和偏置，αi是Softmax函數(shù)，它代表第i個(gè)高斯核函數(shù)φi的權(quán)重.高斯核函數(shù)φi定義為

φi(m|d;w)=

(4)

式中c為模型參數(shù)m的維度，μi和σi分別為第i個(gè)高斯核函數(shù)的均值及標(biāo)準(zhǔn)差，μi、σi和αi可以分別控制第i個(gè)高斯核函數(shù)的中心位置、寬窄、高低(權(quán)重).不同于傳統(tǒng)線性反演方法中以單個(gè)高斯分布模型作為概率模型，本研究月球內(nèi)部結(jié)構(gòu)模型參數(shù)所用的概率模型(見(jiàn)公式(3))是M個(gè)高斯分布的線性組合.理論上隨著高斯核函數(shù)個(gè)數(shù)M趨近于無(wú)窮大，公式(3)給出的概率模型可以近似地表示任意概率分布函數(shù)(Bishops,1995).利用公式(3)給出的概率模型，可實(shí)現(xiàn)從一個(gè)觀測(cè)值到多個(gè)可能模型值(平均值μi)之間的映射關(guān)系，從而解決了反演的“不唯一性”問(wèn)題.此外通過(guò)混合密度神經(jīng)網(wǎng)絡(luò)方法，還實(shí)現(xiàn)了由觀測(cè)數(shù)據(jù)到概率模型參數(shù)μi、σi和αi的函數(shù)映射，而不僅限于線性映射，因此該方法解決了傳統(tǒng)最小二乘方法只能處理觀測(cè)值與模型參數(shù)為線性映射的問(wèn)題，這也是本研究選用該方法的主要原因.

在使用MDN反演方法的過(guò)程中，構(gòu)建了一個(gè)以固體潮勒夫數(shù)、平均密度及平均轉(zhuǎn)動(dòng)慣量的觀測(cè)值為輸入層，以低速層半徑、內(nèi)外核的厚度、P波和S波速度的聯(lián)合后驗(yàn)概率分布(即公式(3))的參數(shù)μi、σi和αi為輸出層的兩層神經(jīng)網(wǎng)絡(luò).該神經(jīng)網(wǎng)絡(luò)輸入層節(jié)點(diǎn)個(gè)數(shù)為I=3，其激活函數(shù)為“Rule”函數(shù)；隱藏層節(jié)點(diǎn)個(gè)數(shù)為50，其激活函數(shù)為g(a)=a；而對(duì)于輸出層，設(shè)高斯核個(gè)數(shù)M為10，節(jié)點(diǎn)個(gè)數(shù)為K=170.各項(xiàng)超參數(shù)(訓(xùn)練步數(shù)、批訓(xùn)練大小、優(yōu)化器類型、激發(fā)函數(shù)類型以及節(jié)點(diǎn)個(gè)數(shù)等)的具體選擇方法將在下面進(jìn)行敘述.

神經(jīng)網(wǎng)絡(luò)的訓(xùn)練過(guò)程是一個(gè)不斷更新節(jié)點(diǎn)權(quán)值和每一層的偏置(即公式(1)和(2)中的w與b)，從而使得損失函數(shù)接近最小的過(guò)程.MDN方法的損失函數(shù)loss是后驗(yàn)分布函數(shù)(公式(3))的似然函數(shù)(Bishop，1995)，它可以用來(lái)衡量所輸出的概率模型與訓(xùn)練數(shù)據(jù)之間的相似程度，其定義如下：

(5)

其中N為訓(xùn)練樣本總數(shù).訓(xùn)練所采用的最優(yōu)化方法是隨機(jī)梯度下降法，通過(guò)這種方法可將訓(xùn)練集數(shù)據(jù)分批逐次地代入損失函數(shù)當(dāng)中，而不必一次性計(jì)算所有訓(xùn)練集數(shù)據(jù)損失函數(shù)的和(Goodfellow et al.,2016).本研究采用的隨機(jī)梯度下降算法是“Adam”算法.為了避免最優(yōu)化解陷入局部最優(yōu)解，需要對(duì)神經(jīng)網(wǎng)絡(luò)進(jìn)行多次訓(xùn)練，從中選擇損失函數(shù)最小的解，并且在每次訓(xùn)練之前需要對(duì)權(quán)值進(jìn)行不同的初始化(Bishop,1995).此外，為了避免輸出的月球內(nèi)部結(jié)構(gòu)模型中各層半徑出現(xiàn)負(fù)值，本文采用log(ri)、log(ro-ri)和log(rl-ro)代替內(nèi)核半徑ri、外核半徑ro以及低速區(qū)半徑rl作為神經(jīng)網(wǎng)絡(luò)的輸出.最后通過(guò)z-socre標(biāo)準(zhǔn)化方法對(duì)觀測(cè)值與模型參數(shù)進(jìn)行預(yù)處理可以有效避免訓(xùn)練過(guò)程中出現(xiàn)“梯度飽和”(De Wit et al.,2013).

在神經(jīng)網(wǎng)絡(luò)訓(xùn)練完后，還要對(duì)其泛化能力進(jìn)行評(píng)估，神經(jīng)網(wǎng)絡(luò)的泛化能力是指其能正確預(yù)測(cè)未參與訓(xùn)練數(shù)據(jù)的能力.為了度量神經(jīng)網(wǎng)絡(luò)的泛化能力，本研究以8∶1∶1的比例將100000組觀測(cè)數(shù)據(jù)與對(duì)應(yīng)模型參數(shù)成對(duì)組成的數(shù)據(jù)集分成訓(xùn)練集、測(cè)試集以及驗(yàn)證集，三者中只有訓(xùn)練集參與神經(jīng)網(wǎng)絡(luò)的訓(xùn)練，后面兩者作為檢驗(yàn)標(biāo)準(zhǔn)分別在訓(xùn)練過(guò)程中超參數(shù)的選擇階段和訓(xùn)練結(jié)束后結(jié)果評(píng)價(jià)階段對(duì)神經(jīng)網(wǎng)絡(luò)的泛化能力進(jìn)行評(píng)估(Goodfellow et al.,2016).與傳統(tǒng)最優(yōu)化問(wèn)題不同，為了提高訓(xùn)練得到神經(jīng)網(wǎng)絡(luò)的泛化能力，避免“過(guò)擬合”的發(fā)生，神經(jīng)網(wǎng)絡(luò)最優(yōu)化的目標(biāo)函數(shù)并非參與訓(xùn)練的訓(xùn)練集損失函數(shù)，而是未參與訓(xùn)練的驗(yàn)證集與測(cè)試集損失函數(shù)，訓(xùn)練集的損失函數(shù)僅僅為最優(yōu)化過(guò)程提供了方向(Goodfellow et al.,2016).在訓(xùn)練過(guò)程中，利用驗(yàn)證集損失函數(shù)的大小來(lái)對(duì)訓(xùn)練過(guò)程中的各種超參數(shù)進(jìn)行抉擇.以訓(xùn)練步數(shù)的選擇為例，圖2中展示了1000個(gè)步長(zhǎng)學(xué)習(xí)過(guò)程中訓(xùn)練集與驗(yàn)證集損失函數(shù)隨訓(xùn)練步數(shù)的變化曲線.在上文給出的各種超參數(shù)的條件下，經(jīng)過(guò)1000個(gè)步長(zhǎng)時(shí)訓(xùn)練集和驗(yàn)證集的損失函數(shù)均趨近于收斂，且無(wú)明顯過(guò)擬合現(xiàn)象.而經(jīng)過(guò)1000個(gè)步長(zhǎng)后，驗(yàn)證集的損失函數(shù)不再繼續(xù)減小，因此我們將訓(xùn)練步長(zhǎng)選為1000是合理的.在訓(xùn)練結(jié)束階段，用測(cè)試集損失函數(shù)的大小對(duì)訓(xùn)練好的神經(jīng)網(wǎng)絡(luò)的泛化能力進(jìn)行評(píng)估，最終得到的神經(jīng)網(wǎng)絡(luò)訓(xùn)練集的平均損失函數(shù)值為-1.88，驗(yàn)證集的平均損失函數(shù)值為-1.84，測(cè)試集的平均損失函數(shù)值為-1.82，三個(gè)值相差很小，說(shuō)明未參與訓(xùn)練數(shù)據(jù)的損失函數(shù)與參與訓(xùn)練數(shù)據(jù)的損失函數(shù)相當(dāng)，足以證明訓(xùn)練得到的神經(jīng)網(wǎng)絡(luò)模型具有足夠強(qiáng)的泛化能力.

圖2 學(xué)習(xí)曲線(訓(xùn)練集與驗(yàn)證集損失函數(shù)隨訓(xùn)練步數(shù)的變化曲線)Fig.2 Learning curves (loss function′s variation curves of training set and verification set with the number of training steps)

2 計(jì)算結(jié)果與分析

為了驗(yàn)證混合密度神經(jīng)網(wǎng)絡(luò)(MDN)方法的有效性，需要將MDN方法的結(jié)果其與其他獨(dú)立方法的結(jié)果進(jìn)行比較，本研究將其與隨機(jī)采樣方法的結(jié)果進(jìn)行比較.并且根據(jù)MDN方法得到模型參數(shù)后驗(yàn)概率密度分布對(duì)月球內(nèi)部結(jié)構(gòu)模型進(jìn)行采樣，并正演計(jì)算月球內(nèi)部結(jié)構(gòu)模型參數(shù)樣本對(duì)應(yīng)的模擬觀測(cè)值，從而分析了MDN方法的結(jié)果與實(shí)際觀測(cè)值誤差的符合程度．然后基于模型參數(shù)的后驗(yàn)概率密度分布，計(jì)算平均月球內(nèi)部結(jié)構(gòu)模型(Mean模型)、最大后驗(yàn)概率密度對(duì)應(yīng)的模型(MAP模型)以及基于1-σ準(zhǔn)則給出的月球內(nèi)部結(jié)構(gòu)模型(1-σ模型).最后依據(jù)得到的不同觀測(cè)值與模型參數(shù)間的函數(shù)關(guān)系，并利用模型參數(shù)在不同觀測(cè)值誤差范圍內(nèi)的變化幅度為度量，研究了不同觀測(cè)值的觀測(cè)精度對(duì)月球內(nèi)部結(jié)構(gòu)模型的影響.

2.1 MDN方法計(jì)算結(jié)果有效性的獨(dú)立驗(yàn)證

由于似然函數(shù)不是一個(gè)絕對(duì)的評(píng)價(jià)指標(biāo)，其對(duì)應(yīng)的損失函數(shù)值為-1.82，并不能說(shuō)明MDN方法的反演結(jié)果能有效地反映觀測(cè)數(shù)據(jù)中所包含月球內(nèi)部結(jié)構(gòu)的全部信息.因此為了驗(yàn)證MDN方法結(jié)果的有效性，需要將MDN方法計(jì)算的結(jié)果與其他獨(dú)立方法的結(jié)果進(jìn)行比較.本研究用來(lái)對(duì)比的獨(dú)立方法是隨機(jī)采樣方法，其主要步驟如下：(1)利用隨機(jī)采樣方法獲得100000個(gè)月球結(jié)構(gòu)模型樣本，它們參數(shù)的變化范圍如表2所示；(2)利用1.1節(jié)所述正演方法計(jì)算每個(gè)樣本對(duì)應(yīng)的觀測(cè)值；(3)假設(shè)這些觀測(cè)值在各自誤差范圍內(nèi)滿足均勻分布，對(duì)于步驟(2)中觀測(cè)值落在誤差范圍內(nèi)的樣本加以保留，對(duì)于落入誤差范圍外的樣本則拒絕，從而獲得符合觀測(cè)誤差的月球結(jié)構(gòu)模型樣本.所采用的觀測(cè)值及其誤差范圍如表1所示.

圖3給出了MDN方法計(jì)算結(jié)果與隨機(jī)采樣方法獲得的月球內(nèi)部結(jié)構(gòu)模型樣本的比較.圖3中藍(lán)色為MDN方法的結(jié)果，其右側(cè)和上側(cè)是單個(gè)月球內(nèi)部結(jié)構(gòu)模型參數(shù)的邊緣概率密度分布.由于隨機(jī)采樣方法受限于“維度詛咒”，最終100000個(gè)樣本中符合觀測(cè)誤差的樣本只有36個(gè)，在圖中用紅色“×”表示，并且這36個(gè)樣本大部分都落在MDN方法聯(lián)合概率密度分布的高密度區(qū)域，因此MDN方法的計(jì)算結(jié)果能夠有效地表達(dá)隨機(jī)采樣方法結(jié)果的主要特征.以內(nèi)核半徑與外核半徑聯(lián)合概率密度分布為例，MDN方法得到二者的皮爾遜相關(guān)系數(shù)(Pearsonr)為-0.95，p值(顯著性水平)為0，隨機(jī)采樣方法得到二者皮爾遜相關(guān)系數(shù)為-0.91，p值為5.4×10-24；由于二者的皮爾遜相關(guān)系數(shù)均接近-1，說(shuō)明二者存在線性負(fù)相關(guān)關(guān)系，而p值均遠(yuǎn)小于0.05說(shuō)明這種線性相關(guān)是顯著的.MDN方法與隨機(jī)采樣方法的比較結(jié)果可證明MDN方法能有效解決月球內(nèi)結(jié)構(gòu)模型的反演問(wèn)題，得到的解是非平凡解.

圖3 MDN方法結(jié)果與隨機(jī)采樣方法結(jié)果的比較藍(lán)色為MDN方法的結(jié)果，紅色叉號(hào)為隨機(jī)采樣方法的結(jié)果．(a)內(nèi)核半徑與外核半徑聯(lián)合概率分布；(b)外核半徑與低速區(qū)半徑聯(lián)合概率分布；(c)低速區(qū)P波與S波速度聯(lián)合概率密度分布．Fig.3 Comparison of results obtained by the MDN method and the random sampling methodThe blue areas are results of the MDN method,and the red crosses are results of the random sampling method.(a)Joint probability distribution of the inner core radius and the outer core radius.(b)Joint probability distribution of radii of the outer core and the LVZ (Low-Velocity Zone).(c)Joint probability distribution of P-wave and S-wave velocities of LVZ.

為了評(píng)價(jià)MDN方法得到月球內(nèi)部結(jié)構(gòu)模型參數(shù)的精度，需要研究這些月球內(nèi)部結(jié)構(gòu)模型對(duì)應(yīng)的模擬觀測(cè)值與實(shí)際觀測(cè)值的符合程度.對(duì)MDN方法得到月球內(nèi)部結(jié)構(gòu)模型后驗(yàn)概率分布進(jìn)行1000次獨(dú)立采樣，同樣利用正演方法計(jì)算它們對(duì)應(yīng)的模擬觀測(cè)值.圖4給出了這些模擬觀測(cè)值與實(shí)際觀測(cè)值的比較，其中橙色豎線是實(shí)際觀測(cè)值，橙色橫線是其誤差范圍.圖4a中平均密度模擬觀測(cè)值的平均值為3345.64 kg·m-3,標(biāo)準(zhǔn)差為0.31 kg·m-3，與實(shí)際觀測(cè)值3345.56±0.4 kg·m-3相比差別非常小，模擬觀測(cè)值的平均值偏差與標(biāo)準(zhǔn)差的和(0.39 kg·m-3)在實(shí)際觀測(cè)值的誤差范圍以內(nèi)，說(shuō)明反演結(jié)果對(duì)應(yīng)的模擬觀測(cè)值相對(duì)于與實(shí)際觀測(cè)值的偏離程度以及反演結(jié)果的離散程度在誤差允許范圍內(nèi).圖4b中平均轉(zhuǎn)動(dòng)慣量MOI的平均模擬觀測(cè)值為0.393100，標(biāo)準(zhǔn)差為0.000029，其實(shí)際觀測(cè)值為0.393112±0.000012，因此反演結(jié)果對(duì)應(yīng)的平均值相對(duì)于實(shí)際觀測(cè)值的偏差(0.000012)在誤差范圍內(nèi)，然而反演結(jié)果對(duì)應(yīng)的標(biāo)準(zhǔn)差要大于1倍中誤差.圖4c中勒夫數(shù)k2模擬觀測(cè)值的平均值為0.0228，標(biāo)準(zhǔn)差為0.0008，其實(shí)際觀測(cè)值為0.0220±0.0017；與實(shí)際觀測(cè)值相比，模擬觀測(cè)值平均值有0.0008的偏離，但是這個(gè)偏離加上標(biāo)準(zhǔn)差(0.0016)小于實(shí)際觀測(cè)值的一倍中誤差.造成三個(gè)模擬觀測(cè)值精度與實(shí)際精度不一致的主要原因是MDN方法的目標(biāo)函數(shù)是月球內(nèi)部結(jié)構(gòu)模型參數(shù)的似然函數(shù)，如公式(5)所示，并未利用任何觀測(cè)值的精度信息.

圖4 MDN方法獲得的月球內(nèi)部結(jié)構(gòu)模型對(duì)應(yīng)模擬觀測(cè)值的統(tǒng)計(jì)直方圖橙色實(shí)線表示實(shí)際觀測(cè)值以及其誤差范圍．(a)平均密度;(b)平均轉(zhuǎn)動(dòng)慣量;(c)勒夫數(shù)k2.Fig.4 The statistical histogram of the simulated observation values corresponding to the lunar model obtained by the MDN methodThe orange solid lines are the real observations and their error range.(a)Mean density;(b)Mean solid MOI;(c)Love number k2.

MDN方法的反演結(jié)果是月球內(nèi)部結(jié)構(gòu)模型的后驗(yàn)概率密度分布，對(duì)于所有滿足該分布的月球內(nèi)部結(jié)構(gòu)模型的模擬觀測(cè)值而言，它們與觀測(cè)值中的平均密度以及勒夫數(shù)k2的符合程度較好；而對(duì)于觀測(cè)值中的平均轉(zhuǎn)動(dòng)慣量而言，它們的平均值的符合程度較好，而標(biāo)準(zhǔn)差則稍大于實(shí)際誤差范圍.但是本文給出的平均月球內(nèi)部結(jié)構(gòu)模型(Mean模型)和最大后驗(yàn)概率密度模型(MAP模型)對(duì)三個(gè)實(shí)際觀測(cè)值觀測(cè)誤差都符合，進(jìn)一步說(shuō)明本文給出的最優(yōu)化月球內(nèi)部結(jié)構(gòu)模型滿足了所有觀測(cè)值的觀測(cè)誤差.

2.2 月球內(nèi)部結(jié)構(gòu)模型的反演結(jié)果

表2給出了月球內(nèi)部結(jié)構(gòu)模型參數(shù)先驗(yàn)變化區(qū)間以及反演結(jié)果，其中Mean模型是服從后驗(yàn)概率密度分布月球內(nèi)部結(jié)構(gòu)模型的期望，它是所有反演得到月球內(nèi)部結(jié)構(gòu)模型的平均值；MAP模型是最大后驗(yàn)概率密度對(duì)應(yīng)的月球結(jié)構(gòu)模型，它是所有反演得到月球內(nèi)部結(jié)構(gòu)模型當(dāng)中可能性最大的模型.兩個(gè)模型相差很小，原因是反演得到的后驗(yàn)概率密度分布是一個(gè)簡(jiǎn)單的單峰分布，因此可以將MAP模型當(dāng)作反演得到的最優(yōu)化月球結(jié)構(gòu)模型.“1-σ”一列是在1-σ準(zhǔn)則下月球內(nèi)部結(jié)構(gòu)模型參數(shù)的變化范圍，它對(duì)應(yīng)著在Mean模型參數(shù)附近概率密度積分為68.26%的區(qū)域內(nèi)所有月球內(nèi)部結(jié)構(gòu)模型.

表2 月球內(nèi)部結(jié)構(gòu)模型參數(shù)先驗(yàn)變化區(qū)間和反演結(jié)果Table 2 The priori variation intervals and inversion results of the lunar internal structure model parameters

圖5是月球內(nèi)部結(jié)構(gòu)模型中各層半徑、P波和S波速度的垂直剖面圖，其中灰色線條表示W(wǎng)eber等(2011)的月球內(nèi)部結(jié)構(gòu)模型，藍(lán)色虛線表示表2中的MAP模型.在圖5a中的紅色誤差棒表示表2中1-σ模型中的低速區(qū)、外核以及內(nèi)核半徑的變化范圍，它們值分別為212～251 km、305～336 km、334～467 km.Matsuyama等(2016)給出1-σ準(zhǔn)則下的半徑變化范圍分別為6～275 km，194～362 km，345～505 km.本文結(jié)果與他們的結(jié)果相比，給出的變化區(qū)間更窄，這是由于本文假設(shè)月球結(jié)構(gòu)模型中的密度以及低速區(qū)LVZ以外圈層的其他模型參數(shù)均為定值，而在Matsuyama等(2016)的工作中這些參數(shù)都是變量，因此模型的不確定性更大.圖5(b，c)表示P波速度VP以及S波速度VS的垂直剖面，其中紅色虛線的區(qū)間表示表2中1-σ準(zhǔn)則下P波速度與S波速度的變化范圍.Matsumoto等(2015)給出低速區(qū)內(nèi)的P波及S波速度分別為7.1±1.0 km·s-1、2.9±0.5 km·s-1；本文1-σ準(zhǔn)則下給出的6.49～7.70 km·s-1、2.61～3.21 km·s-1與Matsumoto等(2015)的結(jié)果相比區(qū)間更窄，原因與上面半徑類似.MDN方法的反演結(jié)果與前人的反演結(jié)果基本吻合，但是由于簡(jiǎn)化了反演參數(shù)個(gè)數(shù)，因此本文反演的模型不確定度小于前人的結(jié)果.此外Matsuyama等(2016)的結(jié)果顯示低速區(qū)LVZ存在的可能性約為50%，本文結(jié)果給出的低速區(qū)S波速度為2.63 ～ 3.21 km·s-1，明顯低于Weber等(2011)的月幔S波速度4.50 km·s-1，因此本文結(jié)果支持月幔底部存在低速區(qū)的觀點(diǎn).

圖5 月球內(nèi)部結(jié)構(gòu)模型各層半徑(a)、P波(b)及S波(c)速度的垂直剖面灰色實(shí)線為Weber等(2011)的月球內(nèi)部結(jié)構(gòu)模型，左圖中紅色的誤差棒和右邊兩圖的紅色虛線為本研究中1-σ準(zhǔn)則下月球內(nèi)部結(jié)構(gòu)模型變化區(qū)間，藍(lán)色虛線為MAP模型.Fig.5 Vertical profiles of the radii (a),P-wave (b)and S-wave (c)velocities of lunar modeled layersThe gray solid lines are the lunar model given by Weber et al.(2011).The red error bars in the left figure and the red dotted lines in the two right figures are the ranges of the lunar model under the 1-σ criterion in this study.

利用模型參數(shù)間的聯(lián)合概率密度分布可以對(duì)模型參數(shù)之間的相關(guān)性進(jìn)行研究.圖6給出了利用MDN方法得到的內(nèi)外核半徑以及外核半徑與低速區(qū)半徑的聯(lián)合概率密度分布，其中黃色菱形表示表2中MAP模型,紅色五角星則對(duì)應(yīng)于Williams等(2014)的月球內(nèi)部結(jié)構(gòu)模型GPM 1—5，其數(shù)值見(jiàn)表2.如2.1節(jié)所述，內(nèi)外核半徑之間存在明顯的線性負(fù)相關(guān)關(guān)系，這是平均密度、平均轉(zhuǎn)動(dòng)慣量以及勒夫數(shù)k2的觀測(cè)值與誤差范圍共同作用的結(jié)果.從2.3節(jié)的研究結(jié)果來(lái)看，在三個(gè)觀測(cè)值現(xiàn)有的觀測(cè)誤差范圍內(nèi)，k2允許模型參數(shù)在較大的范圍內(nèi)變化，因此它對(duì)模型參數(shù)的約束較弱，因此月球內(nèi)部結(jié)構(gòu)模型的變化范圍主要受到平均密度和平均轉(zhuǎn)動(dòng)慣量約束.在本研究問(wèn)題的假設(shè)背景下，低速區(qū)密度(3.4 kg·m-3)與月幔相同，因此低速區(qū)半徑的改變對(duì)平均密度和平均轉(zhuǎn)動(dòng)慣量沒(méi)有任何影響，而月核中內(nèi)核密度(8.0 kg·m-3)大于外核密度(5.1 kg·m-3).因此為了保持平均密度不變，二者的半徑必然存在線性負(fù)相關(guān)關(guān)系；此外對(duì)于平均轉(zhuǎn)動(dòng)慣量，由于采用的是歸一化的數(shù)值，因此其變化趨勢(shì)與平均密度一致.內(nèi)外核半徑聯(lián)合概率密度分布的結(jié)果與GPM 2—5模型符合程度較好，并且Mean模型和MAP模型幾乎與GPM 3模型重合.而月球內(nèi)部結(jié)構(gòu)模型中的低速區(qū)半徑低于GPM 1—5模型，造成這種差異的主要原因是所采用的固體潮勒夫數(shù)k2為Matsuyama等(2016)給出的0.0220±0.0017，而非Williams等(2014)采用的0.02422±0.00022,該結(jié)果進(jìn)一步說(shuō)明了固體潮勒夫數(shù)k2的觀測(cè)值及其精度對(duì)月球內(nèi)部結(jié)構(gòu)模型影響較為顯著.

圖6 月球內(nèi)部結(jié)構(gòu)模型參數(shù)聯(lián)合概率密度分布(a)內(nèi)核半徑與外核半徑；(b)外核半徑與低速區(qū)半徑.黃色菱形為MAP模型,紅色五角星分別對(duì)應(yīng)了月球內(nèi)部結(jié)構(gòu)模型GPM 1—5．Fig.6 Joint probability density distribution (JPDD)of lunar model parameters(a)The JPDD of the inner core radius and the outer core radius；(b)The JPDD of the outer core radius and the LVZ radius.The yellow diamond refers to the MAP model and the red stars are lunar models GPM 1—5,respectively.

2.3 觀測(cè)值精度對(duì)反演的月球內(nèi)部結(jié)構(gòu)模型的影響

與其他反演方法相比，MDN方法除了可以得到月球內(nèi)部結(jié)構(gòu)模型的后驗(yàn)概率密度分布之外，還可以得到觀測(cè)值與月球內(nèi)部結(jié)構(gòu)模型參數(shù)的函數(shù)關(guān)系，這種函數(shù)關(guān)系為我們?cè)u(píng)估觀測(cè)值精度對(duì)所反演月球內(nèi)部結(jié)構(gòu)模型的影響提供了可能.

圖7給出了在平均密度與平均轉(zhuǎn)動(dòng)慣量為定值(取值如表1所示)條件下，勒夫數(shù)k2與月球內(nèi)部結(jié)構(gòu)模型參數(shù)之間的關(guān)系.橙色實(shí)線為MDN方法得到的Mean模型與勒夫數(shù)k2的函數(shù)關(guān)系，藍(lán)色實(shí)線是利用正演數(shù)值計(jì)算方法得到二者函數(shù)關(guān)系.MDN方法與正演數(shù)值計(jì)算方法結(jié)果的差異反映了MDN方法作為一種函數(shù)近似方法的近似誤差.在一倍觀測(cè)誤差(0.0017)以內(nèi)，函數(shù)的梯度有著明顯的變化，模型參數(shù)與勒夫數(shù)k2之間函數(shù)關(guān)系呈現(xiàn)出明顯的非線性．在Williams等(2014)給出的固體潮勒夫數(shù)k2=0.02422±0.00022附近，內(nèi)核半徑與外核半徑的斜率趨近于0，而低速區(qū)半徑的斜率則是一個(gè)較大的正數(shù)，因此勒夫數(shù)k2在誤差范圍內(nèi)的變化對(duì)低速區(qū)半徑的精度影響較大，而對(duì)內(nèi)外核半徑的精度影響較小.而對(duì)于Matsuyama等(2016)給出的固體潮勒夫數(shù)0.0220±0.0017而言，勒夫數(shù)k2的觀測(cè)誤差對(duì)上述三個(gè)模型參數(shù)精度的影響程度相當(dāng).此外由于內(nèi)外核半徑以及低速區(qū)半徑存在ri

圖7 月球內(nèi)部結(jié)構(gòu)模型參數(shù)與勒夫數(shù)k2的函數(shù)關(guān)系(a)內(nèi)核半徑；(b)外核半徑；(c)低速區(qū)半徑.Fig.7 Function relationship between lunar model parameters and the Love number k2(a)Radius of the inner core；(b)Radius of the outer core；(c)Radius of the LVZ.

圖8給出了觀測(cè)誤差為0時(shí)，利用MDN方法以及表1中的觀測(cè)數(shù)據(jù)得到月球內(nèi)部結(jié)構(gòu)模型的后驗(yàn)概率密度分布，其中紅色五角星為GPM1～5模型.圖8中內(nèi)外核半徑以及低速區(qū)半徑均在275 km附近，而2.2節(jié)中的平均月球內(nèi)部結(jié)構(gòu)模型Mean內(nèi)核半徑為231 km，外核半徑為321 km，低速區(qū)半徑為414 km.造成二者差異的原因主要是因?yàn)橛^測(cè)值k2=0.0220趨近于其取值范圍的下邊界，因此對(duì)于Matsuyama等(2016)給出的固體潮勒夫數(shù)k2=0.0220±0.0017,其有效變化范圍為k2=0.0220～0.0237，如果k2為0.0220時(shí)反演結(jié)果表明內(nèi)核半徑為275 km而外核與低速區(qū)厚度為0，而如果k2為0.0237時(shí)反演得到內(nèi)核半徑為218 km，外核半徑為318 km,低速區(qū)半徑為480 km；在假設(shè)k2服從均勻分布的條件，其平均值為0.02285，則反演結(jié)果便如圖6中所給出.這也是圖4反演結(jié)果中勒夫數(shù)k2模擬觀測(cè)值平均值相對(duì)于實(shí)際觀測(cè)值有明顯偏離的主要原因.

表3 觀測(cè)值在其誤差范圍內(nèi)變化時(shí)導(dǎo)致的月球內(nèi)部結(jié)構(gòu)模型參數(shù)變化幅度
Table 3 Variation range of lunar model parameters caused by variations of observations within their errors

觀測(cè)值精度對(duì)月球內(nèi)部結(jié)構(gòu)模型的影響主要由兩個(gè)方面決定：第一個(gè)方面是問(wèn)題的物理背景，月球內(nèi)部結(jié)構(gòu)模型參數(shù)通過(guò)正演過(guò)程得到對(duì)應(yīng)的模擬觀測(cè)值，因此模擬觀測(cè)值與模型參數(shù)之間通過(guò)控制方程建立起如圖7所示的函數(shù)關(guān)系.在不同觀測(cè)值附近函數(shù)斜率不同，因此月球內(nèi)部結(jié)構(gòu)模型對(duì)模擬觀測(cè)值的敏感程度也不同.第二個(gè)方面是實(shí)際觀測(cè)值的概率密度分布.如果實(shí)際觀測(cè)值觀測(cè)誤差為0，那么反演得到的模型參數(shù)會(huì)是一個(gè)準(zhǔn)確值；而當(dāng)用概率模型描述觀測(cè)值的精度時(shí)，實(shí)際觀測(cè)值并不是常數(shù)而是一個(gè)隨機(jī)變量.實(shí)際觀測(cè)值作為隨機(jī)變量在誤差范圍內(nèi)變化時(shí)，對(duì)應(yīng)的模型參數(shù)會(huì)在函數(shù)的某一段區(qū)間內(nèi)上下浮動(dòng)，其變化幅度的大小取決于誤差范圍的大小以及函數(shù)在各點(diǎn)上的斜率.

從圖7可以看出勒夫數(shù)k2與模型參數(shù)的函數(shù)在其誤差范圍內(nèi)具有明顯的非線性，因此利用函數(shù)在觀測(cè)值處的斜率分析模型參數(shù)可能的變化范圍是不合理的.此外考慮到其他模擬觀測(cè)值與模型參數(shù)的函數(shù)可能是非單調(diào)性的，因此用在觀測(cè)值誤差范圍內(nèi)模型參數(shù)最大值與最小值的差來(lái)度量模型參數(shù)的變化范圍是一種更為合理的方法.為了比較不同觀測(cè)值對(duì)模型精度的影響，在其他兩個(gè)觀測(cè)值為確定值的條件下，計(jì)算單個(gè)觀測(cè)值在其誤差范圍內(nèi)變化時(shí)對(duì)應(yīng)月球內(nèi)部結(jié)構(gòu)模型參數(shù)的變化幅度.表3給出了三個(gè)觀測(cè)值在各自一倍中誤差范圍內(nèi)對(duì)應(yīng)月球內(nèi)部結(jié)構(gòu)模型參數(shù)的變化幅度.表3結(jié)果顯示在觀測(cè)值的誤差范圍內(nèi)，勒夫數(shù)k2對(duì)應(yīng)的月球內(nèi)部結(jié)構(gòu)模型參數(shù)變化幅度遠(yuǎn)遠(yuǎn)大于其他兩個(gè)觀測(cè)值，例如對(duì)于低速區(qū)S波速度而言，勒夫數(shù)k2的影響幅度為202.79 m·s-1，約是MOI影響幅度1.91 m·s-1的200倍；對(duì)于其他月球內(nèi)部結(jié)構(gòu)模型參數(shù)而言，k2的影響同樣遠(yuǎn)大于MOI和平均密度的影響.

內(nèi)核外核低速區(qū)半徑ri(km)P波速度VPi(m·s-1)S波速度VSi(m·s-1)半徑ro(km)P波速度VPo(m·s-1)半徑rl(km)P波速度VPl(m·s-1)S波速度VSl(m·s-1)勒夫數(shù)k260.5952.9735.3648.4319.46198.50105.42202.79平均轉(zhuǎn)動(dòng)慣量MOI0.357.681.430.281.030.354.691.91平均密度ρ(kg·m-3)2.5927.66.072.183.62.3821.03.21

3 結(jié)論

利用混合密度神經(jīng)網(wǎng)絡(luò)(MDN)方法以Williams等(2014)基于GRAIL任務(wù)給出的月球平均密度和平均轉(zhuǎn)動(dòng)慣量、Matsuyama等(2016)給出的月球彈性勒夫數(shù)k2作為觀測(cè)約束，對(duì)月球內(nèi)外核以及低速區(qū)的半徑、P波和S波速度共計(jì)8個(gè)參數(shù)進(jìn)行了反演，得到了月球內(nèi)部結(jié)構(gòu)模型的后驗(yàn)概率密度分布，獲得了最優(yōu)化月球內(nèi)部結(jié)構(gòu)模型，即最大后驗(yàn)概率密度對(duì)應(yīng)的MAP模型.除此之外還計(jì)算了后驗(yàn)概率密度分布的平均月球內(nèi)部結(jié)構(gòu)模型(Mean模型)以及基于1-σ準(zhǔn)則給出的月球內(nèi)部結(jié)構(gòu)模型范圍.Mean模型與MAP模型差別較小且后驗(yàn)概率密度分布是單峰分布說(shuō)明以MAP模型作為最優(yōu)化月球內(nèi)部結(jié)構(gòu)模型是合理的.反演得到的月球內(nèi)部結(jié)構(gòu)模型中，內(nèi)核半徑與外核半徑呈現(xiàn)出明顯的線性負(fù)相關(guān)關(guān)系，這是因?yàn)樵谌齻€(gè)觀測(cè)值觀測(cè)精度內(nèi)勒夫數(shù)k2允許模型參數(shù)在較大的區(qū)間內(nèi)變化，而與之相比平均密度和平均轉(zhuǎn)動(dòng)慣量在誤差允許范圍內(nèi)幾乎可以認(rèn)為是不變的.反演結(jié)果中低速區(qū)S波速度(2.63～3.2 km·s-1)明顯低于月幔S波速度(4.50 km·s-1)，因此本文結(jié)果支持月幔底部存在低速區(qū).

相比于其他反演方法，MDN方法的優(yōu)勢(shì)主要體現(xiàn)在三個(gè)方面：(1)與傳統(tǒng)的最小二乘法反演方法相比，MDN方法可以解決反演的不唯一性與非線性性；(2)與反演中經(jīng)常采用的蒙特卡洛方法相比，MDN方法避免了“維度詛咒”，因此計(jì)算效率更高；(3)MDN方法不僅可以獲得一個(gè)實(shí)際觀測(cè)值對(duì)應(yīng)的月球內(nèi)部結(jié)構(gòu)模型參數(shù)，還可以獲得實(shí)際觀測(cè)值與月球內(nèi)部結(jié)構(gòu)模型參數(shù)之間的函數(shù)關(guān)系，為進(jìn)一步分析不同觀測(cè)值觀測(cè)誤差對(duì)月球內(nèi)部結(jié)構(gòu)模型的影響提供了可能.

通過(guò)研究勒夫數(shù)k2與模型參數(shù)間的函數(shù)關(guān)系，數(shù)值結(jié)果表明勒夫數(shù)k2存在一個(gè)約為0.022的下限值.利用在觀測(cè)誤差范圍內(nèi)模型參數(shù)的變化分布作為度量，研究了不同觀測(cè)值精度對(duì)模型的影響，結(jié)果顯示現(xiàn)有觀測(cè)精度條件下勒夫數(shù)k2的精度對(duì)月球內(nèi)部結(jié)構(gòu)模型的影響明顯要大于其他兩個(gè)觀測(cè)值.在將來(lái)的月球觀測(cè)任務(wù)中，期望能首先考慮提高固體潮勒夫數(shù)的觀測(cè)精度，這樣就能對(duì)月球內(nèi)部結(jié)構(gòu)進(jìn)行更加準(zhǔn)確的約束.