專利名稱：太陽活動經(jīng)度帶預(yù)報(bào)方法
技術(shù)領(lǐng)域：
本發(fā)明涉及太陽物理和空間天氣預(yù)報(bào)領(lǐng)域，尤其涉及太陽活動預(yù)報(bào)。
背景技術(shù)：
太陽活動預(yù)報(bào)是太陽物理和空間天氣預(yù)報(bào)的主要研究領(lǐng)域之一，研究歷史 悠久。太陽活動預(yù)報(bào)根據(jù)預(yù)報(bào)的提前時間長度可分為長期、中期和短期預(yù)報(bào)等 三類。在相鄰的兩類之間沒有公認(rèn)的嚴(yán)格界限。 一般認(rèn)為，長期預(yù)報(bào)是指提前 量為一年或幾年以上至幾十年或更長時間的預(yù)報(bào)，中期預(yù)報(bào)是指提前量為半個 或一個太陽自轉(zhuǎn)周以上至幾個月的預(yù)報(bào)，短期預(yù)報(bào)是指提前量為幾十分鐘至幾 天的預(yù)報(bào)。太陽活動預(yù)報(bào)的這三種分類主要依據(jù)是太陽上各種活動過程的時間 尺度，預(yù)報(bào)技術(shù)方法的差別以及預(yù)報(bào)使用者的要求。長期預(yù)報(bào)的預(yù)報(bào)內(nèi)容以太 陽黑子相對數(shù)極大值及其時間為主，中期預(yù)報(bào)主要包括太陽耀斑、質(zhì)子事件、
黑子相對數(shù)、太陽10cm射電流量、太陽總體活動水平等多項(xiàng)內(nèi)容，短期預(yù)報(bào) 內(nèi)容主要包括太陽耀斑、質(zhì)子事件、太陽10cm射電流量等。太陽活動對無線 電通訊、宇宙空間環(huán)境、氣象、水文、地震等都有影響。因此，預(yù)報(bào)太陽活動 對于宇宙航行、空間研究、國防、國民經(jīng)濟(jì)等各個方面都具有極其重要的意 義。 太陽活動經(jīng)度帶是指太陽活動發(fā)生相對頻繁的經(jīng)度區(qū)域。太陽活動經(jīng)度帶 研究始于20世紀(jì)60年代。研究發(fā)現(xiàn)了一些太陽活動經(jīng)度帶的存在現(xiàn)象。由于 太陽表面較差自轉(zhuǎn)對太陽活動經(jīng)度帶的影響，使得這些經(jīng)度帶在不同時間出現(xiàn) 的位置不同。因此，雖然早在19世紀(jì)末人們就發(fā)現(xiàn)了太陽活動經(jīng)度帶的存在證 據(jù)，但是這些經(jīng)度帶的具體位置不是完全一致。早在20世紀(jì)70年代北京天文 臺在對太陽活動進(jìn)行預(yù)報(bào)時就開始考慮了太陽活動經(jīng)度帶的一些特征。但是由 于觀測資料的不連續(xù)，此時考慮的只是太陽活動經(jīng)度帶的一些短期特征，即太陽活動經(jīng)度帶在未來一到三天(最多半月)時間內(nèi)的運(yùn)動趨勢。并且此時考慮 太陽活動經(jīng)度帶的運(yùn)動趨勢也是為其他預(yù)報(bào)內(nèi)容(如耀斑、質(zhì)子事件等)服務(wù) 的，并沒有建立一種可操作的太陽活動經(jīng)度帶預(yù)報(bào)模型，只是預(yù)報(bào)員根據(jù)太陽 活動經(jīng)度帶的以往記錄來判斷其運(yùn)動趨勢，很少甚至沒有數(shù)據(jù)計(jì)算，因此有很 大的主觀性，不同預(yù)報(bào)員的判斷結(jié)果會有所不同。
近幾年人們發(fā)現(xiàn)如果消除了太陽表面較差自轉(zhuǎn)的影響，太陽黑子出現(xiàn)的
位置長期U20年)穩(wěn)定，在南北半球上均存在兩個相距約180度的太陽黑子
經(jīng)度帶；太陽x射線耀斑的經(jīng)度爆發(fā)位置也長期(最近三個太陽活動周)穩(wěn)定，
在南北半球上也同樣存在兩個相距約180度的耀斑經(jīng)度帶。因此，利用這些規(guī)律 人們能夠建立一種方法，對太陽活動經(jīng)度帶進(jìn)行預(yù)報(bào)。

發(fā)明內(nèi)容
本發(fā)明的主要目的在于根據(jù)太陽活動長期統(tǒng)計(jì)規(guī)律建立一種可搡作的太陽 活動經(jīng)度帶預(yù)報(bào)方法，該方法不僅豐富了太陽活動中長期預(yù)報(bào)的內(nèi)容，也為短 期預(yù)報(bào)提供一定的參考價值。
為達(dá)到上述目的，本發(fā)明的技術(shù)方案是這樣實(shí)現(xiàn)的
A、 選取觀測數(shù)據(jù)。首選太陽X射線耀斑數(shù)據(jù)；在太陽活動極小期間，沒 有太陽耀斑時就用太陽黑子數(shù)據(jù)；
B、 根據(jù)步驟A選取的觀測數(shù)據(jù)計(jì)算出太陽自轉(zhuǎn)的相關(guān)參數(shù)和太陽活動經(jīng) 度帶的起始位置(卡林頓坐標(biāo))，并在此基礎(chǔ)上計(jì)算出太陽活動經(jīng)度帶的當(dāng)前位 置(卡林頓坐標(biāo))。由于太陽南北半球的轉(zhuǎn)動情況不完全相同，所以在確定太陽 自轉(zhuǎn)的相關(guān)參數(shù)和經(jīng)度帶初始位置時對南北半球分別操作；
C、 根據(jù)步驟B計(jì)算出的太陽自轉(zhuǎn)的相關(guān)參數(shù)、太陽活動經(jīng)度帶的起始位 置和當(dāng)前位置計(jì)算出(預(yù)報(bào))一定時間內(nèi)的太陽活動經(jīng)度帶中心位置(卡林頓坐
標(biāo))；
D、 將步驟C預(yù)報(bào)的太陽活動經(jīng)度帶中心位置的卡林頓坐標(biāo)換算成日面坐 標(biāo)，預(yù)報(bào)太陽活動經(jīng)度帶的日面位置。說明本預(yù)報(bào)方法所用到的時間是世界時(UT)。
下面分別介紹步驟B中確定太陽自轉(zhuǎn)參數(shù)和太陽活動經(jīng)度帶初始位置的原
理和具體實(shí)現(xiàn)方法。
步驟B的實(shí)現(xiàn)原理太陽表面較差自轉(zhuǎn)規(guī)律可以簡單表示為Q-Q。-5sin、 (①)，其中，Q。(度/天)為太陽赤道自轉(zhuǎn)角速度，￡ (度/天)為緯向較差修 正常數(shù)，這兩個參數(shù)是太陽自轉(zhuǎn)參數(shù)，p (度)是緯度，Q (度/天)是太陽上 緯度^處的轉(zhuǎn)動角速度。
根據(jù)耀斑級別加權(quán)法計(jì)算出第f個卡林頓自轉(zhuǎn)周(會合周期，27.2753天) 內(nèi)的太陽耀斑平均纟韋度仍(Z表示卡林頓自轉(zhuǎn)周序號)，則該自轉(zhuǎn)周內(nèi)太陽上a處 的轉(zhuǎn)動角速度n,可以表示為Q,-Q。-Bsin>,(②)。
設(shè)兩個經(jīng)度帶的初始位置為AM和A。p根據(jù)兩經(jīng)度帶相距180度可得 |八01-A。2| = 180° (③)。當(dāng)時間到第!'個自轉(zhuǎn)周的第*天時(A可以是小數(shù))'經(jīng)度
帶相對于卡林頓坐標(biāo)系的經(jīng)度漂移為^^ ;li(Q,-A)+WQi-^)(④)，其
中，7;為卡林頓自轉(zhuǎn)周周期(27.2753天)，iV。和TV,—,分別代表起始周和第(/-1)周， ^為卡林頓坐標(biāo)系的轉(zhuǎn)動角速度(恒星速度，14.1844度沃)。此時太陽活動經(jīng) 度帶的卡林頓坐標(biāo)為Afl =(入。1+^01110(1360°和八。=(A。2+M)mod360° ( )，其中， mod360。表示劃入0°-360°范圍內(nèi)。
第/周觀測到的第/個耀斑的經(jīng)度位置&和經(jīng)度帶中心位置的偏離記為 ~=min(k-A》r^-A,,|)(⑥)。耀斑位置和經(jīng)度帶中心位置的平均偏離為
e-丄SZ 2 (⑦)，其中，AT是耀斑總數(shù)。綜合方程②③④⑤⑥⑦得到的f的
表達(dá)式為
rnin(;^. —[A01 + ; 2 (Q。 —Qc —5sin2 p,.) + A:(Q。 —Qc —5sin2 ^)〗mod360'
一 y
o" & — [A01 + ;藝(Q0 — Qc — 5 sin2 ^) + A:(Q。 — Qc — 5 sin2仍)± 180° ] mod 360。 }其中，'土，當(dāng)A。^18(T時取'+，，反之取'-，，該表達(dá)式中只有Q。， B和A。,是待定 參數(shù)，A.為觀測量，7;和^為已知常量，^是由觀測計(jì)算得到的量?？梢栽谝?定范圍內(nèi)搜索出滿足e為最小的一組參數(shù)Q。、 B和A。" 步驟B的實(shí)現(xiàn)方法如下
Bl、計(jì)算數(shù)據(jù)起始點(diǎn)的自轉(zhuǎn)周序號iV。；從觀測數(shù)據(jù)中提取每個耀斑的經(jīng)
度，將第z自轉(zhuǎn)周發(fā)生的第/個耀斑的經(jīng)度記為;i"如果不是卡林頓經(jīng)度要將其
換算成卡林頓經(jīng)度；將該耀斑爆發(fā)的時間單位換算成以天為單位，記為"該 耀斑對應(yīng)的卡林頓自轉(zhuǎn)周序號記為iv,;以x射線耀斑級別為權(quán)重計(jì)算每個自轉(zhuǎn) 周的太陽活動平均煒度^，如果觀測數(shù)據(jù)是太陽黑子，則以黑子面積為權(quán)重計(jì) 算平均緯度仍。
X射線耀斑按X射線峰值流量從低到高分為A、 B、 C、 M和X五大級別。 由于A和B級耀斑級別太低，它們的X射線峰值流量和背景X射線流量區(qū)別 不大，因此忽略不計(jì)，即此處的耀斑觀測數(shù)據(jù)選擇C級及其以上的所有耀斑。 相鄰兩級別耀斑的峰值流量相差一個量級，耀斑M(jìn)l.O的峰值流量人，是耀斑
ci.o峰值流量y;的io倍，即厶-io義，如果以耀斑ci.o的強(qiáng)度計(jì)為單位i,則
耀斑M(jìn)1.0的強(qiáng)度為10，耀斑X1.0的強(qiáng)度為100,在一個自轉(zhuǎn)周內(nèi)如果共發(fā)生 了一個C1.0級、2個M1.0級和1個X1.0級耀斑，它們的緯度分別為10、 8、 12、和10,則該自轉(zhuǎn)周發(fā)生的耀斑總強(qiáng)度為y^l.0 + 1.0xl0x2 + 1.0xl00-121.0,每 個耀斑的權(quán)重為1.0/121.0、 10.0/121.0、 10.0/121.0和100.0/121.0,平均纟爭度
a=10x_H_+8x!+12xi^+10x!M=10.0 (度)。黑子面積的權(quán)重計(jì)算方法
'121.0 121.0 121.0 121.0
簡單，黑子觀測數(shù)據(jù)給出了黑子面積數(shù)值，用一個黑子的面積除以一個自轉(zhuǎn)周 內(nèi)所有黑子的總面積即是該黑子的面積權(quán)重，每個黑子面積權(quán)重乘以黑子緯度， 再相加就得到了黑子平均緯度。
B2、把步驟B1得到的數(shù)據(jù)&, yfc， ^， AA。和7V,，以及已知常量7;和a.,
代入表達(dá)式⑧。B3、在給定區(qū)間內(nèi)以給定步長搜索滿足f為最小的一組參數(shù)Q。， S和A。p
并相應(yīng)記為Q一， ^和A，。其中，Q。的搜索區(qū)間為[13.50, 15.00],步長為 0.01; B的搜索區(qū)間為
,步長為0.01; A。,的搜索區(qū)間為
， 步長為1°。
B4、根據(jù)步驟B3中得到的參數(shù) ,，~和八。1/;,，以及方程②④⑤計(jì)算出 太陽活動經(jīng)度帶的當(dāng)前位置A,.,。 步驟C的具體執(zhí)行方法如下
Cl、預(yù)報(bào)太陽活動經(jīng)度帶的轉(zhuǎn)動角速度Q(,+^。設(shè)當(dāng)前卡林頓自轉(zhuǎn)周為第/周，
計(jì)算當(dāng)前周的太陽活動經(jīng)度帶的轉(zhuǎn)動角速度Q,Q。力,—A,sinV,，其中，r已在
步驟B1中計(jì)算出。如果考慮太陽活動經(jīng)度帶在一個自轉(zhuǎn)周的角速度的變化可 以忽略，則可以將下一周(第(/ + 1)周)的經(jīng)度帶轉(zhuǎn)動速度0(,+1),取作經(jīng)度帶的當(dāng)
前速度Q,.，即q,+^-Q,，其中戶代表預(yù)報(bào)。
C2、預(yù)報(bào)太陽活動經(jīng)度帶中心位置的卡林頓坐標(biāo)。將第i+l周的起始時間定 為該周起始日的00:00時刻；根據(jù)公式Ane rf ={An-Qc)}mod360°計(jì)算出第i 周最后一天的太陽活動經(jīng)度帶位置的卡林頓坐標(biāo)A,^，其中yt代表第i周最后一 個耀斑的爆發(fā)時間到第i+l周起始日的00:00時刻經(jīng)過的時間，以天為單位；
根據(jù)公式A(叫一 ={A,w+A:(^+1)p-Qe)}mod36(t ，計(jì)算第i+l周 一個太陽活動經(jīng)度 帶中心位置的卡林頓坐標(biāo)，A:e
，當(dāng)A:取整數(shù)O, 1, 2…時對應(yīng)的是第 i+l周的第l， 2, 3，....曰00:00時刻太陽活動經(jīng)度帶的卡林頓坐標(biāo)。根據(jù)公式 A(,，2/;十18(T)mod36(r ,得到第i+l周另 一個太陽活動經(jīng)度帶中心位置的
卡林頓坐標(biāo)。
步驟D的具體執(zhí)行方法如下
Dl、將預(yù)報(bào)的兩個太陽活動經(jīng)度帶的卡林頓坐標(biāo)換算成日面坐標(biāo)。 本發(fā)明中的日面坐標(biāo)的緯度定義與通常所用的日面坐標(biāo)相同，不再敘述，而本發(fā)明中的日面坐標(biāo)經(jīng)度是這樣定義的日輪中心經(jīng)度為0°，從中心沿0緯
度向東到太陽背面的中心為0°~-180°，從中心沿0緯度向西至太陽背面的中心 為0°~180°,即本發(fā)明中的曰面經(jīng)度范圍為[-180°，180°],負(fù)為日面中心東側(cè)，正 為曰面中心西側(cè)。
計(jì)算出每天世界時00:00時的太陽活動經(jīng)度帶中心位置的曰面經(jīng)度八,41和 八,,2,其中，角標(biāo)s代表日面坐標(biāo)，6代表第A天。具體計(jì)算方法如下
Dl.l 根據(jù)步驟C2 的公式A(w)一HA,w+A:(Q(,+w-Qjmod36(r和 A(,+,)^^A(,+D^+18(T)mod36(r計(jì)算出第i+l周每曰00:00時太陽活動經(jīng)度帶的卡 林頓坐標(biāo)A^和A^;
D1.2計(jì)算出第i+l周每日00:00時曰面中心的卡林頓坐標(biāo)A雙；
D1.3計(jì)算出第i+l周每曰00:00時太陽活動經(jīng)度帶的曰面坐標(biāo)為 A," = A出-Aa和A,42 = Aci2 - Aa ,并將其劃入日面經(jīng)度范圍[-18(T,18(T]內(nèi)。
D2、確定太陽活動經(jīng)度帶在一天內(nèi)的曰面經(jīng)度變化范圍由于卡林頓坐標(biāo) 系是時刻轉(zhuǎn)動的，每天轉(zhuǎn)過的角度約14.1844度，而日面坐標(biāo)基本是固定的， 所以卡林頓經(jīng)度坐標(biāo)相對于日面經(jīng)度坐標(biāo)每天有14度左右的漂移，又由于地球 繞太陽公轉(zhuǎn)每天約l度，最后的視覺效果就是卡林頓坐標(biāo)系相對于日面坐標(biāo)系 每天轉(zhuǎn)動約13度，此范圍與由步驟D1得到的連續(xù)兩個世界時對應(yīng)的曰面經(jīng)度
和A,(一之間的變化范圍|A,41 - A#+1)11基本相同，預(yù)報(bào)的太陽活動經(jīng)度帶中心
位置的每天日面經(jīng)度變化范圍為、~A 和A,)2 。
D3、確定太陽活動經(jīng)度帶日面煒度變化范圍由于太陽活動出現(xiàn)的緯度范 圍也是比較集中的，并不是全日面的，所以本預(yù)報(bào)方法中的緯度范圍取為該自 轉(zhuǎn)周的平均緯度仍的上下10度，即^-1(r ^+l(T,但不能變符號，也就是說當(dāng)
—,|，時,煒度范圍相應(yīng)為0。 p,+l(T ( ^e
),
伊,的正負(fù)表示太陽的北半球和南半球。
D4、預(yù)報(bào)太陽活動經(jīng)度帶的日面位置。由步驟D2和D3得出北半球兩個經(jīng)度帶中心位置的每天曰面變化范圍為(A^ A^部，1(t 仍+1(t )和 (ai42~ai(i+1)2, a—10° a+10° ) ( a>1(T );或(a^ a卓部，(T 仍+l(r )和 (A,42~A,,)2，(T w+l(r )( 0、仍"0° )。同理可以得出南半球的兩個經(jīng)度帶中心 位置的每天變化范圍,即(A:、w A',闊,仍-1(Ti,+1(T )和(A',4廣A:拜, 10° ~ a+10° ) ( ^<-10° );或(A:、w ~ A,(4+1)1 ，仍—10° ~ 0° ) 和 (A',a廣A,(一2， a-10° 0° ) ( -10、仍《(t )。這樣，每天預(yù)報(bào)的全曰面360度范 圍內(nèi)共有四個活動區(qū)域。
說明1、每過一個自轉(zhuǎn)周，加入新太陽活動觀測數(shù)據(jù)重新搜索合適的一組 參數(shù)Q。， S和A。,，再對下一個自轉(zhuǎn)周進(jìn)行預(yù)報(bào)。在當(dāng)前自轉(zhuǎn)周如果有耀斑數(shù)據(jù)
就選耀斑數(shù)據(jù)，如果沒有耀斑數(shù)據(jù)就選黑子數(shù)據(jù)，如果連黑子數(shù)據(jù)也沒有，則 不用重新計(jì)算參數(shù)Q。， S和A。,，即步驟a和b省略。此時，用一年內(nèi)的太陽 活動平均緯度代替該自轉(zhuǎn)周的平均緯度a,只作步驟c和d即可。2、本預(yù)報(bào) 方法也可以對未來幾個月甚至一年的太陽活動經(jīng)度帶進(jìn)行預(yù)報(bào)，因?yàn)樘栕赞D(zhuǎn) 參數(shù)在一年內(nèi)的變化幅度不大，但需要估測好在預(yù)報(bào)時間段內(nèi)太陽活動發(fā)生的 平均緯度。
本發(fā)明所提供的太陽活動經(jīng)度帶的預(yù)報(bào)方法，具有以下優(yōu)點(diǎn)
1) 可搡作性強(qiáng)，具有客觀性，預(yù)報(bào)結(jié)果不會因預(yù)報(bào)員而不同；
2) 預(yù)報(bào)時間即可以是短期，也可以是中長期。 本發(fā)明在太陽活動經(jīng)度帶預(yù)報(bào)方面屬于開拓性的發(fā)明。


圖1為太陽活動經(jīng)度帶預(yù)報(bào)方法的流程示意圖。
具體實(shí)施例方式
本發(fā)明中所用的數(shù)據(jù)主要是由美國宇航局環(huán)境應(yīng)用同步衛(wèi)星 (Geostationary Operational Environmental Satellites,即GOES)觀測的太陽X射線耀斑數(shù)據(jù)。其次是由美國宇航局空間天氣預(yù)報(bào)中心(Space Weather Prediction Center,即SWPC)上公布的太陽事件(包括X射線耀斑)及活動區(qū)數(shù)據(jù)。在太 陽活動極小期間，有時連續(xù)幾個月也觀測不到太陽耀斑，此時就用太陽黑子數(shù) 據(jù)。從太陽耀斑觀測數(shù)據(jù)中提取耀斑級別，爆發(fā)時間和位置。從黑子觀測數(shù)據(jù) 中提取黑子的面積，出現(xiàn)時間和位置。
本發(fā)明所使用的程序tim2carr.pro來源于SOHO (Solar and Heliospheric Observatory,歐洲空間局和美國宇航局共建的太陽和太陽大氣觀測衛(wèi)星)網(wǎng)站 (http:〃sohowww.nascom.nasa.gov/solarsoft/gen/idl/solar/)上。
以2005年8月的一個太陽自轉(zhuǎn)周(第2033周)為當(dāng)前自轉(zhuǎn)周，對下一個 自轉(zhuǎn)周(第2034周)作太陽活動經(jīng)度帶預(yù)報(bào)，舉例說明本發(fā)明的具體實(shí)現(xiàn)方法， 圖l是本發(fā)明的該具體實(shí)施方式
的流程圖，由步驟101~106組成。
GOES數(shù)據(jù)網(wǎng)址http:〃www.ngdc.noaa.gov/stp/SOLAR/ftpsolarflares.html, 每年一個文件，數(shù)據(jù)詳細(xì)說明見該網(wǎng)站的format文件，原始數(shù)據(jù)中按耀斑發(fā)生 時間為序，南北半球的記錄都在一個文件中，此處是將南北半球分開后的，只 取北半球的列在此，原始數(shù)據(jù)中包括A和B級耀斑，此處只選擇C級以上耀斑， 其它任何符號沒有改動，因此，可將初步篩選的數(shù)據(jù)視為原始數(shù)據(jù)。
數(shù)據(jù)選擇從GOES觀測的第一個太陽活動周的第 一個C級耀斑直到第2033 自轉(zhuǎn)周的全部C級以上耀斑。GOES從1975年開始觀測，由于1975和1976 兩年為太陽活動低年，幾乎沒有C級耀斑，所以數(shù)據(jù)選擇從1977年開始。在 此給出了 1977年開始有記錄的兩個自轉(zhuǎn)周和2005年8月的記錄，中間數(shù)據(jù)略。 下邊是GOES的原始記錄。
31777770105 1457 1510 1508 N30W901N C 0 GOES 766 761229.6
31777770131 1726 1735 1730 N27E601B M 2 GOES 779 7702 5.4
31777770202 0215 0218 0216 N26E39SF C 0 GOES 779 7702 5.131777770202 0322 0326 0323 N26E39SF 7702 5.2
31777770202 0425 0615 0505 N28E361N 7702 5.0
31777770209 1458 1521 1505 N12E42SN 770212.8
31777770217 1339 1418 1350 N19W54SF 770213,4
31777770218 2344 2358 2350 N23W77SN 770213.0
C 0 GOES
C
GOES
C 0 GOES
C 0 GOES
C 1 GOES
779
779
784
787
787
31777050823 2027 2036 2031 N17W06SX C 19
1.1E+03 3.0E+03 31777050823 2.7E+03 1.2E+04
31777050823 1.6E+03 4.1E+03
31777050824 2.9E+03 l.OE+04
31777050824 1.7E+03 3.8E+03
31777050825 1.2E+04 1.2E+05
31777050826 2.3E+03 7.6E+03 31777050826 2.1E+03 9.9E+03 31777050826 2152
2158 2211 2205 N12E85SX
2310 2319 2315 N13E85SX
0155 0207 0201 N12E85SX
0407 0418 0411 N12E84SX
0431 0445 0440 N09E80SX
0241 0338 0312 N12E68SX
1139 1210 1155 N10E61SX
C 34
C 12
C 22
C 12
M 64
C 13
C 21
2208 2202 N10E50SX
14
C 23
GOES
GOES
GOES
GOES
GOES
GOES
GOES
GOES
GOES2.9E+03 9.8E+03
31777050827 2036 1.7E+03 6.5E+03
31777050828 1017 5.9E+03 5.3E+04 31777050831 1026 4.8E+03 UE+04
2105 2051 N10E37SX 1037 1028 N09E36SX 1251 1151 N13W13SX
C 14 GOES M 16 GOES C 20 GOES
首先以第一行數(shù)據(jù)為例，對數(shù)據(jù)進(jìn)行簡單說明
31777770105 1457 1510 1508 N30W901N C 0 GOES 766 761229.6
317775個字符表示GOES觀測站代碼 770105 6個字符表示年月日，1977年1月5曰 14574個字符表示耀斑開始時刻14:57 ) l510 4個字符表示耀斑結(jié)束時刻(15:10) 1508 4個字符表示耀斑最大時刻(15:08) N l個字符表示北半球，或S南半球
30 2個字符表示緯度
W 1個字符表示曰面中心西側(cè)，或E曰面東側(cè)
90 2個字符表示中經(jīng)距，即耀斑到日面中心的經(jīng)度間隔
1N 4個字符表示兩個字符，耀斑面積和亮度，不考慮此因素
C l個字符表示耀斑級別，或M、 X
0 兩個字符表示X射線密度，從1.0至U9.9，省略了小數(shù)點(diǎn)，O表示在O.O
到0.1之間，我們在計(jì)算時將C0級，即CO.O級耀斑的峰值流量計(jì)為 0.1,而不計(jì)為0.0,因此如果耀斑峰值流量記錄值為f,則計(jì)算時其峰 值流量為(f+l)/10.0,即C5為C0.6級，C23為C2.4級。
GOES 4個字符表示觀測衛(wèi)星GOES符號之后的數(shù)據(jù)，本方法不用這幾列數(shù)據(jù)，不再解釋。 1、提取數(shù)據(jù)(對應(yīng)步驟101)。
l.l從原始記錄中提取有用的幾列數(shù)據(jù)(耀斑爆發(fā)時間，位置和級別)，并將耀斑 級別轉(zhuǎn)換為可計(jì)算的形式，數(shù)據(jù)記錄形式變?yōu)?br> 7701051457N30W90CO.l
7701311726N27E60M0.3
7702020215N26E39CO.l
7702020322N26E39CO.l
7702020425N28E36C0.6
7702091458N12E42CO.l
7702171339N19W54CO.l
7702182344N23W77C0.2
0508232027N17W06C2.0
0508232158N12E85C3.5
0508232310N13E85C 1.3
0508240155N12E85C2.3
0508240407N12E84C 1.3
0508250431磨E80M6.5
0508260241N12E68C 1.3
0508261139N10E61C2.2
0508262152N10E50C2.4
0508272036廳E37C 1.5
0508281017N0犯36M 1.7
0508311026N13W13C2.12、計(jì)算每個耀斑爆發(fā)時間對應(yīng)的自轉(zhuǎn)周序號《、爆發(fā)位置的卡林頓經(jīng)度&和每 個自轉(zhuǎn)周的太陽活動平均綿度p,等方程⑧要用到的參數(shù)(對應(yīng)步驟102 )。
2.1計(jì)算JV,和^:將步驟101提取的耀斑爆發(fā)時間和爆發(fā)位置參數(shù)輸入程序 tim2carr.pro，執(zhí)行該程序，輸出耀斑爆發(fā)時間的太陽自轉(zhuǎn)周序號A^和耀斑爆發(fā) 位置的卡林頓經(jīng)度&，具體的執(zhí)行結(jié)果如下，其中下邊倒數(shù)第二列(10個字符) 和最后一列(9個字符)分別為自轉(zhuǎn)周序號和卡林頓經(jīng)度。
7701051457N30W90C0.11650.1749727.01227
7701311726N27E60M0.21651.12972253.30064
7702020215N26E39C0.11651.17973256.29756
7702020322N26E39C0.11651.18143255.68496
7702020425N28E36C0.61651.18303258.10894
7702091458N12E42C0.11651.45513154.15489
7702171339N19W54CO,l1651.74574145.53523
7702182344N23W77C0.21651.79768149.83363
0508232027N17W06C2.02033.57556204.53967
0508232158N12E85C3.52033.57788112.70466
0508232310N13E85C 1.32033.57971112.04400
0508240155N12E85C2.32033,58392110.52999
0508240407N12E84C 1.32033.58710110.38302
0508250431N09E80M6.52033.62460100.88577
0508260241N12E68C 1.32033.65849100.68275
0508261139N10E61C2.22033.67220102.74663
0508262152N10E50C2.42033.68783108.12247
0508272036畫E37C 1.52033.72259108.60839
0508281017歸E36M 1.72033.74351102.07632050831 1026N13W13C2.1 2033.85382 111.36379
因此根據(jù)以上數(shù)據(jù)，計(jì)算數(shù)據(jù)起始點(diǎn)的自轉(zhuǎn)周序號記為乂， iV。=1650;當(dāng)前 自轉(zhuǎn)周為第2033周，預(yù)報(bào)的自轉(zhuǎn)周為第2034周。 2.2計(jì)算每個自轉(zhuǎn)周的太陽活動平均緯度p,
第1650周只有一個耀斑，因此p，就是該耀斑的緯度，即仍65。=30°
第1651周的平均緯度^51根據(jù)步驟Bl中介紹的計(jì)算方法得到，
—0.2x10x27。 +0.1x1x26° +0.1x1x26。 +0.6x1x28。 +0.1x1x12。 +0.1x1x19° +0.2x1x23。 "1651— 0.2x10 + 0.1x1 + 0.1x1 + 0.6x1 + 0.1x1 + 0.1x1 + 0.2x1 同理
=26.16。
得到第2033周的平均緯度％。33 = 9.sr
2.3根據(jù)步驟2.2得到的平均緯度，可以計(jì)算出每個自轉(zhuǎn)周對應(yīng)的sin2 ^
2.4計(jì)算方程⑧中的疋值，即把每個耀斑的爆發(fā)時間換算成以天為單位。因?yàn)橐?個太陽自轉(zhuǎn)周2;為27.2753天，所以每個耀斑在所屬自轉(zhuǎn)周的天數(shù)為自轉(zhuǎn)周的 小數(shù)部分乘以27.2753天。
即，* = [M-*(iV,)]x7;,力x(iV,)表示自轉(zhuǎn)周取整數(shù)。所以第1650.17497周對應(yīng)的 耀斑爆發(fā)在第1650周的第0.17497x27.2753-4.77236天。同理可以得到每個耀斑對 應(yīng)的先值。
3、計(jì)算太陽赤道自轉(zhuǎn)角速度Q。、緯向較差修正常數(shù)S、以及太陽活動經(jīng)度帶初 始位置A。,和太陽活動經(jīng)度帶當(dāng)前位置A;1 (對應(yīng)步驟103 )。
3.1計(jì)算Q。、 5和A。,。此時，方程⑧中只有A。,、 Q。和S為未知數(shù)，其他均已計(jì) 算出或?yàn)橐阎?shù)(1^=14.1844, 7; = 27.2753 )。按步驟B3中的搜索方法搜索到 的一組最佳參數(shù)為Q, =14.59 (度/天)，丑,=1.43 (度/天)，A。1/;,=132 (度)。 3.2計(jì)算太陽活動經(jīng)度帶當(dāng)前位置八,,。根據(jù)方程②④⑤，Q，、 A,和A一,三個 參數(shù)可以計(jì)算出爆發(fā)最后一個耀斑時的經(jīng)度帶中心位置A,^A^一^283 (度)，另一條太陽活動經(jīng)度帶中心位置為103度。
4、 從步驟2.2中提取當(dāng)前太陽活動平均緯度p，i， =9.87° (對應(yīng)步驟104)。
5、 預(yù)報(bào)太陽活動經(jīng)度帶在下一自轉(zhuǎn)周的卡林頓位置(對應(yīng)步驟105)。 5.1計(jì)算太陽活動經(jīng)度帶在最近一 自轉(zhuǎn)周的轉(zhuǎn)動速度Q，，
3 = , -萬/"sin2^3 = 14.59-1.43sin2(9.87°) = 14.59-1.43x0，38 = 14.548 (度/天)， 太陽活動經(jīng)度帶在 一個自轉(zhuǎn)周的角速度的變化可以忽略，可以將此值視為下一 個自轉(zhuǎn)周的轉(zhuǎn)動速度，即Q2一-Q鵬。
5.2計(jì)算第2033周最后一天的經(jīng)度帶卡林頓位置A,w。太陽活動經(jīng)度帶在卡林 頓坐標(biāo)系中每天的漂移d鄉(xiāng),《婦=Q2。34/;-Qc =QM,-QC =14.548-14.1844 = 0.3636(度
/天)。第2033周最后一個耀斑的爆發(fā)時間是2005年8月31日10時，是第 2(B3,85382周，第2044.0周從2005年9月4日14時開始，由于后邊計(jì)算的都 是世界時00:00時刻的位置，所以我們將每個自轉(zhuǎn)周的起始時間定為該周起始 曰的00:00時刻，即第2044周從9月4曰00:00時刻起，從8月31曰10時到 9月4曰00:00時刻，經(jīng)過了3天14小時，在9月4曰00:00時刻，經(jīng)度帶的 卡林頓坐標(biāo)八,w為，A^"3 + 14/24)x0.3636 +八a =284.303 (度)。預(yù)報(bào)第2044周 一條經(jīng)度帶的卡林頓坐標(biāo)A阿一 =A2。4, ="0.3636°+284.303°, A:e
, 另一條經(jīng)度帶的卡林頓坐標(biāo)八(,+,=八2。44^為
A204442/; =(A2M4(Hp +180°) mod 360° =(A:x0.3 6° +284.303° +180°)mocB60° 0.3636。 +104.303° ，"
.
當(dāng)&取整數(shù)0， 1， 2...時對應(yīng)的是9月4， 5, 6....日00:00時刻經(jīng)度帶的卡 林頓坐標(biāo)。即在4日00:00時刻，太陽活動經(jīng)度帶位于284.303和104.303度處， 在5日00:00時刻,位于284.667和104.667度處。
6、將預(yù)報(bào)的太陽活動經(jīng)度帶中心位置的卡林頓坐標(biāo)換算成日面坐標(biāo)(對應(yīng)步驟106)。
6.1根據(jù)步驟5,2計(jì)算出在第2044周9月4日到10月1日每日世界時00:00的 太陽活動經(jīng)度帶中心位置卡林頓坐標(biāo)A^,和A^，其中在9月4日00:00時刻， 太陽活動經(jīng)度帶位于Aa,284.303度和A^,104.303度處；在9月5日00:00時 刻，太陽活動經(jīng)度帶位于八^+1)1=284.667度和八^+1)2=104.667度處。 6.2計(jì)算在第2044周9月4曰到IO月1曰每曰世界時00:00的曰面中心的卡林 頓坐標(biāo)。執(zhí)行程序tim2carr.pro，時間參數(shù)為每日的00:00時刻。以9月4曰為 例，執(zhí)行L=TIM2CARR('4-Sep-2005 00:00'),得到的1^=5.55491(度)，就是4曰 00:00時刻日面中心的卡林頓坐標(biāo)A^ =5.55491° 。執(zhí)行L=TIM2CARR('5-Sep-2005 00:00')，得到的1^=352.347(度)，就是5日00:00時刻日面中心的卡林頓坐標(biāo) A，)-352.347。。
6.3計(jì)算相鄰兩日00:00時刻太陽活動經(jīng)度帶的日面坐標(biāo)，確定太陽活動經(jīng)度帶
每天的日面經(jīng)度變化范圍。
4曰00:00時刻兩經(jīng)度帶的曰面坐標(biāo)為
Aw = AClU - Aa = 284.303° - :5.5，1° = 278.748° A'42 =~2-Aa =104.303°-5.55491° =98.748°
其中將278.748°化入[-180°, 180°] ，AiA1 =278.748° -360° =-81.252° 。
5曰00:00時刻兩經(jīng)度帶的曰面坐標(biāo)為 = Ac(4+I)1 -As(,+1) =284.667° -352.347° =-67.680° A*+i)2 =八*+1)2 -As("!) =104.667° - 352.347° = —247.68° =112.320°
因此，從4日00:00時刻到5日00:00時刻，太陽活動經(jīng)度帶中心位置的曰 面經(jīng)度變化范圍分別是(-81.252' ~-67.680°)和(98.748° -112.320°)。 6.4確定緯度變化范圍。平均緯度％。34=9.87。<10°，根據(jù)步驟D3中所述的方法， 得緯度變化范圍為0° ~ 9.87° +10° — 0° ~ 19.87° 。
6.5確定經(jīng)度帶的日面位置，或預(yù)報(bào)日面上最活躍的區(qū)域。根據(jù)步驟6.3和6.4 得到北半球在9月4日00:00時刻到5日00:00時刻，最活躍的區(qū)域?yàn)?—81.252° ~—67.680。， 0° ~ 19.87° )和(98.748° ~ 112.320°, 0°~19.87°)。同理可以預(yù) 報(bào)該自轉(zhuǎn)周每相鄰兩天的經(jīng)度帶日面變化范圍。
每過一個自轉(zhuǎn)周，在該自轉(zhuǎn)周內(nèi)如果有太陽活動發(fā)生，則加入新觀測數(shù)據(jù) 重復(fù)上述步驟作預(yù)報(bào)，如果沒有太陽活動發(fā)生，則估計(jì)該自轉(zhuǎn)周的太陽活動平 均緯度，只需從步驟104開始執(zhí)行，步驟IOI， 102，和103用上個自轉(zhuǎn)周執(zhí)行 的結(jié)果即可。
以上所述，僅為本發(fā)明的較佳實(shí)施例而已，并非用于限定本發(fā)明的保護(hù)范圍。
權(quán)利要求
1、一種太陽活動經(jīng)度帶預(yù)報(bào)方法，其特征是所述預(yù)報(bào)方法包括以下步驟，A、選取太陽X射線耀斑觀測數(shù)據(jù)；B、根據(jù)所述觀測數(shù)據(jù)計(jì)算出太陽自轉(zhuǎn)的相關(guān)參數(shù)和太陽活動經(jīng)度帶起始位置的卡林頓坐標(biāo)7，進(jìn)而計(jì)算出太陽活動經(jīng)度帶當(dāng)前位置的卡林頓坐標(biāo)；所述太陽自轉(zhuǎn)的相關(guān)參數(shù)包括太陽赤道自轉(zhuǎn)角速度Ω0和緯向較差修正常數(shù)B；C、根據(jù)所述太陽自轉(zhuǎn)的相關(guān)參數(shù)、太陽活動經(jīng)度帶的起始位置和當(dāng)前位置計(jì)算出下一個太陽自轉(zhuǎn)周或更長時間內(nèi)的太陽活動經(jīng)度帶中心位置的卡林頓坐標(biāo)；D、將所述預(yù)報(bào)的太陽活動經(jīng)度帶中心位置的卡林頓坐標(biāo)換算成日面坐標(biāo)，得到所述預(yù)報(bào)的太陽活動經(jīng)度帶的日面位置。
2、 根據(jù)權(quán)利要求l所述的太陽活動經(jīng)度帶預(yù)報(bào)方法，其特征在于所述步 驟A中，從太陽X射線耀斑觀測數(shù)據(jù)中提取耀斑級別、爆發(fā)時間和爆發(fā)位置參 數(shù)；所述步驟B中計(jì)算太陽活動經(jīng)度帶的當(dāng)前位置卡林頓坐標(biāo)的方法為 Bl、根據(jù)從所述步驟A中提取的每個耀斑的爆發(fā)時間和爆發(fā)位置參數(shù)， 計(jì)算該耀斑爆發(fā)時間對應(yīng)的自轉(zhuǎn)周序號和爆發(fā)位置的卡林頓經(jīng)度，其中將第/自 轉(zhuǎn)周發(fā)生的第y個耀斑的卡林頓經(jīng)度記為& ,該耀斑爆發(fā)時間對應(yīng)的自轉(zhuǎn)周序 號記為W,.,該耀斑的爆發(fā)時間換算成以天為單位記為"根據(jù)從所述步驟A中 提取的耀斑級別參數(shù)，以X射線耀斑級別為權(quán)重，.計(jì)算每個自轉(zhuǎn)周的太陽活動平均緯度仍；將計(jì)算數(shù)據(jù)起始點(diǎn)的自轉(zhuǎn)周序號記為 ；B2、把步驟B1得到的數(shù)據(jù)a, l a， iV。和iV,，以及已知常量7;和a.，代入表達(dá)式<formula>formula see original document page 3</formula>min(耳—[A^ +TC 2 (Qo —Qc —Ssin2仍)+—Qc —_Ssin2仍)]mod360。& 一[VV +rc ^(Qo —Qc —5sin2仍)+—Qc -丑sin2^)±180' ]mod360該式中，s代表耀斑位置和太陽活動經(jīng)度帶中心位置的平均偏離程度，A。,為太 陽活動經(jīng)度帶起始位置的卡林頓坐標(biāo)，Q。和S是待定常數(shù)，iV是耀斑總數(shù)，'±' 當(dāng)A。^180。時取'+，，反之取'-，，mod360。表示劃入0。-360。范圍內(nèi)；B3、在給定搜索區(qū)間內(nèi)以給定步長搜索滿足e為最小的一組參數(shù)Q。， S和 A。,，并相應(yīng)記為Q一， A,和A,;其中，Q。的搜索區(qū)間為[13.50, 15.00]，步 長為0.01; ￡的搜索區(qū)間為
，步長為0.01; AQ1的搜索區(qū)間為
，步長為r;B4、根據(jù)步驟B3中得到的參數(shù)Q季，A,和A，，以及方程fi,Q。-SsinV,、 M = rc Z (Q, —Qc) + — 。c)、八,.,=(A。,十M)mod360。計(jì)算出太陽活動經(jīng)度帶的,=A'0當(dāng)前位置的卡林頓坐標(biāo)A,."其中A是第/周某個耀斑爆發(fā)時間，換算成以天為單 位；所述步驟C中計(jì)算下一個太陽自轉(zhuǎn)周內(nèi)太陽活動經(jīng)度帶中心位置卡林頓坐 標(biāo)的方法為Cl、設(shè)當(dāng)前卡林頓自轉(zhuǎn)周為第/周，計(jì)算當(dāng)前周的太陽活動經(jīng)度帶的轉(zhuǎn)動角 速度Q,Q。,-^sin2^,將下一周即第i+l周的太陽活動經(jīng)度帶轉(zhuǎn)動角速度A,取作當(dāng)前周的太陽活動經(jīng)度帶的轉(zhuǎn)動角速度Q,，得到第i+l周的太陽活動 經(jīng)度帶的轉(zhuǎn)動角速度Q,+^，即Q(^p-Q,，其中p代表預(yù)報(bào)；C2、將第i+l周的起始時間定為該周起始日的00:00時刻；根據(jù)公式 A^nd =仏 + - A)} mod 360°計(jì)算出第i周最后一天的太陽活動經(jīng)度帶位置的卡 林頓坐標(biāo)A,一，其中yt代表第i周最后一個耀斑的爆發(fā)時間到第i+l周起始曰的 00:00時刻經(jīng)過的時間，以天為單位；根據(jù)公式A(,+，r仏,w+^VD廣Qc》mod36(r,得到第i+l周一個太陽活動經(jīng)度帶中心位 置的卡林頓坐標(biāo)，its
,當(dāng)A取整數(shù)0， 1， 2...時對應(yīng)的是第i+l周的 第1， 2, 3，....日00:00時刻太陽活動經(jīng)度帶的卡林頓坐標(biāo)；根據(jù)公式 A(,+1)i2p = (A(,孝p +18(T)m0d36(r ，得到第i+l周另 一個太陽活動經(jīng)度帶中心位置的 卡林頓坐標(biāo)；所述步驟D中將所述太陽活動經(jīng)度帶中心位置的卡林頓坐標(biāo)換算成日面坐 標(biāo)的方法為Dl、記日輪中心經(jīng)度為0°,從中心沿0諱度向東到太陽背面的中心為 0° -180°，從中心沿O緯度向西至太陽背面的中心為0°~180°;根據(jù)步驟C2所述 的公式A, ={A,w"(Q(,+1);)_Qc)}m0d36(X和A(孝尸=(A( +18(T)m0d36(r可以計(jì)算出第i+l周每曰00:00時太陽活動經(jīng)度帶的卡林頓坐標(biāo)A^和A^;計(jì)算出第 i+l周每日00:00時日面中心的卡林頓坐標(biāo)A說；根據(jù)公式Ai41 = Acw-Aa和 A.、,2=At,2-A雙計(jì)算出每天世界時00:00時的太陽活動經(jīng)度帶中心位置的曰面經(jīng) 度A^和A^，其中，角標(biāo)^代表曰面坐標(biāo)，^代表第ifc天；D2 、根據(jù)連續(xù)兩個世界時對應(yīng)的日面經(jīng)度A,,,和A^+d,之間的變化范圍 l八^ -八一wl ，得到 一條太陽活動經(jīng)度帶中心位置的每天曰面經(jīng)度變化范圍為 A, ~A#+1)1;同理，另一條太陽活動經(jīng)度帶中心位置的每天曰面經(jīng)度變化范圍為A^2 A卓+i)2 ，D3、緯度范圍取為該自轉(zhuǎn)周的平均緯度R的上下10度，即當(dāng)l仍 10。時， 緯度范圍相應(yīng)為0°~^.+10°或^-10°~0°，其中，爐,的正、負(fù)分別表示太陽的北半球和南半球；D4、根據(jù)經(jīng)度帶中心位置的每天日面經(jīng)度變化范圍及緯度范圍，得到北半 球兩個經(jīng)度帶中心位置的每天日面變化范圍，經(jīng)度變化范圍為A^ A,(一,和 A_rt2~A.，)2;緯度變化范圍分兩種情況為在p,〉l(T時，兩經(jīng)度帶中心位置的緯 度變化范圍均為^-10° ^+10°; 0、A《10。時，緯度變化范圍均為0°~^+10° ;這樣得出北半球兩個經(jīng)度帶中心位置的每天曰面變化范圍為(A,H ~ A#+1)l , a—10° ~ a+10° )和(~ A,(i+1)2 ,a—10° 仍+10° ) ( r>1(t )，或 (A,h A卓,，0。 r+10。)和(A,A2~A,(t+I)2，0°~a+10° ) ( 0° 《10° );以及得出南 半球的兩個經(jīng)度帶中心位置的每天變化范圍為(A',w A',(w,， r-1(t a+10。)和(A:、i2 A',,)2， ^一i(t a+i(r ) ( "〈一i(t ),或(A必 A',(w),wi(T (r)和(A'ii2~A',(t+1)2，A-10° 0° ) ( -10 0° )。
3、 根據(jù)權(quán)利要求2所述的太陽活動經(jīng)度帶預(yù)報(bào)方法，其特征在于所述步 驟A中，若在太陽活動極小期間沒有太陽耀斑，則選取太陽黑子觀測數(shù)據(jù)。
4、 根據(jù)權(quán)利要求3所述的太陽活動經(jīng)度帶預(yù)報(bào)方法，其特征在于所述步 驟A中，選取太陽黑子觀測數(shù)據(jù)時,從太陽黑子觀測數(shù)據(jù)中提取黑子的面積、出 現(xiàn)時間和出現(xiàn)位置參數(shù)。
5、 根據(jù)權(quán)利要求4所述的太陽活動經(jīng)度帶預(yù)報(bào)方法，其特征在于，所述步 驟A中，選取太陽黑子觀測數(shù)據(jù)時,所述步驟B1中以太陽黑子面積為權(quán)重計(jì)算 平均緯度a。
6、 根據(jù)權(quán)利要求5所述的太陽活動經(jīng)度帶預(yù)報(bào)方法，其特征在于:每過一 個自轉(zhuǎn)周，加入新太陽活動觀測數(shù)據(jù)，并重復(fù)步驟A D,對下一個自轉(zhuǎn)周進(jìn)行 預(yù)報(bào)。
7、 根據(jù)權(quán)利要求6所述的太陽活動經(jīng)度帶預(yù)報(bào)方法，其特征在于如果當(dāng) 前自轉(zhuǎn)周有太陽X射線耀斑觀測數(shù)據(jù)，則所述加入的新太陽活動觀測數(shù)據(jù)為太 陽X射線耀斑觀測數(shù)據(jù)；如果當(dāng)前自轉(zhuǎn)周沒有太陽X射線耀斑觀測數(shù)據(jù)，則所 述加入的新太陽活動觀測數(shù)據(jù)為太陽黑子觀測數(shù)據(jù)；如果當(dāng)前自轉(zhuǎn)周既沒有太 陽X射線耀斑觀測數(shù)據(jù)，也沒有太陽黑子觀測數(shù)據(jù)，則不加入的新太陽活動觀 測數(shù)據(jù)，省略步驟A和B;此時，用一年內(nèi)的太陽活動平均緯度代替該自轉(zhuǎn)周 的平均緯度a，重復(fù)步驟C和D，獲得預(yù)報(bào)數(shù)據(jù)。
全文摘要
本發(fā)明公開了一種太陽活動經(jīng)度帶的預(yù)報(bào)方法。該方法主要包括以下步驟首先選取太陽X射線耀斑數(shù)據(jù)(在太陽活動極小期間，沒有太陽耀斑時就用太陽黑子數(shù)據(jù))；根據(jù)選取的觀測數(shù)據(jù)確定太陽自轉(zhuǎn)的相關(guān)參數(shù)和太陽活動經(jīng)度帶的當(dāng)前位置(卡林頓坐標(biāo))；根據(jù)確定的太陽自轉(zhuǎn)參數(shù)和太陽活動經(jīng)度帶的當(dāng)前位置預(yù)報(bào)一定時間內(nèi)的太陽活動經(jīng)度帶中心位置(卡林頓坐標(biāo))；將預(yù)報(bào)的太陽活動經(jīng)度帶中心位置的卡林頓坐標(biāo)換算成日面坐標(biāo)，預(yù)報(bào)太陽活動經(jīng)度帶的日面位置。該預(yù)報(bào)不僅豐富了太陽活動中長期的預(yù)報(bào)內(nèi)容，也為太陽活動短期預(yù)報(bào)提供一定的參考價值。此外，這種預(yù)報(bào)方法的可操作性強(qiáng)，降低了預(yù)報(bào)員的主觀性。
文檔編號G01W1/10GK101430388SQ20081018045
公開日2009年5月13日 申請日期2008年11月27日 優(yōu)先權(quán)日2008年11月14日
發(fā)明者張麗云, 杜占樂, 王華寧, 晗 賀 申請人:中國科學(xué)院國家天文臺