• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    “嫦娥5號”登陸候選地Mons Rümker的光照與溫度特征分析*

    2020-06-30 12:14:10鐘振張騰張杰陳世國
    物理學(xué)報 2020年11期
    關(guān)鍵詞:風(fēng)化層熱流月球

    鐘振 張騰 張杰 陳世國

    1) (貴州師范大學(xué)物理與電子科學(xué)學(xué)院, 貴陽 550001)2) (貴州師范學(xué)院地理與資源學(xué)院, 貴陽 550018)(2020 年1 月16日收到; 2020 年3 月15日收到修改稿)

    即將開展的“嫦娥5號”探月任務(wù), 將使我國在月球上首次實現(xiàn)無人鉆井取樣. 考慮到實際探測活動與探測區(qū)域的光照和溫度有關(guān), 有必要對研究區(qū)域的光照和溫度特征進(jìn)行分析, 為此, 本文利用SPICE系統(tǒng)對“嫦娥5號”候選登陸區(qū)Mons Rümker高原的實時光照進(jìn)行計算. 發(fā)現(xiàn)求解的相對光強(qiáng)度分布與日本SELENE衛(wèi)星提供的早晨光照影像一致, 驗證了光照算法及計算程序的合理性. 以此為基礎(chǔ), 利用1維熱傳導(dǎo)模型, 對候選登陸區(qū)風(fēng)化層不同深度的溫度進(jìn)行仿真分析. 結(jié)果表明風(fēng)化層溫度在近表面區(qū)域受光照的影響較大, 隨著深度的增加, 光照影響逐漸減弱. 到達(dá)0.57 m深度時, 風(fēng)化層溫度不再變化. 為確保鉆井任務(wù)的開展, 實際鉆井作業(yè)應(yīng)考慮風(fēng)化層內(nèi)外溫度差異引起的應(yīng)力不均. 考慮“嫦娥5號”的實際鉆井深度遠(yuǎn)大于0.57 m, 應(yīng)能測量到常溫層的熱流值, 后續(xù)探月任務(wù)可考慮搭載熱流探測設(shè)備, 以促進(jìn)月球科學(xué)研究的發(fā)展.

    1 引 言

    根據(jù)我國探月工程計劃, 繼“嫦娥”1號、2號、3號和4號的成功實施, 我國即將開展探月工程第三期“回”[1-3]. 該階段計劃發(fā)射“嫦娥5號”月球車,在月表進(jìn)行軟著陸, 對月表風(fēng)化層進(jìn)行鉆井, 提取相關(guān)巖石樣品并返回地球, 屆時, 我國將成為繼蘇聯(lián)和美國之后, 全球第3個實現(xiàn)月球采樣并返回的國家[4]. 有關(guān)“嫦娥5號”登陸區(qū)的選擇, 文獻(xiàn)[5]認(rèn)為風(fēng)暴洋(oceanus procellarum)北面的Rümker區(qū)域, 曾經(jīng)擁有較長的火山活動, 形成了多重地質(zhì)單元, 具有不同的元素組成, 采集該區(qū)域的巖石樣本具有重要的地質(zhì)和科學(xué)意義. 該區(qū)域的Mons Rümker位于風(fēng)暴洋北部, 中心坐標(biāo)為(40°N, 58°W),擁有古火山活動時遺留的月溪、穹隆和熔巖流管道. 其地形起伏在1000—2000 m之間, 文獻(xiàn)[6-8]認(rèn)為它屬于火山中心, 文獻(xiàn)[9]認(rèn)為它屬于大型盾型火山. 文獻(xiàn)[10]對Mons Rümker的地形地貌,以及地質(zhì)特征進(jìn)行分析, 提出將該區(qū)域作為“嫦娥5號”登陸候選區(qū), 采集該區(qū)域的巖石樣品, 有助于理解月球火山作用特征及其演化歷史.

    文獻(xiàn)[11]對Mons Rümker高原的累積光照條件進(jìn)行研究, 發(fā)現(xiàn)3年時間擁有30%至53%的光照率, 且地基觀測系統(tǒng)的通信條件較好. 優(yōu)良的光照條件可以快速地提升月表溫度, 但由于月表風(fēng)化層具有良好的絕熱特性, 使得光照對風(fēng)化層內(nèi)部溫度的貢獻(xiàn)有限[12-15]. 盡管文獻(xiàn)[12-15]對風(fēng)化層溫度分布進(jìn)行了估算, 但太陽光照方向主要依據(jù)理論公式求得, 與現(xiàn)實情況存在一定的偏差. 另外,他們在估算光照引起的溫度分布時, 并沒有考慮地形對太陽光的遮擋, 而地形遮擋效應(yīng)對實際探測活動至關(guān)重要. 為此, 本文利用美國NASA (National Aeronautics and Space Administration)開發(fā)的SPICE(Spacecraft, Planet, Instrument, Cameramatrix, Events)系統(tǒng), 精確確定不同時刻的太陽光照方向. 結(jié)合文獻(xiàn)[16,17]提供的高分辨率月球數(shù)字地形模型DTM (digital terrain model), 估算顧及地形遮擋的實時太陽光照. 同時, 考慮月球內(nèi)部熱流, 對Mons Rümker高原風(fēng)化層的溫度分布進(jìn)行計算, 以期為“嫦娥5號”鉆井采樣活動, 以及后繼探月活動的開展提供一定程度的參考.

    2 地形模型和光照條件

    Mons Rümker高原的平均面積約為4000 km2,地形圖如圖1白色方框區(qū)域所示. 本文采用的數(shù)字地形模型, 來自LRO (lunar reconnaissance orbiter)的激光測高數(shù)據(jù)LOLA (lunar orbiter laser altimeter). LRO自2009年發(fā)射以來, 其激光測距載荷LOLA取得了高精度的月球全球地形數(shù)據(jù),促進(jìn)了月表永久陰影區(qū)及水冰的研究, 同時也確保了月球全球高精度大地參考框架的建立, 為后續(xù)載人登月提供了安全保障[18]. LRO目前仍處于在軌狀態(tài), 盡管后期拓展任務(wù)階段因調(diào)軌操作, 無法采集月球北半球地形數(shù)據(jù), 但激光測距載荷LOLA仍然在不斷地更新數(shù)據(jù), 地形格網(wǎng)數(shù)據(jù)已更新至2019年(https://pds-geosciences.wustl.edu/missions/lro/lola.htm). 本文采用的地形模型, 來自LOLA的格網(wǎng)數(shù)據(jù)LDEM_512_00N_45N_270_360 (1/512° × 1/512°), 其分辨率約為59.2 m ×59.2 m. 其中, 黑色方塊表示文獻(xiàn)[10]建議的“嫦娥5號”候選登陸點(303.34°E, 40.11°N).

    圖 1 Mons Rümker區(qū)域地形圖, 如白色方框所示, 其中黑色方塊表示文獻(xiàn)[10]建議的“嫦娥5號”登陸點(303.34°E,40.11°N)Fig. 1. Topography around Mons Rümker region, which is figured out with a white box. The black box indicates the candidate landing site of CE-5 proposed by reference [10],and this site is centered at (303.34°E, 40.11°N).

    如圖2所示, 為了判斷目標(biāo)點A是否有太陽光照射, 以及判斷A點是否被鄰區(qū)地形點B遮擋,可以通過點A的太陽高度角, 以及比較BC和DC的大小來判斷. 過點A作等效球面(圓弧型虛線所示), 與OB相交于點C, 點A和點B間的中心角為θ. 假定點A的太陽高度角為α, 它表示太陽入射光線與目標(biāo)點水平地面間的夾角. 若α ≤ 0,文獻(xiàn)[11]的研究表明, 對于高海拔區(qū)域, 太陽光有可能照射至目標(biāo)點A. 如圖1所示, 本文研究區(qū)域Mons Rümker高原, 其最大高程不超過—1.5 km,不屬于高海拔區(qū)域, 因此, 本文不考慮文獻(xiàn)[11]的高海拔特殊情況, 認(rèn)為α ≤ 0時, 太陽光不能照射至目標(biāo)點A. 當(dāng)α > 0時, 若DC < BC, 則目標(biāo)點A被鄰區(qū)地形點B遮擋, 太陽光無法到達(dá)點A; 當(dāng)DC≥BC時, 則不受點B的地形遮擋. 考慮月球表面地形高低的最大落差, 文獻(xiàn)[11]表明最大中心角θ = 8°, 通過不斷地減小中心角, 可以依次判斷距離目標(biāo)點A遠(yuǎn)近不同的地形遮擋效應(yīng). BC的大小可以根據(jù)數(shù)字地形模型求得, 依據(jù)文獻(xiàn)[11], 可得

    圖 2 太陽光照條件示意圖Fig. 2. Schematic of illumination condition.

    其中AC可以根據(jù)中心角θ和OA求得; OA可以根據(jù)數(shù)字地形模型DTM求得; α為實時太陽光高度角, 可以根據(jù)SPICE系統(tǒng)由行星歷表求得.

    3 熱傳導(dǎo)方程及邊界條件

    已知光照強(qiáng)度, 根據(jù)熱傳導(dǎo)方程, 可以估算風(fēng)化層的溫度變化. 月表風(fēng)化層的厚度, 就整個研究區(qū)域而言相對較小, 因此, 通常采用一維熱傳導(dǎo)模型來計算溫度分布[12-15]. 就溫度T和風(fēng)化層深度z而言, 一維模型的控制方程為

    其中, ρ和cp分別表示風(fēng)化層的密度和比熱容,k表示熱傳導(dǎo)系數(shù), 這些參數(shù)通常與風(fēng)化層的深度和溫度分布有關(guān). 根據(jù)文獻(xiàn)[13,14,19-21], 有

    式中相關(guān)參數(shù)的取值如表1所列.

    表 1 相關(guān)參數(shù)取值Table 1. Values of parameters used in study.

    為了求解(2)式, 需要指定兩個邊界條件, 一個在月表, 另一個在風(fēng)化層底部. 月表邊界取決于太陽光照及月表紅外輻射, 風(fēng)化層底部取決于月球內(nèi)部熱流分布. 對于深度為zb的底部邊界z = zb,有如下關(guān)系[14]:

    其中, Q表示風(fēng)化層底部的熱流值, 根據(jù)Apollo任務(wù)測量值[22], 取Q = 0.018 W·m—2.

    對于月表邊界溫度Ts, 假定太陽光照加熱率為Qs, 考慮紅外輻射后的月表溫度梯度有[13-15]:

    其中a, b為系數(shù), 參考文獻(xiàn)[13], 取a = 0.06, b = 0.25;A0表示月球表面的邦德反照率(又稱球面反照率),該參數(shù)用于描述太陽于正中天, 太陽光垂直入射時, 天體表面的反射系數(shù), 參考文獻(xiàn)[12], 本文取A0= 0.12. 太陽實時入射流Fsun和入射角β通常與月表計算點和太陽之間的實時距離r有關(guān), 該距離可以通過SPICE系統(tǒng)精確求得. 假定地球與太陽間的平均距離為r0, 對應(yīng)的太陽常數(shù)為S, 則入射流Fsun為

    其中, r0為日-地平均間距. 參考文獻(xiàn)[25], 取r0=1.49598261 × 1011m, S = 1361 W·m—2. 由(8)式和(9)式可知, 求取太陽光照加熱率, 關(guān)鍵在于入射角β的求取. 參考文獻(xiàn)[26], 入射角β與研究點的坡度λ、方位角γ、太陽天頂距θz、太陽方位角γs有關(guān), 關(guān)系為

    其中, 方位角γ可由研究點的坡向求出; 太陽天頂距θz及其方位角γs可由SPICE系統(tǒng)提供的實時太陽星下點坐標(biāo)以及文獻(xiàn)[26]的方法求出; 研究點的坡度λ可由文獻(xiàn)[17,18]提供的月球全球坡度分布圖求出.

    4 結(jié)果與分析

    為了測試算法與程序的合理性, 圖3給出了Mons Rümker高原的早晨光照分布. 該圖對應(yīng)月球地方時tm= 06:30:30, Mons Rümker高原在早晨相應(yīng)時刻, 地表對太陽光反射的相對光強(qiáng)度分布. 其中圖3(a)來自日本SELENE衛(wèi)星提供的早晨光照影像(http://darts.isas.jaxa.jp/planet/pdap/selene/), 由 影 像TCO_MAPm04_N42E 300N39E303SC和TCO_MAPm04_N42E303N39 E306SC拼接得到, 其中白色和黑色分別代表光照強(qiáng)度最強(qiáng)和最弱, 其他顏色表示光照強(qiáng)度在兩者間的變化. 圖3(a)表明太陽光自東向西, 東邊光照強(qiáng)度強(qiáng), 西邊光照弱; 迎光方向出現(xiàn)光亮區(qū)域, 而背光方向則由于地形遮擋出現(xiàn)陰暗區(qū)域. 采用與圖3(a)相同的光照時刻, 圖3(b)表示基于LOLA地形模型、(1)式和(10)式求解的相對光強(qiáng)度RII(relative intensity of illumination). 圖3(b)中等高線表示Mons Rümker高原, 相對半徑為1737.151 km的參考球面的地形高. 由圖3(b)可知: 1)東邊相對光強(qiáng)度強(qiáng), 西邊相對光強(qiáng)度弱;2)東邊迎光方向部分區(qū)域, 由于坡向正對太陽光,相對光強(qiáng)度較強(qiáng), 呈現(xiàn)出白色; 3)西邊背光方向部分區(qū)域, 由于地形遮擋, 相對光強(qiáng)度接近于零, 呈現(xiàn)出黑色. 對比圖3(a)和圖3(b), 可以發(fā)現(xiàn)除少部分小區(qū)域外, 兩者光照強(qiáng)度的分布總體上一致. 部分小區(qū)域出現(xiàn)偏差, 主要是由于光照計算采用的LOLA地形模型的分辨率, 低于日本SELENE衛(wèi)星提供的早晨光照影像所致. 本文計算的光強(qiáng)度分布, 在總體上與觀測結(jié)果一致, 說明整個算法及計算程序具有一定的合理性, 可以進(jìn)一步應(yīng)用于Mons Rümker高原風(fēng)化層溫度的估算.

    圖 3 Mons Rümker區(qū)域早晨光照圖, 對應(yīng)月球地方時tm =06:30:30 (a)日本SELENE衛(wèi)星提供的早晨光照圖; (b)本文計算的與圖3(a)相同時刻的實時光照圖Fig. 3. Morning map of illumination over Mons Rümker at the lunar local time tm = 06:30:30: (a) Japan’s SELENE morning map of illumination; (b) our estimated relative intensity of illumination at the same time of Fig. 3(a).

    利用本文計算的實時光照強(qiáng)度, 根據(jù)(2)—(7)式, 可以求出Mons Rümker高原風(fēng)化層的實時溫度分布. 為了計算風(fēng)化層不同深度處的溫度,利用文獻(xiàn)[14]提供的標(biāo)準(zhǔn)化有限差分方法對(2)式進(jìn)行數(shù)值計算. 圖4給出了Mons Rümker高原不同時刻的表面溫度, 其中圖4(a)對應(yīng)協(xié)調(diào)世界時UTC(universal time coordinated): 2020年10月28日11點30分00秒, 簡寫為UTC 2020-10-28T 11:30:00, 后文所有時間依此格式表述. 此時, 由于沒有太陽光照, 大部分區(qū)域的溫度在80 K左右. 至圖4(b)時(UTC 2020-10-29T 06:45:00),隨著太陽光的到來, 東邊溫度逐漸升高, 特別是坡向正對太陽光的區(qū)域, 溫度一度接近200 K, 西邊由于背光, 溫度維持在120 K左右. 隨著太陽高度角的增大, 類似情況也出現(xiàn)在圖4(c)中, 大部分區(qū)域的溫度升至200 K左右, 高溫區(qū)域一度接近280 K, 該圖對應(yīng)的時刻為UTC 2020-10-30T 06:45:00. 到達(dá)正午時分(如圖4(d)所示, 對應(yīng)時刻為UTC 2020-11-02T 04:45:00), 大部區(qū)域的溫度升至360 K左右, 少部分區(qū)域接近或超過400 K. 至圖4(e)時, 由于太陽高度角的下降, 光照減少, 表面溫度也逐漸下降. 此時, Mons Rümker高原處于月球地方時的下午時分, 太陽光照至西向東, 西邊迎光方向光照強(qiáng), 東邊背光方向光照弱,對應(yīng)時刻為UTC 2020-11-12T 02:45:00. 如圖4(f)所示(對應(yīng)時刻為UTC 2020-11-12T 17:00:00), 隨著太陽高度角的進(jìn)一步降低, 大部分區(qū)域的溫度降至120 K左右, 西邊迎光方向部分區(qū)域的溫度維持在200 K左右. 盡管溫度下降, 部分背光區(qū)域的溫度一度高于200 K, 這些區(qū)域恰好對應(yīng)圖4(b)—(d)中的高溫區(qū)域. 在表面光照強(qiáng)時, 這些高溫區(qū)域的熱量不斷地向風(fēng)化層底部傳遞, 至圖4(f)時, 盡管表面光照減弱, 已傳遞至風(fēng)化層內(nèi)部的熱量反向傳遞至風(fēng)化層表面, 使得表面溫度升高. 由于這些區(qū)域內(nèi)部溫度高于其他區(qū)域, 至圖4(f)時, 盡管光照減弱, 得到內(nèi)部熱量的補(bǔ)充, 表面溫度仍然高于其他區(qū)域.

    圖 4 Mons Rümker區(qū)域表面溫度分布圖 (a)對應(yīng)時刻為UTC 2020-10-28T 11:30:00; (b)對應(yīng)時刻為UTC 2020-10-29T 06:45:00; (c)對應(yīng)時刻為UTC 2020-10-30T 06:45:00; (d)對應(yīng)時刻為UTC 2020-11-02T 04:45:00; (e)對應(yīng)時刻為UTC 2020-11-12T 02:45:00; (f)對應(yīng)時刻為UTC 2020-11-12T 17:00:00Fig. 4. Surface temperature distribution with time over Mons Rümker plateau: (a) Time at UTC 2020-10-28T 11:30:00; (b) time at UTC 2020-10-29T 06:45:00; (c) time at UTC 2020-10-30T 06:45:00; (d) time at UTC 2020-11-02T 04:45:00; (e) time at UTC 2020-11-12T 02:45:00; (f) time at UTC 2020-11-12T 17:00:00.

    這種結(jié)論也可以由圖5得到佐證. 圖5表示Mons Rümker高原風(fēng)化層5 cm深度的溫度分布,對應(yīng)時刻與圖4各分圖一致. 由圖5(a)—5(c)可知, 隨著太陽光照增加, 盡管表面溫度在升高(如圖4(a)—4(c)所示), 但風(fēng)化層5 cm深度處的溫度始終維持在200 K左右, 至正午時分圖5(d)時, 底部溫度才有所升高. 這說明月球風(fēng)化層的導(dǎo)熱能力較弱, 具有一定的絕熱性能, 該結(jié)果與文獻(xiàn)[14,15]的一致. 隨著下午時分太陽光照的減弱, 盡管表面溫度下降至120 K左右(如圖4(e)和圖4(f)所示), 但5 cm深度的溫度始終維持在260 K左右,部分區(qū)域甚至達(dá)到300 K. 此時, 隨著表面溫度下降, 內(nèi)部溫度高于月表溫度, 熱量由內(nèi)向外傳遞,使得圖4(f)部分背光區(qū)域的溫度高于鄰近區(qū)域.

    為了進(jìn)一步探究風(fēng)化層的溫度變化, 圖6給出了溫度隨深度變化的剖面圖, 剖面方向沿圖1藍(lán)色直線的方向, 橫軸表示經(jīng)度, 縱軸表示深度, 對應(yīng)時刻與圖4和圖5一致. 由圖6(a)—(c)可知, 隨著太陽光照的增加, 風(fēng)化層5 cm深度以內(nèi)的溫度不斷地增加, 5 cm以下的溫度變化較弱, 這與圖5(a)—(c)的結(jié)果一致. 至正午時刻圖6(d)時,表面溫度升至最大值, 熱量向風(fēng)化層內(nèi)部傳遞, 內(nèi)部溫度不斷升高. 圖6(d)中黑色曲線表示沿剖面方向的表面地形輪廓, 其幅度表示地形的相對變化, 其大小與圖6縱軸刻度范圍無關(guān), 該圖表明溫度的變化與地形有關(guān). 至圖6(e)和圖6(f)時, 內(nèi)部儲存的熱量向表面?zhèn)鬟f, 使得表面對應(yīng)區(qū)域的溫度高于鄰近區(qū)域.

    為了進(jìn)一步研究0.2 m深度以下的溫度變化,圖7給出了兩個參考點不同深度處溫度隨時間的變化關(guān)系. 其中, 圖7(a)表示“嫦娥5號”候選登陸點(圖1所示黑色方塊)的溫度變化, 圖7(b)表示另一參考點(圖1所示黑色五角星)的溫度變化.由圖7可知, 兩個參考點的溫度在0.27 m深度的溫度最小值在239 K左右. 溫度變化幅度隨著深度的降低逐漸減小, 到達(dá)0.57 m深度時, 溫度幾乎不再變化, 保持在241.5 K左右. 結(jié)合圖4—圖7,保守估計風(fēng)化層的常溫層深度在0.6 m. 參考文獻(xiàn)[4], “嫦娥5號”的鉆井深度接近2 m, 在進(jìn)行鉆井作業(yè)及仿真分析時, 有必要考慮風(fēng)化層內(nèi)外溫度的差異; 另外, 考慮到常溫層深度在0.6 m左右,“嫦娥5號”的鉆井深度達(dá)2 m, 能探測到常溫層的熱流. 結(jié)合“嫦娥5號”的鉆井經(jīng)驗, 后期探月工程可考慮搭載探測月球內(nèi)部熱流值的載荷, 以促進(jìn)月球內(nèi)部結(jié)構(gòu)及熱演化研究的發(fā)展.

    圖 5 Mons Rümker區(qū)域底部5 cm深度的溫度分布圖 (a)對應(yīng)時刻為UTC 2020-10-28T 11:30:00; (b)對應(yīng)時刻為UTC 2020-10-29T 06:45:00; (c)對應(yīng)時刻為UTC 2020-10-30T 06:45:00; (d)對應(yīng)時刻為UTC 2020-11-02T 04:45:00; (e)對應(yīng)時刻為UTC 2020-11-12T 02:45:00; (f)對應(yīng)時刻為UTC 2020-11-12T 17:00:00Fig. 5. Subsurface temperature distribution at the depth of 5 cm over Mons Rümker plateau: (a) Time at UTC 2020-10-28T 11:30:00; (b) time at UTC 2020-10-29T 06:45:00; (c) time at UTC 2020-10-30T 06:45:00; (d) time at UTC 2020-11-02T 04:45:00;(e) time at UTC 2020-11-12T 02:45:00; (f) time at UTC 2020-11-12T 17:00:00.

    圖 6 剖面溫度沿圖1所示經(jīng)度方向的分布 (a)對應(yīng)時刻為UTC 2020-10-28T 11:30:00; (b)對應(yīng)時刻為UTC 2020-10-29T 06:45:00; (c)對應(yīng)時刻為UTC 2020-10-30T 06:45:00; (d)對應(yīng)時刻為UTC 2020-11-02T 04:45:00; (e)對應(yīng)時刻為UTC 2020-11-12T 02:45:00; (f)對應(yīng)時刻為UTC 2020-11-12T 17:00:00. 圖6(d)所示黑色曲線表示圖1藍(lán)線方向的表面地形輪廓Fig. 6. Temperature variation along the longitude direction shown in Fig.1: (a) Time at UTC 2020-10-28T 11:30:00; (b) time at UTC 2020-10-29T 06:45:00; (c) time at UTC 2020-10-30T 06:45:00; (d) time at UTC 2020-11-02T 04:45:00; (e) time at UTC 2020-11-12T 02:45:00; (f) time at UTC 2020-11-12T 17:00:00. The black carves in Fig. 6(d) represents the surface topography along the same blue line direction displayed in Fig. 1.

    5 討 論

    前文在估算Mons Rümker高原風(fēng)化層的溫度分布時, 控制風(fēng)化層密度分布的H參數(shù)取值為0.06 m, 熱流值Q選用的是Apollo 15號和Apollo 17號測量的平均值(Q = 0.018 W·m—2), 不同參數(shù)取值有可能引起溫度波動, 為此, 有必要討論相關(guān)參數(shù)的適用性. 參考文獻(xiàn)[14], 相關(guān)參數(shù)如熱傳導(dǎo)系數(shù)、風(fēng)化層密度和比熱容等參數(shù)主要受H參數(shù)的調(diào)控, 因此, 僅討論H參數(shù)對溫度的影響. 參考文獻(xiàn)[14], H參數(shù)的最小值為0.02 m, 最大值為0.09 m. 參考文獻(xiàn)[22], 熱流測量的最大值在Apollo 15號登陸點, 約為0.021 W·m—2. 文獻(xiàn)[27]利用數(shù)值方法, 估算的全球熱流值約為0.012 W·m—2, 小于Apollo任務(wù)的測量值. 為了直觀方便, 圖8和表2給出了H參數(shù)和熱流值Q取極值時, 參考點(如圖1黑色五角星所示)溫度隨深度的變化. 其中黑色實線和紅色虛線, 表示Q =0.018 W·m—2, H參數(shù)分別取0.02 m和0.09 m時,參考點溫度隨深度的變化; 綠色和藍(lán)色虛線表示H = 0.06 m, 熱流值分別取Q = 0.012 W·m—2和Q = 0.021 W·m—2時的溫度分布. 圖8(a)—8(c)分別對應(yīng)月球地方時tm= 06:30:30, tm= 12:30:30和tm= 18:30:30, 對應(yīng)參考點的早晨、正午和傍晚時刻. 1)由圖8(a)—8(c)可知, H參數(shù)和熱流值Q取值的不同, 對溫度的變化只有細(xì)微的影響, 絕對誤差不超過10 K; 2)如圖8(c)以及表2所示,H參數(shù)對溫度的影響超過熱流值Q, H參數(shù)主要對風(fēng)化層0.2 m以內(nèi)區(qū)域產(chǎn)生較大的影響, 對底層溫度的影響較弱; 3)如表2所列, 熱流值Q取值的不同, 主要對風(fēng)化層基底溫度產(chǎn)生較小的影響, 絕對誤差不超過1 K. 綜上所述, 同時考慮實際鉆進(jìn)深度超過1 m, H參數(shù)和熱流值Q的差異, 對溫度分布的影響可以忽略不計, 因此, 本文計算結(jié)果具有一定的參考價值.

    圖 8 參考點(圖1中黑色五角星)溫度隨深度的變化 (a)對應(yīng)月球地方時tm = 06:30:30; (b)對應(yīng)月球地方時tm = 12:30:30;(c)對應(yīng)月球地方時tm = 18:30:30Fig. 8. Subsurface temperature variations with depth for the point of black star in Fig. 1: (a) Temperature variations at the lunar local time tm = 06:30:30; (b) temperature variations at the lunar local time tm = 12:30:30; (c) temperature variations at the lunar local time tm = 18:30:30.

    表 2 參考點(圖1中黑色五角星)溫度(單位為K)在不同時刻隨深度的變化Table 2. Temperature (in K) variations with depth for the point of black star in Fig. 1 at various lunar local time.

    6 結(jié) 論

    針對“嫦娥5號”候選登陸區(qū)Mons Rümker高原風(fēng)化層的溫度變化, 利用SPICE系統(tǒng)對實時光照進(jìn)行計算, 在此基礎(chǔ)上, 結(jié)合1維熱傳導(dǎo)模型對風(fēng)化層的溫度進(jìn)行仿真分析. 結(jié)果表明: 1)本文計算的實時相對光照強(qiáng)度分布與日本SELENE衛(wèi)星提供的早晨光照影像總體一致, 說明本文的算法及計算程序具有一定的合理性; 2)依據(jù)求解的光照作為邊界條件, 發(fā)現(xiàn)0.5 cm深度以內(nèi)的風(fēng)化層溫度, 受表面光照的影響較大, 隨著深度的增加, 光照影響逐漸減弱; 3)風(fēng)化層底部0.57 m深度以下的溫度不隨時間變化, 保守估計常溫層的深度在0.6 m左右; 4)決定風(fēng)化層物理特征的H參數(shù), 其大小的改變并不會引起溫度的大幅波動, 對常溫層產(chǎn)生的影響較弱; 5)熱流值Q的不同, 僅對基底邊界溫度產(chǎn)生較小的影響, 實際應(yīng)用中可以不考慮熱流值Q選取的差異.

    猜你喜歡
    風(fēng)化層熱流月球
    到月球上“飆車”
    軍事文摘(2022年12期)2022-07-13 03:11:44
    巖質(zhì)邊坡風(fēng)化層的關(guān)鍵滑動面及其穩(wěn)定性研究
    近地表速度模型精度影響分析
    陪我去月球
    飛碟探索(2022年1期)2022-03-12 14:11:20
    月球上的另一個我
    軍事文摘(2021年18期)2021-12-02 01:28:54
    內(nèi)傾斜護(hù)幫結(jié)構(gòu)控釋注水漏斗熱流道注塑模具
    空調(diào)溫控器上蓋熱流道注塑模具設(shè)計
    流域風(fēng)化層雨洪調(diào)蓄生態(tài)機(jī)理
    聚合物微型零件的熱流固耦合變形特性
    中國塑料(2017年2期)2017-05-17 06:13:24
    坡度在巖石風(fēng)化層解譯中的應(yīng)用
    河北遙感(2015年2期)2015-07-18 11:11:14
    久久精品久久精品一区二区三区| 最新的欧美精品一区二区| 久久免费观看电影| 久久国产精品大桥未久av| 国产不卡av网站在线观看| 伊人久久国产一区二区| 丁香六月欧美| 不卡视频在线观看欧美| 久久久国产欧美日韩av| 一级爰片在线观看| 成人午夜精彩视频在线观看| 久久这里只有精品19| av福利片在线| 熟女少妇亚洲综合色aaa.| 午夜免费鲁丝| 七月丁香在线播放| 黄片播放在线免费| 9热在线视频观看99| 日日啪夜夜爽| 免费黄频网站在线观看国产| 中文字幕人妻丝袜一区二区 | 大码成人一级视频| 毛片一级片免费看久久久久| 女性被躁到高潮视频| 一级黄片播放器| 亚洲精品国产av蜜桃| 午夜免费男女啪啪视频观看| 看免费成人av毛片| 午夜福利免费观看在线| 中文字幕人妻丝袜制服| 亚洲第一区二区三区不卡| 免费观看人在逋| 精品一区二区免费观看| 中文字幕精品免费在线观看视频| 午夜影院在线不卡| 中文天堂在线官网| 免费av中文字幕在线| 下体分泌物呈黄色| 18禁国产床啪视频网站| 69精品国产乱码久久久| 看片在线看免费视频| 老司机靠b影院| 国产av一区二区精品久久| 日本a在线网址| 亚洲成国产人片在线观看| 香蕉久久夜色| 亚洲美女黄片视频| cao死你这个sao货| www.精华液| 色老头精品视频在线观看| 国产成人系列免费观看| 久久久久久久久久久久大奶| 成人18禁高潮啪啪吃奶动态图| 少妇 在线观看| 色综合亚洲欧美另类图片| 国产成人精品无人区| 91麻豆av在线| 成年版毛片免费区| 久久中文看片网| a在线观看视频网站| 天堂影院成人在线观看| 国产成人欧美| 长腿黑丝高跟| 中文字幕人妻熟女乱码| 欧美乱码精品一区二区三区| avwww免费| 欧美在线一区亚洲| 久久午夜综合久久蜜桃| 老司机靠b影院| 黄色视频不卡| 午夜精品久久久久久毛片777| 国产人伦9x9x在线观看| 国产99白浆流出| 91精品三级在线观看| 国产亚洲av嫩草精品影院| 9色porny在线观看| 好男人电影高清在线观看| 亚洲av日韩精品久久久久久密| 三级毛片av免费| 亚洲中文日韩欧美视频| 老熟妇仑乱视频hdxx| 琪琪午夜伦伦电影理论片6080| av网站免费在线观看视频| 久99久视频精品免费| 国产aⅴ精品一区二区三区波| 久久久久久久久中文| 欧美+亚洲+日韩+国产| 亚洲av电影在线进入| 日本撒尿小便嘘嘘汇集6| 欧美日本视频| 午夜影院日韩av| 一夜夜www| 亚洲av电影不卡..在线观看| 国产精品1区2区在线观看.| 两个人免费观看高清视频| www.www免费av| 国产在线精品亚洲第一网站| 好男人电影高清在线观看| 免费高清在线观看日韩| 国产精品 欧美亚洲| 侵犯人妻中文字幕一二三四区| 欧美午夜高清在线| 欧美人与性动交α欧美精品济南到| 亚洲精品av麻豆狂野| 自线自在国产av| 极品教师在线免费播放| 亚洲第一电影网av| 日本免费一区二区三区高清不卡 | 欧美国产精品va在线观看不卡| 国产野战对白在线观看| 久久天躁狠狠躁夜夜2o2o| 欧美大码av| 99久久久亚洲精品蜜臀av| 国产一区二区三区在线臀色熟女| www.精华液| 一区二区三区激情视频| 亚洲性夜色夜夜综合| 亚洲成国产人片在线观看| 国产三级在线视频| 一级,二级,三级黄色视频| 黄片播放在线免费| 亚洲国产精品999在线| 国产精品久久久久久人妻精品电影| 午夜成年电影在线免费观看| 亚洲熟女毛片儿| 国产精品免费一区二区三区在线| 久久九九热精品免费| 啦啦啦韩国在线观看视频| 日韩欧美国产一区二区入口| 色精品久久人妻99蜜桃| 午夜日韩欧美国产| av视频免费观看在线观看| 嫁个100分男人电影在线观看| 高潮久久久久久久久久久不卡| 午夜福利高清视频| 日本免费a在线| 国产精品久久电影中文字幕| 老熟妇仑乱视频hdxx| 香蕉丝袜av| 国产成人啪精品午夜网站| 国产亚洲精品久久久久久毛片| 国产精品一区二区三区四区久久 | 多毛熟女@视频| 夜夜看夜夜爽夜夜摸| 女性生殖器流出的白浆| 国产av一区在线观看免费| 色哟哟哟哟哟哟| 亚洲av电影在线进入| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲精品一卡2卡三卡4卡5卡| 欧美激情 高清一区二区三区| 可以在线观看的亚洲视频| 国产精品亚洲美女久久久| 亚洲片人在线观看| 人人妻人人澡人人看| 欧美乱色亚洲激情| 亚洲av成人av| 99久久精品国产亚洲精品| 97人妻天天添夜夜摸| 男女床上黄色一级片免费看| 久久久久国产一级毛片高清牌| 在线观看66精品国产| 免费高清视频大片| 欧洲精品卡2卡3卡4卡5卡区| 一进一出抽搐gif免费好疼| 久久人妻熟女aⅴ| 亚洲国产精品久久男人天堂| 久久欧美精品欧美久久欧美| 高清毛片免费观看视频网站| 中文字幕精品免费在线观看视频| 国产黄a三级三级三级人| 狂野欧美激情性xxxx| 国产亚洲av高清不卡| 国产精品九九99| 亚洲第一欧美日韩一区二区三区| 不卡一级毛片| 日韩欧美免费精品| 亚洲欧美日韩高清在线视频| 天天躁狠狠躁夜夜躁狠狠躁| 99在线视频只有这里精品首页| 国产免费男女视频| 久久久水蜜桃国产精品网| 亚洲免费av在线视频| 亚洲欧美精品综合一区二区三区| 亚洲熟妇熟女久久| 多毛熟女@视频| 纯流量卡能插随身wifi吗| 制服丝袜大香蕉在线| 最新美女视频免费是黄的| 人人妻,人人澡人人爽秒播| 免费搜索国产男女视频| 精品高清国产在线一区| 久久中文看片网| 两人在一起打扑克的视频| 女同久久另类99精品国产91| 99国产综合亚洲精品| 丰满人妻熟妇乱又伦精品不卡| 午夜视频精品福利| 中亚洲国语对白在线视频| 日日夜夜操网爽| 美女免费视频网站| 18禁黄网站禁片午夜丰满| 色老头精品视频在线观看| 一二三四社区在线视频社区8| 精品久久蜜臀av无| 国产一区二区三区视频了| 性色av乱码一区二区三区2| 黄片播放在线免费| 免费少妇av软件| 如日韩欧美国产精品一区二区三区| 精品卡一卡二卡四卡免费| 免费av毛片视频| 无人区码免费观看不卡| 久久香蕉精品热| 免费高清在线观看日韩| 久久青草综合色| 国产又爽黄色视频| 性少妇av在线| 亚洲中文字幕日韩| 夜夜爽天天搞| 中文字幕高清在线视频| 亚洲激情在线av| 久久久久久亚洲精品国产蜜桃av| 亚洲国产欧美日韩在线播放| 1024视频免费在线观看| 亚洲av片天天在线观看| 一进一出好大好爽视频| 亚洲精品中文字幕在线视频| 久久精品91无色码中文字幕| 不卡av一区二区三区| 国产一区二区三区综合在线观看| 亚洲欧美精品综合久久99| 亚洲精品国产区一区二| 麻豆av在线久日| 香蕉久久夜色| 久久久水蜜桃国产精品网| 午夜久久久在线观看| 神马国产精品三级电影在线观看 | 99国产精品一区二区蜜桃av| 一本综合久久免费| 亚洲无线在线观看| 亚洲欧美日韩高清在线视频| 国产91精品成人一区二区三区| 女人被狂操c到高潮| 亚洲国产欧美日韩在线播放| 午夜福利视频1000在线观看 | 成人手机av| 在线观看舔阴道视频| 88av欧美| 精品久久久久久成人av| 午夜影院日韩av| 麻豆国产av国片精品| 国产精品99久久99久久久不卡| 免费高清在线观看日韩| 免费一级毛片在线播放高清视频 | 亚洲午夜理论影院| 亚洲精品粉嫩美女一区| 我的亚洲天堂| 两个人免费观看高清视频| 在线观看日韩欧美| 国产亚洲精品综合一区在线观看 | 成年版毛片免费区| 国语自产精品视频在线第100页| 国产av在哪里看| 国产一区二区在线av高清观看| √禁漫天堂资源中文www| 97碰自拍视频| 国产一卡二卡三卡精品| 中文字幕人妻熟女乱码| 一级毛片高清免费大全| 亚洲国产高清在线一区二区三 | 91九色精品人成在线观看| 久久精品国产清高在天天线| 色综合欧美亚洲国产小说| 身体一侧抽搐| 久久精品91无色码中文字幕| 国产精品99久久99久久久不卡| 欧美日韩亚洲综合一区二区三区_| 夜夜躁狠狠躁天天躁| 久久婷婷人人爽人人干人人爱 | 在线免费观看的www视频| 可以免费在线观看a视频的电影网站| 欧美日韩瑟瑟在线播放| 日本免费一区二区三区高清不卡 | 在线观看舔阴道视频| 中文字幕人妻丝袜一区二区| 欧美大码av| 国产精品香港三级国产av潘金莲| 亚洲第一欧美日韩一区二区三区| 亚洲自拍偷在线| 人妻丰满熟妇av一区二区三区| 啦啦啦免费观看视频1| 无人区码免费观看不卡| 999精品在线视频| 亚洲精品久久国产高清桃花| 国产精品99久久99久久久不卡| 国产精品自产拍在线观看55亚洲| 久久久久久亚洲精品国产蜜桃av| 波多野结衣巨乳人妻| 亚洲专区国产一区二区| 久久天堂一区二区三区四区| 一级,二级,三级黄色视频| 级片在线观看| 91麻豆av在线| 成人国产一区最新在线观看| 人人妻人人澡人人看| 不卡一级毛片| 91国产中文字幕| 久热爱精品视频在线9| 满18在线观看网站| 一本大道久久a久久精品| 亚洲熟妇中文字幕五十中出| 国产精品一区二区免费欧美| 欧美黑人欧美精品刺激| 国产成人av激情在线播放| 他把我摸到了高潮在线观看| 伦理电影免费视频| 午夜免费观看网址| 亚洲第一电影网av| 亚洲视频免费观看视频| 国产人伦9x9x在线观看| 国产精品,欧美在线| 欧美另类亚洲清纯唯美| 亚洲片人在线观看| 欧美一级毛片孕妇| 亚洲伊人色综图| bbb黄色大片| 黄网站色视频无遮挡免费观看| 天堂动漫精品| 亚洲五月婷婷丁香| 夜夜躁狠狠躁天天躁| 精品电影一区二区在线| 波多野结衣一区麻豆| 亚洲成人免费电影在线观看| 一本久久中文字幕| 日本黄色视频三级网站网址| 两性夫妻黄色片| 黄网站色视频无遮挡免费观看| 中国美女看黄片| 国产精品99久久99久久久不卡| 高清黄色对白视频在线免费看| 国产亚洲精品一区二区www| 久久精品人人爽人人爽视色| 91在线观看av| av有码第一页| 一区二区三区国产精品乱码| 成人亚洲精品一区在线观看| 亚洲情色 制服丝袜| 亚洲一区二区三区色噜噜| 国产一区二区在线av高清观看| 两人在一起打扑克的视频| 女人被躁到高潮嗷嗷叫费观| 别揉我奶头~嗯~啊~动态视频| 妹子高潮喷水视频| 高清毛片免费观看视频网站| 一边摸一边做爽爽视频免费| 午夜福利视频1000在线观看 | 亚洲一卡2卡3卡4卡5卡精品中文| 少妇粗大呻吟视频| 久久久国产成人免费| 久久婷婷成人综合色麻豆| 在线观看免费日韩欧美大片| 国产成人欧美| 国产极品粉嫩免费观看在线| 久久欧美精品欧美久久欧美| 久久久久久亚洲精品国产蜜桃av| 日本欧美视频一区| 亚洲国产中文字幕在线视频| 丰满的人妻完整版| 伦理电影免费视频| 国产精品二区激情视频| 日日摸夜夜添夜夜添小说| 午夜福利,免费看| 精品久久久久久久久久免费视频| 亚洲熟妇熟女久久| 亚洲av成人一区二区三| 午夜亚洲福利在线播放| 色婷婷久久久亚洲欧美| cao死你这个sao货| 国产亚洲av高清不卡| 午夜免费激情av| 亚洲人成网站在线播放欧美日韩| 真人一进一出gif抽搐免费| 亚洲中文字幕一区二区三区有码在线看 | 欧美成人免费av一区二区三区| 国产一区在线观看成人免费| √禁漫天堂资源中文www| 久久精品国产亚洲av香蕉五月| 高清在线国产一区| 天堂动漫精品| 精品久久久久久久人妻蜜臀av | 久久久久久久久免费视频了| 国产午夜精品久久久久久| 一级毛片精品| 淫秽高清视频在线观看| 好看av亚洲va欧美ⅴa在| 欧美日韩乱码在线| 国语自产精品视频在线第100页| 亚洲精华国产精华精| 黑人巨大精品欧美一区二区mp4| 一本久久中文字幕| 午夜福利高清视频| 亚洲免费av在线视频| 亚洲成人精品中文字幕电影| www.精华液| 一进一出好大好爽视频| 中文字幕久久专区| 亚洲av电影在线进入| 国产伦一二天堂av在线观看| 免费搜索国产男女视频| 精品国产超薄肉色丝袜足j| 久9热在线精品视频| 午夜视频精品福利| 欧美丝袜亚洲另类 | 久久久水蜜桃国产精品网| 久久久久精品国产欧美久久久| 久久天堂一区二区三区四区| 欧美最黄视频在线播放免费| 久久久精品欧美日韩精品| 久久精品国产综合久久久| 日本免费a在线| 美女国产高潮福利片在线看| 1024视频免费在线观看| 欧美日韩黄片免| 日本一区二区免费在线视频| 村上凉子中文字幕在线| 一卡2卡三卡四卡精品乱码亚洲| 国产蜜桃级精品一区二区三区| 精品日产1卡2卡| 国产av在哪里看| 午夜免费激情av| 亚洲va日本ⅴa欧美va伊人久久| 亚洲第一青青草原| av电影中文网址| 欧美日本中文国产一区发布| 韩国av一区二区三区四区| 久久中文字幕一级| 国语自产精品视频在线第100页| 嫁个100分男人电影在线观看| 99久久久亚洲精品蜜臀av| 免费在线观看亚洲国产| 97人妻精品一区二区三区麻豆 | 日本在线视频免费播放| 久久精品影院6| 91在线观看av| av网站免费在线观看视频| 一区二区三区激情视频| 黄色片一级片一级黄色片| 精品久久久久久,| 欧美日韩亚洲国产一区二区在线观看| 亚洲激情在线av| 日日夜夜操网爽| 叶爱在线成人免费视频播放| 国产精品98久久久久久宅男小说| 日日摸夜夜添夜夜添小说| 亚洲av电影在线进入| 欧美精品啪啪一区二区三区| 午夜福利视频1000在线观看 | 自拍欧美九色日韩亚洲蝌蚪91| 精品国产美女av久久久久小说| 日韩大尺度精品在线看网址 | 在线观看66精品国产| 真人一进一出gif抽搐免费| av天堂在线播放| 一卡2卡三卡四卡精品乱码亚洲| 亚洲,欧美精品.| 日日摸夜夜添夜夜添小说| 免费观看人在逋| 一级黄色大片毛片| 中文字幕人妻熟女乱码| 一级a爱片免费观看的视频| 国产片内射在线| 亚洲精品粉嫩美女一区| 一本大道久久a久久精品| 麻豆av在线久日| 午夜视频精品福利| 大码成人一级视频| 国产成年人精品一区二区| 麻豆一二三区av精品| а√天堂www在线а√下载| 黑人巨大精品欧美一区二区mp4| 久99久视频精品免费| 国产精品久久久av美女十八| 两个人看的免费小视频| or卡值多少钱| 一级a爱视频在线免费观看| 如日韩欧美国产精品一区二区三区| 久久久久亚洲av毛片大全| 午夜福利成人在线免费观看| 一二三四社区在线视频社区8| 香蕉久久夜色| 午夜免费激情av| 日本黄色视频三级网站网址| 99riav亚洲国产免费| 精品少妇一区二区三区视频日本电影| 日韩高清综合在线| 一级黄色大片毛片| 啪啪无遮挡十八禁网站| 好看av亚洲va欧美ⅴa在| 国产成人精品在线电影| 人妻丰满熟妇av一区二区三区| 欧美黑人欧美精品刺激| 搡老熟女国产l中国老女人| 国产高清有码在线观看视频 | 欧美激情 高清一区二区三区| 欧美一级毛片孕妇| 亚洲美女黄片视频| 人人妻人人澡人人看| 嫩草影视91久久| 亚洲人成伊人成综合网2020| 此物有八面人人有两片| 69精品国产乱码久久久| 亚洲一区二区三区色噜噜| 国内精品久久久久久久电影| 日韩欧美一区视频在线观看| 国产精品电影一区二区三区| 天堂√8在线中文| xxx96com| 一区福利在线观看| 亚洲av电影不卡..在线观看| 亚洲第一av免费看| 欧美色欧美亚洲另类二区 | tocl精华| 在线观看www视频免费| 精品一区二区三区四区五区乱码| 亚洲精品久久成人aⅴ小说| 日日夜夜操网爽| 久久久久久久精品吃奶| 纯流量卡能插随身wifi吗| 一卡2卡三卡四卡精品乱码亚洲| АⅤ资源中文在线天堂| 国产精品国产高清国产av| 99精品久久久久人妻精品| 国产又色又爽无遮挡免费看| 啦啦啦韩国在线观看视频| 亚洲一码二码三码区别大吗| 亚洲精品中文字幕在线视频| 亚洲欧美精品综合久久99| 亚洲av成人av| 久久精品亚洲熟妇少妇任你| 一区在线观看完整版| 中文字幕人成人乱码亚洲影| 日日爽夜夜爽网站| av超薄肉色丝袜交足视频| 黄色 视频免费看| 精品欧美一区二区三区在线| 欧美日本中文国产一区发布| 曰老女人黄片| 国产成人一区二区三区免费视频网站| 久久亚洲精品不卡| 麻豆一二三区av精品| 国产精品久久视频播放| 国产男靠女视频免费网站| 正在播放国产对白刺激| 叶爱在线成人免费视频播放| 妹子高潮喷水视频| 成人三级做爰电影| 亚洲中文av在线| 亚洲一区中文字幕在线| 黄色片一级片一级黄色片| 丰满的人妻完整版| av电影中文网址| 香蕉国产在线看| 一a级毛片在线观看| 精品久久久久久久人妻蜜臀av | 亚洲最大成人中文| 可以在线观看的亚洲视频| 国产av一区二区精品久久| 久久精品国产综合久久久| 禁无遮挡网站| 巨乳人妻的诱惑在线观看| 国产99久久九九免费精品| 人成视频在线观看免费观看| 黑人巨大精品欧美一区二区蜜桃| www.熟女人妻精品国产| 精品国产乱码久久久久久男人| 色哟哟哟哟哟哟| 亚洲国产中文字幕在线视频| 精品久久久久久久毛片微露脸| 黑人巨大精品欧美一区二区mp4| 成人国语在线视频| 国产一区二区三区综合在线观看| 国产精品一区二区精品视频观看| 97超级碰碰碰精品色视频在线观看| 国产成年人精品一区二区| 色老头精品视频在线观看| 久久这里只有精品19| 1024香蕉在线观看| 免费在线观看视频国产中文字幕亚洲| 别揉我奶头~嗯~啊~动态视频| av视频免费观看在线观看| av电影中文网址| 亚洲一卡2卡3卡4卡5卡精品中文| av视频免费观看在线观看| 国产伦一二天堂av在线观看| 一级作爱视频免费观看| 色综合亚洲欧美另类图片| 国产免费男女视频| 可以在线观看的亚洲视频| 久久精品国产清高在天天线| 老司机深夜福利视频在线观看| 免费在线观看黄色视频的| 91麻豆精品激情在线观看国产| 不卡av一区二区三区| 黄色片一级片一级黄色片| 日韩欧美三级三区| 一边摸一边抽搐一进一出视频| bbb黄色大片| 亚洲男人天堂网一区| 老熟妇乱子伦视频在线观看| 亚洲欧美激情综合另类| 一二三四社区在线视频社区8| 亚洲欧美精品综合一区二区三区| 在线av久久热| 国产一区二区激情短视频|