地閃回?fù)綦姶艌龇植继卣?/h1>
2020-08-13 07:05:56黃鈺崔蕾張弛唐辟如
科技創(chuàng)新與應(yīng)用 2020年24期
黃鈺 崔蕾 張弛 唐辟如
摘 ?要：建立了空間計算模型，利用MATLAB編程實現(xiàn)了時域有限差分法（FDTD）法對閃電回?fù)艨臻g電磁場的計算和分析。結(jié)果發(fā)現(xiàn)：隨著與閃電回?fù)敉ǖ谰嚯x的增加，垂直電場的幅值有明顯的減弱，波形也有明顯的變化，在較遠(yuǎn)處，波形出現(xiàn)明顯形態(tài)異常變化，隨著高度的增加，垂直電場幅值也變小，波形變緩，但在距離較近的時候表征更為明顯;水平電場幅值比同距離下垂直電場小很多，隨著高度、距離的增加，電場減弱顯著，且曲線形態(tài)趨于平緩;徑向磁場除了隨著距離的增加幅值在減小，不管距離遠(yuǎn)近，高度的高低，波形變化都不明顯。
關(guān)鍵詞：FDTD;地閃回?fù)?電場分布;波形變化
中圖分類號：P427.35 文獻(xiàn)標(biāo)志碼：A ? ? ? ? 文章編號：2095-2945（2020）24-0004-03
Abstract： The spatial calculation model is established， and the finite-difference time-domain （FDTD） method is used to calculate and analyze the space electromagnetic field of lightning return stroke by MATLAB programming. The results show that with the increase of the distance from the lightning return stroke channel， the amplitude of the vertical electric field weakens obviously， and the waveform also has obvious changes; at a distance， the waveform shows obvious abnormal changes; with the increase of height， the amplitude of the vertical electric field becomes smaller and the waveform slows down， but the characterization is more obvious when the distance is closer; the amplitude of the horizontal electric field is much smaller than that of the vertical electric field at the same distance， and the electric field weakens significantly with the increase of height and distance， and the shape of the curve tends to be smooth; except that the amplitude of the radial magnetic field decreases with the increase of distance， regardless of the distance and the height， the waveform change is not obvious.
Keywords： FDTD; ground lightning return stroke; electric field distribution; waveform change
引言
閃電是較為常見的一種強對流天氣現(xiàn)象，其電磁輻射也被稱作閃電電磁脈沖（LEMP）對微電子設(shè)備常常造成損害，研究者們對LEMP在架空電線或電纜的耦合效應(yīng)展開了很多研究[1-2]，并發(fā)現(xiàn)其實閃電回?fù)綦姶艌鲈诳臻g傳播時存在一定的規(guī)律，且受地物的影響[3-5]，這卻是提高閃電探測技術(shù)和分析的重要一環(huán)。傳統(tǒng)的有限差分法[6]是利用網(wǎng)格剖分將求解區(qū)域離散化為網(wǎng)格和節(jié)點，以差分原理為基礎(chǔ)，以各離散點上函數(shù)的差商來近似替代該點的偏導(dǎo)數(shù)，把要求解的邊值問題轉(zhuǎn)化為一組相應(yīng)的差分方程問題，面向高頻電磁場的傳輸、輻射、散射和透入等工程問題的需求，高頻電磁場的傳輸基于麥克斯韋方程組中旋度方程直接轉(zhuǎn)化為了差分方程的時域有限差分（FDTD）。雷電輻射本身是一個相當(dāng)復(fù)雜的電磁問題，F(xiàn)DTD由于是直接建立在電磁問題普遍實用的麥克斯韋方程組之上的，它解決問題的方法具有通用性，各個問題的不同之處僅僅在于其幾何結(jié)構(gòu)和電磁數(shù)的不同，而沒有解決方法上的差異，這樣對于分析電磁輻射問題就具有良好的適用性。本文基于FDTD方法對地閃回?fù)綦姶怒h(huán)境分布特征進(jìn)行分析。
1 建模與計算方法
為了使計算不過于復(fù)雜，假定（1）沿先導(dǎo)通道的電荷分布均勻， 雷擊垂直于大地。雷電通道與地面垂直，地面假設(shè)為水平面，不考慮地形和曲率影響;（2）雷電流沿雷電通道向上傳播到雷云后，消失在云中，不再反射返回地面;（3）大地為良導(dǎo)體，無損耗，即電導(dǎo)率無限大，計算空間模型如圖1所示。
電流模型采用傳輸線（TL）模型：該模型中主放電發(fā)生后，電流從先導(dǎo)通道的底部以一定的速度沿通道無衰減地向上傳播，通道為理想的傳輸線。對于一定的主放電速度，通道中任意點z′處電流的表達(dá)式為：
雷電回?fù)羲俣热=1.3×108m/s，雷電通道高度H為750米。本文采用IEEE推薦的8/20？滋s雙指函數(shù)波作為雷電通道的基電流：
由以上設(shè)定的計算空間模型、電流模型及給定相應(yīng)電磁問題的初始值，采用FDTD方法就可以逐步推進(jìn)地求得以后各個時刻空間電磁場的分布，具體計算基于MATLAB編程完成。
2 地閃回?fù)綦姶艌龇植继卣鞣治?/p>
假定周邊無建筑情況下，計算分析距離累計通道不同距離、或處于不同高度處的電磁場的分布特征。結(jié)合圖2發(fā)現(xiàn)，隨著與雷電通道的距離的增加，垂直電場的幅值有明顯的減弱，100m處最大幅值超10kV/m，而5000m處以下降至0.05kV/m;此外垂直電場波形也有明顯的變化，隨著距離增加整體上波形變緩，時間有個明顯的推移，在較遠(yuǎn)處，波形形態(tài)變化異常（圖2c）。在近距離距離100m時，電場隨著高度的增加有一定變化，隨著距離的增加這種變化越來越不明顯。
對100m處的垂直電場隨高度的變化進(jìn)行進(jìn)一步分析，如圖3所示，隨著高度的增加，垂直電場變化較為明顯，除幅值變小外波形也明顯變緩。對于水平電場，在同一距離情況下（d=100時）幅值比同距離下垂直電場小很多，且隨著高度的增加，電場減弱顯著（圖4a）;計算不同距離處的同一高度水平電場發(fā)現(xiàn)隨著距離的增加，電場值也在減弱，且曲線形態(tài)趨于平緩，變化較為明顯（圖4b）。
對徑向磁場進(jìn)行計算發(fā)現(xiàn)，磁場不管距離遠(yuǎn)近，高度的高低，波形變化都不明顯，除了隨著距離的增加幅值在減小，同一距離處隨著高度的增加磁場的大小也無明顯的變化。
3 結(jié)論
文中的電場的變化主要由衰減引起，但這與設(shè)定的雷擊通道高度也有密切聯(lián)系，且雷擊高處，如樹木，建筑物時，分布特征肯定會有所區(qū)別，且本研究過程中模型較為理想，尤其是閃電回?fù)敉ǖ缹嶋H過程中存在分形，在下一步的研究中將考慮更多的實際因素，獲取更為精確的結(jié)果。
參考文獻(xiàn)：
[1]王浩，陳亞洲，萬浩江.雷電電磁脈沖對雙絞線耦合的仿真研究[J].國外電子測量技術(shù)，2015，34（1）：40-44.
[2]郄秀書，張其林，楊靜.雷電在水平導(dǎo)體中產(chǎn)生感應(yīng)過電壓的觀測及數(shù)值模擬研究[J].物理學(xué)報，2008，57（3）：1968-1975.
[3]賈新民，嚴(yán)文.有限差分法求解拉普拉斯方程[J].昌吉學(xué)院學(xué)報，2009（5）：105-109.
[4]耿雪瑩，張其林，劉明遠(yuǎn).地面建筑物（群）對雷暴云大氣電場影響的模擬研究[J].氣象科技，2012，40（5）：827-833.
[5]郭秀峰，譚涌波，郭鳳霞，等.建筑物尖端對大氣電場畸變影響的數(shù)值計算[J].應(yīng)用氣象學(xué)報，2013，24（2）：189-196.
[6]王琳，陳亞洲，杜思尚，等.地表近區(qū)雷電回?fù)綦姶艌龅挠嬎闩c特征分析[J].高電壓技術(shù)，2012，38（11）：2899-2905.
猜你喜歡
雨天防雷電要選對雨傘中老年保健(2021年5期)2021-08-24 07:08:30 數(shù)列與差分新世紀(jì)智能(數(shù)學(xué)備考)(2021年5期)2021-07-28 06:19:46 外加正交電磁場等離子體中電磁波透射特性空間科學(xué)學(xué)報(2021年6期)2021-03-09 06:20:00 雷電小學(xué)生學(xué)習(xí)指導(dǎo)(低年級)(2019年5期)2019-04-29 09:11:00 任意方位電偶源的MCSEM電磁場三維正演石油地球物理勘探(2017年4期)2017-12-18 07:15:09 電磁場與電磁波課程教學(xué)改革探析湖南城市學(xué)院學(xué)報(自然科學(xué)版)(2016年2期)2016-12-01 04:07:07 計算機機房的雷電防護(hù)中國市場(2016年45期)2016-05-17 05:15:53 基于差分隱私的大數(shù)據(jù)隱私保護(hù)信息安全研究(2015年3期)2015-02-28 20:17:57 相對差分單項測距△DOR太空探索(2014年1期)2014-07-10 13:41:50 差分放大器在生理學(xué)中的應(yīng)用四川生理科學(xué)雜志(2014年2期)2014-02-28 14:09:20
黃鈺 崔蕾 張弛 唐辟如
摘 ?要：建立了空間計算模型，利用MATLAB編程實現(xiàn)了時域有限差分法（FDTD）法對閃電回?fù)艨臻g電磁場的計算和分析。結(jié)果發(fā)現(xiàn)：隨著與閃電回?fù)敉ǖ谰嚯x的增加，垂直電場的幅值有明顯的減弱，波形也有明顯的變化，在較遠(yuǎn)處，波形出現(xiàn)明顯形態(tài)異常變化，隨著高度的增加，垂直電場幅值也變小，波形變緩，但在距離較近的時候表征更為明顯;水平電場幅值比同距離下垂直電場小很多，隨著高度、距離的增加，電場減弱顯著，且曲線形態(tài)趨于平緩;徑向磁場除了隨著距離的增加幅值在減小，不管距離遠(yuǎn)近，高度的高低，波形變化都不明顯。
關(guān)鍵詞：FDTD;地閃回?fù)?電場分布;波形變化
中圖分類號：P427.35 文獻(xiàn)標(biāo)志碼：A ? ? ? ? 文章編號：2095-2945（2020）24-0004-03
Abstract： The spatial calculation model is established， and the finite-difference time-domain （FDTD） method is used to calculate and analyze the space electromagnetic field of lightning return stroke by MATLAB programming. The results show that with the increase of the distance from the lightning return stroke channel， the amplitude of the vertical electric field weakens obviously， and the waveform also has obvious changes; at a distance， the waveform shows obvious abnormal changes; with the increase of height， the amplitude of the vertical electric field becomes smaller and the waveform slows down， but the characterization is more obvious when the distance is closer; the amplitude of the horizontal electric field is much smaller than that of the vertical electric field at the same distance， and the electric field weakens significantly with the increase of height and distance， and the shape of the curve tends to be smooth; except that the amplitude of the radial magnetic field decreases with the increase of distance， regardless of the distance and the height， the waveform change is not obvious.
Keywords： FDTD; ground lightning return stroke; electric field distribution; waveform change
引言
閃電是較為常見的一種強對流天氣現(xiàn)象，其電磁輻射也被稱作閃電電磁脈沖（LEMP）對微電子設(shè)備常常造成損害，研究者們對LEMP在架空電線或電纜的耦合效應(yīng)展開了很多研究[1-2]，并發(fā)現(xiàn)其實閃電回?fù)綦姶艌鲈诳臻g傳播時存在一定的規(guī)律，且受地物的影響[3-5]，這卻是提高閃電探測技術(shù)和分析的重要一環(huán)。傳統(tǒng)的有限差分法[6]是利用網(wǎng)格剖分將求解區(qū)域離散化為網(wǎng)格和節(jié)點，以差分原理為基礎(chǔ)，以各離散點上函數(shù)的差商來近似替代該點的偏導(dǎo)數(shù)，把要求解的邊值問題轉(zhuǎn)化為一組相應(yīng)的差分方程問題，面向高頻電磁場的傳輸、輻射、散射和透入等工程問題的需求，高頻電磁場的傳輸基于麥克斯韋方程組中旋度方程直接轉(zhuǎn)化為了差分方程的時域有限差分（FDTD）。雷電輻射本身是一個相當(dāng)復(fù)雜的電磁問題，F(xiàn)DTD由于是直接建立在電磁問題普遍實用的麥克斯韋方程組之上的，它解決問題的方法具有通用性，各個問題的不同之處僅僅在于其幾何結(jié)構(gòu)和電磁數(shù)的不同，而沒有解決方法上的差異，這樣對于分析電磁輻射問題就具有良好的適用性。本文基于FDTD方法對地閃回?fù)綦姶怒h(huán)境分布特征進(jìn)行分析。
1 建模與計算方法
為了使計算不過于復(fù)雜，假定（1）沿先導(dǎo)通道的電荷分布均勻， 雷擊垂直于大地。雷電通道與地面垂直，地面假設(shè)為水平面，不考慮地形和曲率影響;（2）雷電流沿雷電通道向上傳播到雷云后，消失在云中，不再反射返回地面;（3）大地為良導(dǎo)體，無損耗，即電導(dǎo)率無限大，計算空間模型如圖1所示。
電流模型采用傳輸線（TL）模型：該模型中主放電發(fā)生后，電流從先導(dǎo)通道的底部以一定的速度沿通道無衰減地向上傳播，通道為理想的傳輸線。對于一定的主放電速度，通道中任意點z′處電流的表達(dá)式為：
雷電回?fù)羲俣热=1.3×108m/s，雷電通道高度H為750米。本文采用IEEE推薦的8/20？滋s雙指函數(shù)波作為雷電通道的基電流：
由以上設(shè)定的計算空間模型、電流模型及給定相應(yīng)電磁問題的初始值，采用FDTD方法就可以逐步推進(jìn)地求得以后各個時刻空間電磁場的分布，具體計算基于MATLAB編程完成。
2 地閃回?fù)綦姶艌龇植继卣鞣治?/p>
假定周邊無建筑情況下，計算分析距離累計通道不同距離、或處于不同高度處的電磁場的分布特征。結(jié)合圖2發(fā)現(xiàn)，隨著與雷電通道的距離的增加，垂直電場的幅值有明顯的減弱，100m處最大幅值超10kV/m，而5000m處以下降至0.05kV/m;此外垂直電場波形也有明顯的變化，隨著距離增加整體上波形變緩，時間有個明顯的推移，在較遠(yuǎn)處，波形形態(tài)變化異常（圖2c）。在近距離距離100m時，電場隨著高度的增加有一定變化，隨著距離的增加這種變化越來越不明顯。
對100m處的垂直電場隨高度的變化進(jìn)行進(jìn)一步分析，如圖3所示，隨著高度的增加，垂直電場變化較為明顯，除幅值變小外波形也明顯變緩。對于水平電場，在同一距離情況下（d=100時）幅值比同距離下垂直電場小很多，且隨著高度的增加，電場減弱顯著（圖4a）;計算不同距離處的同一高度水平電場發(fā)現(xiàn)隨著距離的增加，電場值也在減弱，且曲線形態(tài)趨于平緩，變化較為明顯（圖4b）。
對徑向磁場進(jìn)行計算發(fā)現(xiàn)，磁場不管距離遠(yuǎn)近，高度的高低，波形變化都不明顯，除了隨著距離的增加幅值在減小，同一距離處隨著高度的增加磁場的大小也無明顯的變化。
3 結(jié)論
文中的電場的變化主要由衰減引起，但這與設(shè)定的雷擊通道高度也有密切聯(lián)系，且雷擊高處，如樹木，建筑物時，分布特征肯定會有所區(qū)別，且本研究過程中模型較為理想，尤其是閃電回?fù)敉ǖ缹嶋H過程中存在分形，在下一步的研究中將考慮更多的實際因素，獲取更為精確的結(jié)果。
參考文獻(xiàn)：
[1]王浩，陳亞洲，萬浩江.雷電電磁脈沖對雙絞線耦合的仿真研究[J].國外電子測量技術(shù)，2015，34（1）：40-44.
[2]郄秀書，張其林，楊靜.雷電在水平導(dǎo)體中產(chǎn)生感應(yīng)過電壓的觀測及數(shù)值模擬研究[J].物理學(xué)報，2008，57（3）：1968-1975.
[3]賈新民，嚴(yán)文.有限差分法求解拉普拉斯方程[J].昌吉學(xué)院學(xué)報，2009（5）：105-109.
[4]耿雪瑩，張其林，劉明遠(yuǎn).地面建筑物（群）對雷暴云大氣電場影響的模擬研究[J].氣象科技，2012，40（5）：827-833.
[5]郭秀峰，譚涌波，郭鳳霞，等.建筑物尖端對大氣電場畸變影響的數(shù)值計算[J].應(yīng)用氣象學(xué)報，2013，24（2）：189-196.
[6]王琳，陳亞洲，杜思尚，等.地表近區(qū)雷電回?fù)綦姶艌龅挠嬎闩c特征分析[J].高電壓技術(shù)，2012，38（11）：2899-2905.
猜你喜歡
中老年保健(2021年5期)2021-08-24 07:08:30
新世紀(jì)智能(數(shù)學(xué)備考)(2021年5期)2021-07-28 06:19:46
空間科學(xué)學(xué)報(2021年6期)2021-03-09 06:20:00
小學(xué)生學(xué)習(xí)指導(dǎo)(低年級)(2019年5期)2019-04-29 09:11:00
石油地球物理勘探(2017年4期)2017-12-18 07:15:09
湖南城市學(xué)院學(xué)報(自然科學(xué)版)(2016年2期)2016-12-01 04:07:07
中國市場(2016年45期)2016-05-17 05:15:53
信息安全研究(2015年3期)2015-02-28 20:17:57
太空探索(2014年1期)2014-07-10 13:41:50
四川生理科學(xué)雜志(2014年2期)2014-02-28 14:09:20
