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

    中低緯度電離層偶發(fā)E層電波傳播建模?

    2017-08-09 07:34:48郝書吉張文超張雅彬楊巨濤馬廣林
    物理學(xué)報 2017年11期
    關(guān)鍵詞:電子密度電波電離層

    郝書吉 張文超張雅彬 楊巨濤 馬廣林

    (中國電波傳播研究所,電波環(huán)境特性及模化技術(shù)重點實驗室,青島 266107)

    中低緯度電離層偶發(fā)E層電波傳播建模?

    郝書吉 張文超?張雅彬 楊巨濤 馬廣林

    (中國電波傳播研究所,電波環(huán)境特性及?;夹g(shù)重點實驗室,青島 266107)

    (2017年1月4日收到;2017年3月8日收到修改稿)

    基于電離層偶發(fā)E層(sporadic-E,Es)對電波傳播的多條鏈路實測數(shù)據(jù),認(rèn)為電離層Es對入射電波的作用是反射和散射兩種機制,且反射/散射的比例常數(shù)隨著電離層Es臨界頻率的變化而變化,進而建立了包含反射、散射以及反射與散射共同作用的三段式電離層Es電波傳播模型,并與國際電聯(lián)(International Telecommunications Union,ITU)給出的Es層傳播的電波場強預(yù)測模型進行了對比,驗證了本文所建模型的正確性.該模型尤其適用于中低緯地區(qū)甚高頻信號經(jīng)Es層的傳播研究.

    電離層Es,電波傳播,ITU模型

    1 引言

    偶發(fā)E層(sporadic-E,Es)是電離層內(nèi)一種短暫而不規(guī)則的強電離薄層,其高度范圍為90—140 km,厚度從數(shù)百米至1 km不等,水平尺度從幾十千米到幾百千米不等[1].電離層Es電子密度非常高,可達(dá)到常規(guī)E層電子密度的100倍[2].利用Es層可以實現(xiàn)單跳2000 km的超短波超視距通信,這對于建立超短波超視距應(yīng)急通信保障系統(tǒng)具有重要意義.

    研究人員從20世紀(jì)40年代就開始關(guān)注電離層Es對入射電波的影響及其信道特性,主要基于電離層探測儀,通過甚高頻(very high frequency,VHF)后向和斜入射鏈路對不同緯度的Es層進行觀測和統(tǒng)計研究.但由于觀測手段和統(tǒng)計方法的不同,結(jié)論也不盡相同.早在1964年,日本學(xué)者Kobayashi[3]通過長1100 km的VHF鏈路觀測,初步認(rèn)為電離層Es的反射損耗與電波入射角(θ)和電離層Es臨界頻率(foEs)的關(guān)系為Γ=25[f cos(θ)/foEs]2.67,其中θ為高度105 km的鏡反射入射角;1972年,Bramley[4]認(rèn)為電離層Es中VHF信號傳播的反射系數(shù)與(foEs/f)j成正比,其中j的值在7—30之間;1972年,Althouse和Davis[5]研究了Es層的散射特性,并對其后向散射進行了評估;1978年,日本的Miya和Shimizu[6]通過Tokyo-Wakkanai鏈路觀測,得到經(jīng)電離層Es的電波傳播經(jīng)驗公式;1980年,Kerblay和Makarenko[7]研究了Es層的大尺度水平不均勻性對反射系數(shù)的影響,并給出了中間點位置的反射系數(shù);1999年,國際電聯(lián)(International Telecommunications Union,ITU)在ITU-R p.1240建議書中給出了Es的臨界頻率發(fā)生概率分布圖,可以預(yù)測世界各地Es的發(fā)生;2007年,Sherstyukov和Akchurin[8]對中緯度Es層單跳傳播方式建模,用來預(yù)測Es層傳播的可用頻率分布;2012年,ITU發(fā)布了計算經(jīng)過Es層電波場強的建議書,在建議書中給出了經(jīng)過Es傳播的電離層Es衰減經(jīng)驗公式,其衰減值是入射電波頻率與電離層臨界頻率比值的二次方成正比;2015年,Sherstyukov等[9,10]通過對高頻(HF)和VHF頻段的無線電波測量數(shù)據(jù)的統(tǒng)計分析,也獲得了電離層Es的反射損耗.

    上述研究主要從反射角度并基于大量實驗觀測數(shù)據(jù)研究電離層Es對入射電波的作用,從電離層Es對電波的反射系數(shù)和傳播損耗方面進行了分析和建模研究.本文基于多條鏈路實驗數(shù)據(jù),認(rèn)為電離層Es對入射電波的作用很難用單一的反射或散射作用機制進行解釋,結(jié)合中低緯地區(qū)電離層Es電子密度分布特性及其中各向異性結(jié)構(gòu)的尺度關(guān)系,判斷電離層Es對入射電波的作用可能是散射作用和反射作用兩者機制都存在,并且反射作用的比例常數(shù)隨著電離層Es臨界頻率的變化而變化.據(jù)此,根據(jù)不同鏈路上的測量數(shù)據(jù),建立了包含反射、散射以及反射與散射共同作用的三段式電離層Es電波傳播模型,并與ITU模型進行了對比,驗證了本文所建模型的正確性.

    2 電離層Es電波傳播理論

    理論和實驗研究表明:電離層Es電子密度非常高,會對入射的VHF波段信號產(chǎn)生反射作用,同時Es層內(nèi)部存在大量電子密度不規(guī)則體,會對入射的信號產(chǎn)生散射作用.

    2.1 反射理論

    理論和試驗研究表明[2,11]:電離層Es層電子密度梯度在垂直方向上相對比較陡峭,據(jù)此特點,可假定Es層是由水平方向上為有限尺度的、垂直方向上具有不同介電常數(shù)的無數(shù)個薄介電層組成,針對每一個薄介電層,其VHF信號經(jīng)過Es層的反射示意圖見圖1.其中d為發(fā)射點和接收點的距離;θ為掠射角;a和b為兩個方向上的水平尺度,(λ是入射波波長).若相鄰層的介電常數(shù)的增量為dε,且滿足1?θ2?dε,則薄層間的反射系數(shù)為dρ=dε/(4θ2).

    圖1 薄層反射幾何示意Fig.1.Geometry of a thin layer model.

    目前,普遍認(rèn)為“風(fēng)剪切”理論是中緯度Es層形成的物理機制.設(shè)Es層電子密度在“風(fēng)剪切”的作用下其垂直分布滿足如下分布[12]:

    其中,Nmax是電離層Es最大電子密度;?h是Es層半厚度;z′=z?h,h為Es的中心高度,|z′|6?h;n=1,2,3···.通常Es層介電常數(shù)ε可寫為[13]

    其中fN是Es層最大等離子體頻率,fN=√是自由空間介電常數(shù),f是入射波頻率,w=2πf;e是電子電量,me是電子質(zhì)量.根據(jù)反射理論,總反射系數(shù)為

    其中λ為入射波波長;L=4πsinθ?h/λ.根據(jù)雷達(dá)方程,則經(jīng)過Es層反射的VHF信號接收功率為

    其中Pt為發(fā)射功率,Gt和Gr分別為發(fā)射和接收天線在指向Es高度方向的增益,D為傳播路徑.

    2.2 散射理論

    與對流層或常規(guī)F層不規(guī)則體相比[14,15],電離層Es不規(guī)則體尺度體要大得多[16,17],因此可運用其電子密度自相關(guān)函數(shù)來分析其散射現(xiàn)象[18].一般Es層電子密度自相關(guān)函數(shù)為

    其中:C(0)是大氣折射系數(shù)的時間方差,在Es中C(0)=(λ/λN)4(?N/N)/4;Γ為伽馬函數(shù),γ為Es不規(guī)則體尺度,l為一般湍流的尺度,p為常數(shù)因子,Km(p)為二類高階修正貝塞爾函數(shù).由于電離層Es電子密度通常比一般散射體的要大,用高階修正貝塞爾函數(shù)可更好地反映其散射特性[19].

    根據(jù)散射理論[20],電離層Es不規(guī)則體散射截面積為

    其中q=2klsinθ;χ為散射方向與入射電場方向的夾角;L1,L2是Es層不規(guī)則體在水平方向上平行、正交于水平漂移的不同尺度,L3是Es不規(guī)則體在垂直方向的尺度.

    根據(jù)(5)式,Es層不規(guī)則體尺度可表示為γ=pl,p由Es層電子密度自相關(guān)函數(shù)和修正貝塞爾函數(shù)的階數(shù)決定,通常電子密度自相關(guān)函數(shù)[21]取值為0.5,按照中低緯地區(qū)的經(jīng)驗取值,修正貝塞爾函數(shù)的階數(shù)m取5[12],由(5)式可確定出p=3.5,γ=3.5l.則電離層Es不規(guī)則體散射截面進一步表示為

    根據(jù)雷達(dá)方程,則經(jīng)電離層Es散射的VHF信號接收功率為

    3 電離層Es電波傳播模型

    為了深入研究中緯度地區(qū)電離層Es層VHF傳播特性,選擇在夏季Es高發(fā)的時段開展了電離層Es層VHF信號測量試驗,實驗共進行了20余天,有效實驗數(shù)據(jù)約100 h.為覆蓋不同的區(qū)域和距離,試驗采用南北和東西兩種鏈路布局.其中南北布局為昆明-西安的鏈路1,距離為1065 km;東西布局有兩條鏈路,其中德宏-懷化的鏈路2,距離為1240 km,德宏-郴州的鏈路3,距離為1590 km.試驗使用了多個測量頻率,并在每條鏈路中重點使用了一個頻率,分別為50.5,45.5,31 MHz.測量鏈路布局如圖2所示.

    實驗期間典型VHF信號接收功率以及對應(yīng)的頻高圖如圖3所示.由圖3可以看出:在無Es期間VHF接收信號大多淹沒在噪聲中,且信號起伏較大,不穩(wěn)定;隨著電離層Es由弱變強,接收信號逐漸增強.因此判定正常E區(qū)以及F區(qū)不規(guī)則結(jié)構(gòu)相對突發(fā)Es層的散射能力很弱,本文有效數(shù)據(jù)要求:1)明顯出現(xiàn)Es且持續(xù)至少1 min時段;2)信號至少高于背景噪聲10 dB.同時參照信道時延功率譜將較強的F層不均勻體散射的數(shù)據(jù)剔除.經(jīng)上述數(shù)據(jù)篩選后,三條測量鏈路的VHF信號平均接收功率隨著Es層臨界頻率與入射頻率的比值foEs/f的變化如圖4所示.

    由圖4可以看出:在三條測量鏈路中VHF信號平均接收功率隨著foEs/f的增長而上升.但VHF信號接收功率低端增長緩慢,中間段快速上升,高端增長又變緩.利用第2節(jié)給出的反射和散射理論對鏈路1進行仿真,根據(jù)電離層Es的厚度通常在幾百米到1 km,其不規(guī)則體水平方向尺度在幾十米到幾百米,垂直方向尺度幾十米,選擇反射作用機制中的n=5,電離層Es半厚度?h取50和500 m,散射作用機制中垂直尺度L3取30 m,水平方向尺度L1取50和300 m,仿真結(jié)果如圖5所示.

    圖2 實驗測量鏈路布局Fig.2.Link-arrangement of test and measurement.

    圖3 (網(wǎng)刊彩色)實驗期間典型測量結(jié)果(a)VHF接收信號功率;(b)鏈路中點垂測儀觀測結(jié)果Fig.3.(color online)Typical measurement results during the test:(a)Received power of VHF signal;(b)results of ionosonode in midpoint of measurement-link.

    圖5 VHF信號接收功率仿真Fig.5.Simulation of the received power of VHF signal.

    由圖5可以看出,在foEs/f高端,VHF接收信號平均功率與反射作用機制相符,在低端與散射機制相符.而在中間過渡區(qū),VHF信號接收功率快速上升,其斜率和foEs/f較高和較低區(qū)域明顯不同,即在該區(qū)域已經(jīng)不能用單一的反射作用或散射作用機制解釋,據(jù)此可初步判斷在該區(qū)域可能是反射與散射的共同作用機制,并且隨著foEs/f的上升,反射作用機制逐漸增強.

    下面通過實測數(shù)據(jù)并結(jié)合電離層Es對無線電波的反射和散射作用機制,研究并建立電離層Es電波傳播模型.首先定義實測數(shù)據(jù)與仿真結(jié)果(基于反射作用機制或散射作用機制)之間的平均誤差ˉδ為

    其中,N是在考察區(qū)域內(nèi)樣本點的個數(shù);Pr,P′r是實測接收功率和仿真接收功率.為均衡不同鏈路中的誤差,定義本次測量試驗中三條鏈路實測數(shù)據(jù)與仿真結(jié)果之間的總平均誤差ˉδ總為

    3.1 反射作用分析

    根據(jù)(3)式,對于入射到電離層Es單一頻率電波,在收發(fā)兩點距離、Es層高度確定后,影響反射信號強度的主要是Es層的半厚度?h和電離層Es的電子密度分布指數(shù)n.當(dāng)n相同時,?h越大,電波衰減越大,接收功率越小.事實上,相同電離層Es反射系數(shù)具有多組n和?h的組合,本文按照中低緯地區(qū)的經(jīng)驗[22]取值n=5,通過實驗數(shù)據(jù)確定等效半厚度?h.根據(jù)前面論述,經(jīng)Es傳播的VHF接收信號在foEs/f高端主要是Es反射作用,設(shè)這一區(qū)間為foEs/f>fr,fr為完全可用反射作用機制解釋的門限值.圖6給出了不同等效半厚度下,三條鏈路經(jīng)電離層Es反射的實測結(jié)果與基于反射作用機制仿真結(jié)果的平均誤差δˉ隨fr的變化,表1給出了三條鏈路中均滿足實測數(shù)據(jù)與仿真結(jié)果δˉ<2%時的不同鏈路門限值fr和總平均誤差δˉ總.

    通常δˉ<2%在實際測量中已經(jīng)是比較苛刻的要求,由圖6可以看出滿足此要求的條件為:鏈路1中90 m

    圖6 (網(wǎng)刊彩色)不同等效半厚度下平均誤差ˉδ隨fr的變化(a)鏈路1;(b)鏈路2;(c)鏈路3Fig.6.(color online)Average errorˉδ under di ff erent half-thickness changed with fr:(a)Link 1;(b)Link 2;(c)Link 3.

    表1 實測數(shù)據(jù)與基于反射作用仿真結(jié)果誤差ˉδ均小于2%的門限值與平均誤差Table 1.The threshold and total average error under the average errorˉδ between measured power and simulated power based on re fl ecting process less than 2%.

    如果在整個電離層Es持續(xù)時間內(nèi)?h變化不大,上述結(jié)果表明:某一頻率的VHF電波入射Es層,當(dāng)foEs/f達(dá)到某個門限值時,接收端接收到的信號主要來自于Es層對VHF電波的反射.當(dāng)foEs/f大于此門限值時,可用電離層Es反射作用機制進行解釋,此時Es反射系數(shù)等效表示為

    當(dāng)foEs/f低于該值時,需尋求其他機制去解釋.

    3.2 散射作用分析

    根據(jù)(10)式,電離層Es中不規(guī)則體水平尺度相同時,垂直方向上尺度越大,衰減越大;而當(dāng)垂直方向上尺度相同時,水平尺度越大,衰減越小.根據(jù)前面的論述,經(jīng)Es傳播的VHF接收信號在foEs/f低端主要是Es散射作用,設(shè)這一區(qū)間為foEs/f

    根據(jù)3.1節(jié)中δˉ<2%的誤差要求,由圖7可以看出滿足此要求的條件為:鏈路1中160 m

    表2 實測數(shù)據(jù)與基于散射作用仿真結(jié)果誤差ˉδ均小于2%的門限值與總平均誤差Table 2.The threshold and total average error under the average errorˉδ between measured power and simulated power based on scattering process less than 2%.

    圖7 (網(wǎng)刊彩色)不同水平尺度下平均誤差ˉδ隨fs的變化(a)鏈路1;(b)鏈路2;(c)鏈路3Fig.7.(color online)Average errorˉδ under di ff erent scale changed with fs:(a)Link 1;(b)Link 2;(c)Link 3.

    如果在整個電離層Es持續(xù)時間內(nèi)L1,L2,L3變化不大,上述結(jié)果表明:當(dāng)某一頻率的VHF電波入射Es層,當(dāng)foEs/f低于某個門限值時,接收端接收到的信號主要來自于Es層對VHF電波的散射,此時Es散射截面積可以等效表示為

    當(dāng)foEs/f大于該門限值時,需尋求其他機制去解釋.

    3.3 電離層Es電波傳播模型

    綜合分析3.1節(jié)和3.2節(jié)仿真結(jié)果,發(fā)現(xiàn)在鏈路1中當(dāng)foEs/f在0.11—0.14區(qū)間時,鏈路2中foEs/f在0.14—0.21區(qū)間時,鏈路3中foEs/f在0.13—0.25區(qū)間時,利用反射作用機制和散射作用機制都不能有效解釋實測結(jié)果.根據(jù)前面的分析,電離層Es對入射的電波既有反射作用機制,又有散射作用機制,兩種作用隨電離層的背景條件而發(fā)生變化.設(shè)反射作用機制與散射作用機制比例常數(shù)為kr,ks,且kr+ks=1,則此時接收端VHF信號強度Pr可表示為

    其中,Pt是VHF信號發(fā)射功率;Hr,Hs為對應(yīng)Es層對VHF信號反射作用和散射作用機制的傳輸參數(shù).

    根據(jù)上述分析,此區(qū)間即為fs

    圖8 (網(wǎng)刊彩色)過渡區(qū)域內(nèi)反射作用所占比重的擬合結(jié)果Fig.8.(color online)The fi tting result of coefficient re fl ection process in transition section.

    由圖8可以看出,在此區(qū)間內(nèi)反射作用所占比重隨著f′的上升而逐漸增大,并且反射作用所占比重更符合指數(shù)分布,

    確立過渡區(qū)域反射作用和散射作用所占比重后,整個Es存在期間內(nèi)電離層Es對電波的作用隨foEs/f變化的作用機制即確定.根據(jù)上面的分析,將(11),(12)和(14)式代入(13)式,建立由散射作用、散射與反射共同作用、反射作用三段式電離層Es電波傳播模型,表示為

    其中:Pt為發(fā)射功率;Gt和Gr分別為發(fā)射和接收天線在指向Es高度方向的增益;D為傳播路徑;λ是入射波波長;θ為掠射角,L=460πsinθ/λ.kr是反射信號所占比重,符合如下分布:

    其中:fr是完全反射作用機制解釋的門限值,fs是完全散射作用機制解釋的門限值.在本次測量中,鏈路1中fr=0.14,fs=0.11;鏈路2中fr=0.21,fs=0.14;鏈路3中fr=0.25,fs=0.13.根據(jù)該模型的(15)式和(16)式對三條鏈路進行模擬,結(jié)果如圖9所示.

    由圖9可以看出,本文給出的三段式動態(tài)模型很好地解釋了實驗期間接收端VHF信號幅度的變化:foEs/f>fr區(qū)間,VHF信號傳播符合電離層Es反射作用機制;foEs/f

    圖9 (網(wǎng)刊彩色)基于本文Es電波傳播模型的仿真結(jié)果和實測數(shù)據(jù)Fig.9.(color online)The result of simulated on the proposed model of radio wave propagation under Es and the result of measurement.

    3.4 模型對比

    根據(jù)本文所建立的模型,對三條鏈路中電離層Es對入射電波的衰減進行仿真,并與ITU模型比較.設(shè)Es高度110 km,鏈路1中f=50.5 MHz,fr=0.14,fs=0.11;鏈路2中f=45.5 MHz,fr=0.21,fs=0.14;鏈路3中f=31.5 MHz,fr=0.25,fs=0.13.仿真結(jié)果如圖10所示.

    利用本文模型計算的VHF信號衰減值與ITU模型結(jié)果一致,特別是在foEs/f高端時,一致程度更好;但在foEs/f取較小的值時,利用本文模型計算的VHF信號衰減值要小很多.分析其原因,本文所建模型是針對中低緯地區(qū)電離層Es對電波作用的實測數(shù)據(jù),電離層Es特性電離層背景關(guān)系比較密切,不同緯度地區(qū)電離層Es特性有差別,而ITU模型是基于中高緯地區(qū)的實測數(shù)據(jù)而建立,其應(yīng)用在中低緯緯度誤差會增大;其次,ITU模型預(yù)測時采用實驗數(shù)據(jù)統(tǒng)計后的經(jīng)驗值,與實測結(jié)果也存在差異,因此這種差異會導(dǎo)致誤差;最后,ITU預(yù)測模型誤差會隨著foEs/f減小而增大,在特別是對于foEs/f取較小值時,電離層Es對入射電波主要是散射作用機制,而散射作用遠(yuǎn)小于反射作用,進一步增大了誤差.因此,我們認(rèn)為本文所建更適用于中低緯地區(qū).

    圖10 VHF信號接收功率隨foEs/f變化(a)鏈路1;(b)鏈路2;(c)鏈路3Fig.10.Received power of VHF signal changed with foEs/f:(a)Link 1;(b)Link 2;(c)Link 3.

    4 結(jié)論

    本文基于實測數(shù)據(jù),結(jié)合電離層Es電子密度分布特點,建立包含散射、反射與散射共同作用以及反射的三段式動態(tài)電離層Es電波傳播模型,該模型中電離層Es對入射電波的作用隨foEs/f變化如下:

    1)當(dāng)foEs/f高于某一門限(完全可用反射作用機制解釋的門限值fr)時,VHF信號經(jīng)Es層的傳播符合反射作用機制;

    2)當(dāng)foEs/f低于某一門限(完全可用散射作用機制解釋的門限值fs)時,VHF信號經(jīng)Es層的傳播符合散射作用機制;

    3)當(dāng)foEs/f在fs與fr之間的過渡區(qū)域時,VHF信號經(jīng)Es層的傳播是反射和散射共同作用機制,且此區(qū)間隨著foEs/f增大,反射作用機制逐漸增強.

    對三條鏈路實驗數(shù)據(jù)統(tǒng)計可發(fā)現(xiàn),不同鏈路間fr和fs的值存在差異,尤其是鏈路1中fr,fs與鏈路2鏈路3中的fr,fs差異較大,其可能原因是鏈路1與鏈路2、鏈路3地理位置中心點距離較遠(yuǎn),電離層背景參數(shù)差異較大,如E層或F2層的散射強度等;同時由于采用的探測波頻率也不相同,也可能導(dǎo)致各鏈路模型參數(shù)存在差異.

    由于本文在模型建立過程中基于的實驗數(shù)據(jù)相對較少,模型參數(shù)確認(rèn)與優(yōu)化還需要后期通過大量實驗數(shù)據(jù)完成.但本文所建模型是依據(jù)中低緯度地區(qū)測量數(shù)據(jù),在后期模型優(yōu)化后可以用來預(yù)測VHF頻段信號經(jīng)Es傳播能量分布,具有較強的工程應(yīng)用價值.

    [1]Smith L G,Mechtly E A 1972 Radio Sci.7 367

    [2]Beynon W J G,Maude A D 1972 Planet.Space Sci.20 809

    [3]Kobayashi T 1964 Radio Res.11 181

    [4]Bramley E N 1972 J.Atmos.Terr.Phys.34 1495

    [5]Althouse E L,Davis J R 1972 Radio Sci.7 897

    [6]Miya K,Shimizu K 1978 Radio Sci.13 559

    [7]Kerblay T S,Makarenko S F 1980 Geomagn.Aeron.20 449

    [8]Sherstyukov O N,Akchurin A D 2009 Adv.Space Res.43 1835

    [9]Sherstyukov O N,Akchurin A D 2011 General Assembly and Scienti fi c SymposiumIstanbul,Aug.13–20,2011 p1

    [10]Sherstyukov O N,Akchurin A D,Sherstyukov R O 2015 Adv.Space Res.56 1169

    [11]Whitehead J D 1989 J.Atmosph.Solar-Terr.Phys.51 401

    [12]Tao K 1962 Ionospheric Sporadic E(New York:MacMillian Company)p235

    [13]Pan W Y 1981 Acta Phys.Sin.30 661(in Chinese)[潘威炎1981物理學(xué)報30 661]

    [14]Zeng Z,Sokolovskiy S 2010 Geophys.Res.Lett.37 1480

    [15]Huang C S,Li J 1994 Acta Phys.Sin.43 1476(in Chinese)[黃朝松,李鈞1994物理學(xué)報43 1476]

    [16]Zhang Q,Wu X J 2016 Acta Phys.Sin.65 038102(in Chinese)[張卿,武新軍2016物理學(xué)報65 038102]

    [17]Deng F,Zhao Z Y,Shi R,Zhang Y N 2009 Acta Phys.Sin.58 7382(in Chinese)[鄧峰,趙正予,石潤,張援農(nóng)2009物理學(xué)報58 7382]

    [18]Booker H G,Gordon W E 1950 Proceed.Ire.38 401

    [19]Norton K A 1956 Commun.Syst.Ire Trans.4 39

    [20]Malaga A 1986 Rome Air Development Center Air Force Systems Command Griffiss Air Force Base New York,Aug 1985 p13441

    [21]Booker H G,Gordon W E 1958 Proceed.Ire.46 298

    [22]Ovezgeldyev O G 1988 Geomagn.Aeron.28 1024

    (Received 4 January 2017;revised manuscript received 8 March 2017)

    PACS:94.20.Vv,94.30.TzDOI:10.7498/aps.66.119401

    F信號接收功率隨foEs/f的變化 Fig.4.

    power of VHF signal changed with foEs/f.

    Modeling of radio wave propagations under sporadic-E in fl uence at low and middle latitudes?

    Hao Shu-JiZhang Wen-Chao?Zhang Ya-BinYang Ju-TaoMa Guang-Lin
    (National Key Laboratory of Electromagnetic Environment,China Research Institute of Radiowave Propagation,Qingdao 266107,China)

    The sporadic-E(Es)layer is a thin layer of several kilometers existing at an altitude around 100 km and features extremely dense ionized irregularities,which can re fl ect or scatter high frequency(HF)and very high frequency(VHF)radio waves.The most popular theoretical explanation for mid-latitude Es formation is the wind shear theory.Measurements by rocket souding have shown that Es has high electron density and relatively sharp density gradient in the vertical direction.The one-hop propagation of VHF signal in Es can even reach as far as 2000 km.In this paper,we consider incident radio waves in fl uenced by Es via both re fl ecting and scattering processes at low and middle latitudes,the coefficients of which are related to and vary with the critical frequency of Es(foEs).Firstly,with a supposed parabolic density distribution and the autocorrelation function of the electron density given by Booker,HF and VHF radio wave propagations in Es are analyzed according to the re fl ection and scattering theory.Secondly,a numerical model for the combined re fl ecting and scattering processes is developed in the form of piecewise function,the contribution of which can be distinguished by the portion factor of re fl ection(kr).According to the model,there are two threshold ratios of the critical frequency to the wave frequencies frand fsrespectively.The incident radio waves are totally re fl ected by Es when foEs/f is higher than frand mostly scattered when foEs/f is lower than fs.A transition zone exists between two critical points,with the combined processes working together.Thirdly,HF/VHF radio wave propagations in low and middle latitudes of Es are are in the north-southern direction and east-western direction separately.The experiment link in the north-southern direction is from Kunming to Xi’an at distance of 1065 km,and the ionosonde used for Es observation is located at Chongqing.Two east-west links are arranged,one of which is from Dehong to Huaihua and the other is from Dehong to Chenzhou,with the ionosonde located at Kunming and the ground distance as far as 1240 km and 1590 km respectively.The measurement data are treated and parameters of the above mentioned model for wave propagation in Es are experimentally determined.Finally,our model is veri fi ed by comparing with ITU-R model.Our results are consistent with the results from the ITU-R model when the foEs/f is high(i.e.,the re fl ecting process plays a main role).When the scattering process dominates,the attenuation value of VHF signal is far less than that predicted by the ITU-R model,which is closer to actual measurements.It is concluded that our model is more preferable for HF and VHF radio wave propagations in Es at low and middle latitudes.

    sporadic-E,radio wave propagation,International Telecommunications Union model

    10.7498/aps.66.119401

    ?中國電科技術(shù)創(chuàng)新基金(批準(zhǔn)號:A171601C01)資助的課題.

    ?通信作者.E-mail:zhangwenchao0726@126.com

    ?2017中國物理學(xué)會Chinese Physical Society

    http://wulixb.iphy.ac.cn

    *Project supported by the Technology Innovate Fund of China Electronics Technology Group Corporation(Grant No.A171601C01).

    ?Corresponding author.E-mail:zhangwenchao0726@126.com

    猜你喜歡
    電子密度電波電離層
    永遠(yuǎn)的紅色電波
    中國記者(2021年11期)2021-12-11 09:05:04
    一種電離層TEC格點預(yù)測模型
    The Speed of Light
    Kalman濾波估算電離層延遲的一種優(yōu)化方法
    顧及地磁影響的GNSS電離層層析不等像素間距算法*
    不同GPS掩星電離層剖面產(chǎn)品相關(guān)性分析
    測繪通報(2019年11期)2019-12-03 01:47:34
    瞌睡電波
    幽默大師(2019年11期)2019-11-23 08:47:58
    等離子體電子密度分布信息提取方法研究
    一種適用于電離層電子密度重構(gòu)的AMART算法
    “電波衛(wèi)士”在行動
    成人一区二区视频在线观看| 在线免费十八禁| 欧美性猛交黑人性爽| 亚洲五月天丁香| 国产视频内射| 亚洲无线观看免费| 深夜a级毛片| 欧美性猛交黑人性爽| 久9热在线精品视频| 色av中文字幕| 内射极品少妇av片p| 午夜福利在线观看免费完整高清在 | 最近最新免费中文字幕在线| 九九爱精品视频在线观看| 熟妇人妻久久中文字幕3abv| 国产真实伦视频高清在线观看 | 国产精品国产高清国产av| 精品人妻1区二区| 免费高清视频大片| 麻豆国产av国片精品| 天堂影院成人在线观看| 白带黄色成豆腐渣| 国产伦精品一区二区三区视频9| 日本免费a在线| 亚洲人成伊人成综合网2020| 赤兔流量卡办理| 欧美日韩黄片免| 男人舔女人下体高潮全视频| 综合色av麻豆| 丰满乱子伦码专区| 亚洲综合色惰| 神马国产精品三级电影在线观看| 如何舔出高潮| 色av中文字幕| 亚洲精品456在线播放app | 国内少妇人妻偷人精品xxx网站| 日本黄色视频三级网站网址| 麻豆久久精品国产亚洲av| av福利片在线观看| 麻豆久久精品国产亚洲av| 国产高潮美女av| 欧美日韩瑟瑟在线播放| 99久久九九国产精品国产免费| 国产亚洲精品av在线| 韩国av在线不卡| 国产一区二区激情短视频| 亚洲经典国产精华液单| 亚洲欧美日韩高清在线视频| 免费看日本二区| 亚洲av免费高清在线观看| 亚洲国产欧美人成| 国产精品无大码| 欧美三级亚洲精品| 欧美精品啪啪一区二区三区| 淫妇啪啪啪对白视频| 精品99又大又爽又粗少妇毛片 | 搞女人的毛片| 免费av不卡在线播放| 国产精品国产高清国产av| 欧美一区二区精品小视频在线| 一a级毛片在线观看| 九九热线精品视视频播放| 99久久精品国产国产毛片| 婷婷丁香在线五月| 伦精品一区二区三区| 欧美日韩黄片免| 搡女人真爽免费视频火全软件 | 亚洲国产精品成人综合色| 欧美xxxx性猛交bbbb| 黄色欧美视频在线观看| 成人毛片a级毛片在线播放| 九九久久精品国产亚洲av麻豆| 一级毛片久久久久久久久女| 美女高潮的动态| 久久久久久久久中文| 日本黄色视频三级网站网址| 亚洲熟妇中文字幕五十中出| 免费观看的影片在线观看| 亚洲图色成人| 99视频精品全部免费 在线| 国产精品无大码| 91久久精品国产一区二区三区| 嫁个100分男人电影在线观看| 一卡2卡三卡四卡精品乱码亚洲| 18+在线观看网站| 亚洲第一区二区三区不卡| 国产不卡一卡二| 我要看日韩黄色一级片| 久久久久性生活片| 亚洲中文日韩欧美视频| 午夜免费成人在线视频| 我要看日韩黄色一级片| 熟女人妻精品中文字幕| 天堂av国产一区二区熟女人妻| 午夜福利在线观看免费完整高清在 | 99久久精品国产国产毛片| 久久久久精品国产欧美久久久| 国产白丝娇喘喷水9色精品| 精品人妻视频免费看| 听说在线观看完整版免费高清| 国产精品一区www在线观看 | 日本a在线网址| 久久久国产成人精品二区| 午夜免费男女啪啪视频观看 | 国产精品三级大全| 午夜激情欧美在线| 中国美白少妇内射xxxbb| 日韩欧美国产一区二区入口| 日韩欧美免费精品| 亚洲avbb在线观看| 小说图片视频综合网站| 久久精品影院6| 国产亚洲精品久久久com| 国产在视频线在精品| 白带黄色成豆腐渣| 男女下面进入的视频免费午夜| 男女之事视频高清在线观看| 日本三级黄在线观看| 久久久久国产精品人妻aⅴ院| 日本与韩国留学比较| 成人特级av手机在线观看| 亚洲国产色片| 国产亚洲精品综合一区在线观看| 国产男靠女视频免费网站| 露出奶头的视频| 久久久久性生活片| 久久精品国产亚洲av涩爱 | 国产爱豆传媒在线观看| 丰满乱子伦码专区| 精品久久久久久成人av| 欧美区成人在线视频| 国产免费av片在线观看野外av| 毛片一级片免费看久久久久 | 欧美+亚洲+日韩+国产| 国产一区二区激情短视频| 乱系列少妇在线播放| 欧美最黄视频在线播放免费| 午夜影院日韩av| 欧美性感艳星| 亚洲久久久久久中文字幕| h日本视频在线播放| 99久久精品热视频| 99精品久久久久人妻精品| 国产探花在线观看一区二区| 国产av在哪里看| 亚洲内射少妇av| 国产真实乱freesex| 99热精品在线国产| 村上凉子中文字幕在线| 免费搜索国产男女视频| 中文字幕人妻熟人妻熟丝袜美| 白带黄色成豆腐渣| 制服丝袜大香蕉在线| 午夜激情欧美在线| 国产白丝娇喘喷水9色精品| 精品一区二区免费观看| 男人和女人高潮做爰伦理| 国内精品久久久久精免费| 精品一区二区三区人妻视频| 黄色视频,在线免费观看| 国产亚洲欧美98| 国产亚洲精品av在线| 日韩欧美精品v在线| 欧美成人性av电影在线观看| 日韩欧美国产在线观看| 亚洲 国产 在线| 男女边吃奶边做爰视频| 干丝袜人妻中文字幕| 天天一区二区日本电影三级| 97超视频在线观看视频| 51国产日韩欧美| 亚洲成av人片在线播放无| 国产精品福利在线免费观看| 毛片一级片免费看久久久久 | 午夜福利18| 亚洲av一区综合| 国产单亲对白刺激| 网址你懂的国产日韩在线| 亚洲美女视频黄频| 麻豆国产av国片精品| 国产精品野战在线观看| 97人妻精品一区二区三区麻豆| 一个人免费在线观看电影| 亚洲成a人片在线一区二区| 一区二区三区高清视频在线| 亚洲最大成人av| x7x7x7水蜜桃| 欧美日韩国产亚洲二区| 黄色女人牲交| 亚洲四区av| 亚洲av成人精品一区久久| 色哟哟哟哟哟哟| 欧美日韩中文字幕国产精品一区二区三区| 又爽又黄无遮挡网站| .国产精品久久| 亚洲av美国av| 又爽又黄无遮挡网站| 免费av不卡在线播放| 免费高清视频大片| 极品教师在线视频| av在线蜜桃| 露出奶头的视频| 成熟少妇高潮喷水视频| 国产免费一级a男人的天堂| 偷拍熟女少妇极品色| 老司机午夜福利在线观看视频| 看十八女毛片水多多多| 久久久久久大精品| 欧美国产日韩亚洲一区| av中文乱码字幕在线| 日韩av在线大香蕉| 免费一级毛片在线播放高清视频| 亚洲美女视频黄频| 久久久久性生活片| 精品一区二区三区视频在线观看免费| 国产伦人伦偷精品视频| 日韩大尺度精品在线看网址| 男女下面进入的视频免费午夜| 国内精品宾馆在线| 日本熟妇午夜| 琪琪午夜伦伦电影理论片6080| 免费看av在线观看网站| 精品一区二区三区视频在线观看免费| 亚洲七黄色美女视频| 中文在线观看免费www的网站| 亚洲成人久久性| 极品教师在线免费播放| 最近最新中文字幕大全电影3| 国内少妇人妻偷人精品xxx网站| 91在线精品国自产拍蜜月| 国产大屁股一区二区在线视频| 成人毛片a级毛片在线播放| 亚洲av第一区精品v没综合| 级片在线观看| 村上凉子中文字幕在线| 最近在线观看免费完整版| 国产亚洲精品av在线| 午夜影院日韩av| 性色avwww在线观看| 日韩欧美在线乱码| 国产精品一及| 成人二区视频| 99久久久亚洲精品蜜臀av| 嫩草影院新地址| 国产大屁股一区二区在线视频| 日日啪夜夜撸| 一本久久中文字幕| 久久久久国产精品人妻aⅴ院| 波多野结衣高清作品| 久久久国产成人精品二区| 亚洲成人免费电影在线观看| 制服丝袜大香蕉在线| 国产精品一区二区免费欧美| 91av网一区二区| 美女cb高潮喷水在线观看| 免费在线观看日本一区| 色噜噜av男人的天堂激情| av女优亚洲男人天堂| 日本熟妇午夜| 露出奶头的视频| 高清在线国产一区| 村上凉子中文字幕在线| 能在线免费观看的黄片| 欧美最新免费一区二区三区| 偷拍熟女少妇极品色| 精品人妻一区二区三区麻豆 | h日本视频在线播放| or卡值多少钱| 99久久成人亚洲精品观看| 久久欧美精品欧美久久欧美| 成年女人毛片免费观看观看9| 一个人看视频在线观看www免费| 免费av毛片视频| 在线免费观看不下载黄p国产 | 欧美xxxx黑人xx丫x性爽| 波多野结衣巨乳人妻| 最新中文字幕久久久久| 久久久久性生活片| 国产女主播在线喷水免费视频网站 | 亚洲av第一区精品v没综合| 黄色欧美视频在线观看| 国产精品,欧美在线| 女人被狂操c到高潮| 日本黄大片高清| 久久久久久久精品吃奶| 久久久久久久久中文| 国产精品一及| 亚洲久久久久久中文字幕| 日韩欧美 国产精品| 九色成人免费人妻av| 直男gayav资源| 日韩一本色道免费dvd| 亚洲国产精品成人综合色| 成人精品一区二区免费| 热99在线观看视频| 麻豆av噜噜一区二区三区| 12—13女人毛片做爰片一| 日本爱情动作片www.在线观看 | av在线蜜桃| 69av精品久久久久久| 免费在线观看日本一区| 久久中文看片网| 91精品国产九色| 此物有八面人人有两片| 18禁裸乳无遮挡免费网站照片| 国产亚洲精品av在线| 美女高潮喷水抽搐中文字幕| 亚洲男人的天堂狠狠| 男女啪啪激烈高潮av片| 嫁个100分男人电影在线观看| 黄色配什么色好看| 九九久久精品国产亚洲av麻豆| 97超级碰碰碰精品色视频在线观看| 国产老妇女一区| 观看美女的网站| 午夜精品一区二区三区免费看| 一a级毛片在线观看| 夜夜夜夜夜久久久久| 动漫黄色视频在线观看| 欧美xxxx黑人xx丫x性爽| 久久久久久久久久成人| or卡值多少钱| 真人做人爱边吃奶动态| 亚洲精品色激情综合| 亚洲不卡免费看| 亚洲欧美日韩无卡精品| 一进一出抽搐动态| 国产高清视频在线播放一区| 亚洲va在线va天堂va国产| 欧美最黄视频在线播放免费| 日韩精品有码人妻一区| 成年免费大片在线观看| 欧美日韩综合久久久久久 | 亚洲美女搞黄在线观看 | 欧美成人性av电影在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 国产精品精品国产色婷婷| 亚洲欧美日韩高清在线视频| 精品午夜福利在线看| 在线播放国产精品三级| 成人国产麻豆网| 免费看日本二区| 日本五十路高清| 精品一区二区三区视频在线| 三级国产精品欧美在线观看| 男人和女人高潮做爰伦理| 亚洲电影在线观看av| 亚洲av第一区精品v没综合| 日韩大尺度精品在线看网址| 国产精品免费一区二区三区在线| 亚洲国产精品久久男人天堂| 久久精品久久久久久噜噜老黄 | eeuss影院久久| 亚洲精品456在线播放app | 欧美高清性xxxxhd video| 日本色播在线视频| 久久天躁狠狠躁夜夜2o2o| 十八禁国产超污无遮挡网站| 动漫黄色视频在线观看| 啦啦啦韩国在线观看视频| 麻豆久久精品国产亚洲av| 日本爱情动作片www.在线观看 | 欧美性猛交╳xxx乱大交人| 日本一二三区视频观看| 中文字幕人妻熟人妻熟丝袜美| 亚洲不卡免费看| 成人av一区二区三区在线看| 一区二区三区免费毛片| 国产精品野战在线观看| 日韩精品有码人妻一区| 熟女人妻精品中文字幕| 少妇的逼水好多| 国产精品一区二区性色av| 久久天躁狠狠躁夜夜2o2o| 国产精品人妻久久久久久| 精品一区二区免费观看| 久久精品国产亚洲网站| 欧美激情国产日韩精品一区| 国产私拍福利视频在线观看| 精品无人区乱码1区二区| 亚洲国产精品合色在线| 国产私拍福利视频在线观看| 久久人人爽人人爽人人片va| 久99久视频精品免费| 久久国产精品人妻蜜桃| 欧美一区二区精品小视频在线| 亚洲成人久久性| 少妇的逼好多水| 麻豆精品久久久久久蜜桃| 看片在线看免费视频| 3wmmmm亚洲av在线观看| 俺也久久电影网| 免费在线观看日本一区| 一级毛片久久久久久久久女| 午夜精品一区二区三区免费看| 97碰自拍视频| 久久久久国内视频| 亚洲人成网站高清观看| 亚洲精品在线观看二区| 久久九九热精品免费| 老司机深夜福利视频在线观看| 国产黄a三级三级三级人| 日韩欧美在线二视频| 人人妻人人澡欧美一区二区| 男女之事视频高清在线观看| 国产精品野战在线观看| 丰满人妻一区二区三区视频av| 热99re8久久精品国产| 99视频精品全部免费 在线| 国产黄a三级三级三级人| 级片在线观看| 亚洲国产精品sss在线观看| 简卡轻食公司| 国产精品久久久久久久电影| 制服丝袜大香蕉在线| 天堂av国产一区二区熟女人妻| 日日撸夜夜添| 亚洲av电影不卡..在线观看| 亚洲av中文av极速乱 | 在线国产一区二区在线| 99久久精品热视频| 亚洲精品456在线播放app | 无人区码免费观看不卡| 一级黄色大片毛片| 一个人免费在线观看电影| 色综合婷婷激情| 天堂√8在线中文| 女的被弄到高潮叫床怎么办 | 国产中年淑女户外野战色| 亚洲经典国产精华液单| 我要搜黄色片| 国产精品久久久久久久久免| 精品不卡国产一区二区三区| 十八禁网站免费在线| 蜜桃亚洲精品一区二区三区| 一区二区三区高清视频在线| 18禁裸乳无遮挡免费网站照片| 国产免费男女视频| 变态另类丝袜制服| 久久久久精品国产欧美久久久| 五月伊人婷婷丁香| 麻豆av噜噜一区二区三区| 国产高潮美女av| 99久久精品国产国产毛片| 免费av不卡在线播放| 亚洲av熟女| 亚洲性久久影院| 狂野欧美激情性xxxx在线观看| 波多野结衣巨乳人妻| 成人国产麻豆网| 日本撒尿小便嘘嘘汇集6| 午夜福利高清视频| 嫩草影院入口| 国产精品女同一区二区软件 | 亚洲欧美清纯卡通| 欧美最黄视频在线播放免费| 黄色女人牲交| 干丝袜人妻中文字幕| 搡女人真爽免费视频火全软件 | 成人性生交大片免费视频hd| 亚洲色图av天堂| 日韩欧美精品免费久久| 日韩国内少妇激情av| 国模一区二区三区四区视频| 少妇的逼水好多| 成年女人永久免费观看视频| 人妻制服诱惑在线中文字幕| 一级毛片久久久久久久久女| 国产老妇女一区| 亚洲国产色片| 麻豆久久精品国产亚洲av| 97超级碰碰碰精品色视频在线观看| 免费在线观看影片大全网站| 国产日本99.免费观看| 在线天堂最新版资源| 国产成人aa在线观看| 日本与韩国留学比较| 99热这里只有是精品在线观看| 国内揄拍国产精品人妻在线| 蜜桃久久精品国产亚洲av| 国产男人的电影天堂91| 在线国产一区二区在线| 亚洲精品国产成人久久av| av专区在线播放| 国产成年人精品一区二区| netflix在线观看网站| 女生性感内裤真人,穿戴方法视频| 久久久久免费精品人妻一区二区| 听说在线观看完整版免费高清| 国产白丝娇喘喷水9色精品| 观看美女的网站| 亚洲人成网站在线播放欧美日韩| 性插视频无遮挡在线免费观看| 亚洲成人久久性| 最近在线观看免费完整版| 久久人人精品亚洲av| 乱系列少妇在线播放| 精品一区二区三区视频在线观看免费| 精品一区二区三区人妻视频| 波多野结衣巨乳人妻| 97人妻精品一区二区三区麻豆| 日韩高清综合在线| 内射极品少妇av片p| 国产精品免费一区二区三区在线| 精品久久久久久久久久免费视频| 亚洲欧美日韩无卡精品| 噜噜噜噜噜久久久久久91| 亚洲美女搞黄在线观看 | 亚洲人成伊人成综合网2020| 在线播放无遮挡| 国产伦精品一区二区三区视频9| 啦啦啦啦在线视频资源| 午夜福利18| 国产欧美日韩一区二区精品| 国产乱人伦免费视频| 欧美一区二区亚洲| 日韩在线高清观看一区二区三区 | 亚洲美女黄片视频| 桃色一区二区三区在线观看| 日韩精品中文字幕看吧| 国产精品久久视频播放| 日韩,欧美,国产一区二区三区 | a级毛片a级免费在线| 最新在线观看一区二区三区| 国产三级中文精品| 亚洲欧美清纯卡通| 国产精品一区二区性色av| 亚洲av一区综合| 联通29元200g的流量卡| 亚洲成人久久性| 日韩欧美精品免费久久| 亚洲国产高清在线一区二区三| 少妇裸体淫交视频免费看高清| 又爽又黄a免费视频| 一a级毛片在线观看| 99精品久久久久人妻精品| 人妻制服诱惑在线中文字幕| 国产精品久久久久久av不卡| 国产高清三级在线| www日本黄色视频网| 特大巨黑吊av在线直播| 亚洲av中文av极速乱 | 白带黄色成豆腐渣| 亚洲精品在线观看二区| 国产亚洲av嫩草精品影院| 1000部很黄的大片| videossex国产| 国产成人一区二区在线| 欧美激情在线99| 久久久久久久久久成人| 91久久精品电影网| 久久久久久国产a免费观看| 99热精品在线国产| 精品乱码久久久久久99久播| 精品午夜福利视频在线观看一区| 91av网一区二区| 18禁黄网站禁片午夜丰满| 国产亚洲精品av在线| 欧美性感艳星| 欧美极品一区二区三区四区| 日本-黄色视频高清免费观看| 毛片一级片免费看久久久久 | 久久精品国产亚洲av涩爱 | 搡女人真爽免费视频火全软件 | 国产精品1区2区在线观看.| 亚洲国产欧洲综合997久久,| 亚洲av不卡在线观看| 久久精品国产自在天天线| 人妻少妇偷人精品九色| 国产精品一区二区三区四区免费观看 | 亚洲人成伊人成综合网2020| 国语自产精品视频在线第100页| 国产精品日韩av在线免费观看| 国产精品不卡视频一区二区| 在线免费观看的www视频| 99精品在免费线老司机午夜| 观看美女的网站| 搡老熟女国产l中国老女人| 色吧在线观看| 精品久久久久久久久久久久久| 国产精品久久久久久精品电影| 欧美丝袜亚洲另类 | АⅤ资源中文在线天堂| a级一级毛片免费在线观看| 国产不卡一卡二| 中文在线观看免费www的网站| 99热只有精品国产| 99久久成人亚洲精品观看| 天堂√8在线中文| 国产69精品久久久久777片| 国产精品无大码| 日本 欧美在线| 亚洲国产高清在线一区二区三| 高清毛片免费观看视频网站| 中文字幕av成人在线电影| 久久久久九九精品影院| 日韩人妻高清精品专区| 高清毛片免费观看视频网站| 亚洲 国产 在线| 亚洲欧美清纯卡通| 男人舔奶头视频| 一卡2卡三卡四卡精品乱码亚洲| 永久网站在线| 毛片一级片免费看久久久久 | h日本视频在线播放| 久久午夜亚洲精品久久| 国产精品1区2区在线观看.| 亚洲四区av| 国产亚洲精品久久久久久毛片| 在线观看午夜福利视频| 欧美日韩国产亚洲二区| 欧美在线一区亚洲| 91狼人影院| 波多野结衣高清作品| 亚洲欧美精品综合久久99|