地磁低點(diǎn)位移

以甘肅16個地磁臺站的觀測結(jié)果為數(shù)據(jù)的地磁低點(diǎn)位移地震預(yù)報(bào)方法探究目錄我省地磁觀測系統(tǒng)的臺站目錄及數(shù)據(jù)代碼地球磁場的構(gòu)成及變化地磁低點(diǎn)位移預(yù)報(bào)地震方法地磁低點(diǎn)位移法預(yù)報(bào)地震存在的問題地磁低點(diǎn)位移法的程序設(shè)計(jì)1.臺站目錄與數(shù)據(jù)代碼1蘭州'62001_2_31230_分.TXT''62001_8_31230_分.TXT''62001_9_31230_分.TXT'2山丹'62004_1_31230_分.TXT'3天水'62012_4_31230_分.TXT''62012_5_31230_分.TXT''62012_8_31230_分.TXT'4嘉峪關(guān)'62015_1_31230_分.TXT''62015_5_31230_分.TXT''62015_6_31230_分.TXT''62015_8_31230_.TXT'5合作'62052_2_31230_分.TXT''62052_3_31230_分.TXT'6瑪曲'62097_2_31230_分.TXT''62097_3_31230_分.TXT'7舟曲'62098_3_31230_分.TXT''62098_8_31230_分.TXT'8英鴿'62142_3_31230_分.TXT'9寺灘'62144_4_31230_分.TXT'10蘆陽'62146_5_31230_分.TXT'11松山'62147_2_31230_分.TXT'紅砂無數(shù)據(jù)12古豐'62149_3_31230_分.TXT'13黃羊'62150_3_31230_分.TXT'14橫梁'62151_2_31230_分.TXT'15臨夏'62171_1_31230_分.TXT'16肅北'62239_2_31230_分.TXT''62239_3_31230_分.TXT'臺站代碼測點(diǎn)編號垂直分量采樣周期甘肅省地磁觀測臺站分布圖2.地球磁場的構(gòu)成及變化磁場穩(wěn)定磁場，內(nèi)源場占99%，其余是外源場變化磁場，外源場占66%，其余是內(nèi)源場按穩(wěn)定性劃分磁場內(nèi)源場外源磁場按來源劃分異常場基本磁場非偶極子場地心偶極子場局部磁異常，源于地殼淺處區(qū)域磁異常，源于地殼深部地球的地磁軸Sm-Nm與自轉(zhuǎn)軸S-N南北極垂直分量較大赤道附近水平分量較大南北極附近垂直分量較大甘肅，水平地磁北極36N0甘肅地區(qū)垂直分量大于水平分量磁場的變化-擾動變化對觀測有影響平靜變化擾動變化太陽風(fēng)弓激波極尖區(qū)等離子體磁氣圈;

磁層T=24hT=25h[1]太陰日：月球中心連續(xù)兩次通過地球上同一子午線所需的時間，約為24小時50分鐘。這是以月球?yàn)閰⒖键c(diǎn)衡量地球的自轉(zhuǎn)周期。山丹臺磁暴高壓直流輸電影響3.地磁低點(diǎn)位移地震預(yù)報(bào)方法“地磁低點(diǎn)位移異?！爆F(xiàn)象指的是一個大區(qū)域臺站的地磁垂直分量日變化極小值出現(xiàn)的時間（簡稱低點(diǎn)時間）,明顯地與另一個大區(qū)域臺站的低點(diǎn)時間不同,而每個大區(qū)域內(nèi)部低點(diǎn)時間又基本一致,兩個區(qū)域之間有明顯的突變分界線，并且兩個區(qū)域之間的低點(diǎn)時間相差在兩小時以上。地磁日變地點(diǎn)時間正常背景場：垂直分量日變極小值出現(xiàn)時間（簡稱低點(diǎn)時間）出現(xiàn)在當(dāng)?shù)貢r間12點(diǎn)左右。如果偏離此時間過多，則認(rèn)為單臺低點(diǎn)時間發(fā)生位移，從地磁學(xué)理論分析這并不是異常，而是正常狀態(tài)的時間偏離[2]

。丁鑒海對167次震例進(jìn)行統(tǒng)計(jì)之后得出，地震發(fā)生在突變分界線附近,發(fā)震時間在異常發(fā)生日之后第27天或41天的前后4天內(nèi)

，地點(diǎn)在低點(diǎn)位移突變分界線附近。預(yù)報(bào)后再次出現(xiàn)的臨震異常不再預(yù)報(bào)，但可對預(yù)報(bào)參數(shù)進(jìn)行修正。[2]

丁鑒海.地磁日變地震預(yù)報(bào)方法及其震例研究.地震出版社2001年11月14日昆侖山口西8.1級地震日變極小值在空間分布上異常特征—低點(diǎn)位移現(xiàn)象（2001年10月磁場突變分界線在南北帶以西地區(qū)）

據(jù)丁鑒?！瘛瘛?5.00突變界線的性質(zhì)地磁低點(diǎn)位移法一般只預(yù)報(bào)6、7級強(qiáng)震。在少震區(qū)、地震平靜時期且臺網(wǎng)密度較大的地區(qū)，有時可預(yù)報(bào)5級地震。分界線也是地震、地磁、形變、流體短臨異常突出地區(qū)，即屬于短期地球物理場分界線。地震沿地磁低點(diǎn)位移分界線遷移。低點(diǎn)位移反映地質(zhì)構(gòu)(板塊塊體)造微動態(tài)活動。低點(diǎn)位移反映孕震內(nèi)外環(huán)境因素的變化，亦即內(nèi)外機(jī)理耦合作用。每天一個值，基本與觀測臺站地方是午時一致(倫敦時間5點(diǎn)=蘭州時間12點(diǎn))，由東向西按4分鐘/經(jīng)度延遲。以地磁垂直分量(Z)日變化極小時間為分析對象。變化磁場的感應(yīng)場的異常反映到地磁日變上來,這種異常變化在地磁垂直分量上最明顯。當(dāng)日極大值當(dāng)日極小值出現(xiàn)的時間(倫敦時間計(jì)時)采用資料第九頁圖+7+6+0背景場時間約為當(dāng)?shù)貢r間12點(diǎn)正常場在當(dāng)?shù)貢r間12點(diǎn)正常場在當(dāng)?shù)貢r間12點(diǎn)蘭州當(dāng)?shù)貢r間蘭州當(dāng)?shù)貢r間蘭州當(dāng)?shù)貢r間臺站太少，無法繪出低點(diǎn)高點(diǎn)的大區(qū)域。我省轄區(qū)是條狀分布的，繪出的突變界線不準(zhǔn)確格林尼治時間4.基于我省資料的地磁低點(diǎn)位移法預(yù)報(bào)地震存在的問題我省管轄區(qū)域?yàn)殚L條形，不利于突變界線的繪制。按照統(tǒng)計(jì)，臺站數(shù)量>50，才能畫出準(zhǔn)確的突變界線，而我省臺站只有16個，過于稀少。震例不足，導(dǎo)致研究缺乏數(shù)據(jù)。部分臺站維護(hù)不良導(dǎo)致缺數(shù)、無數(shù)據(jù)。地磁低點(diǎn)位移方法不能預(yù)報(bào)準(zhǔn)確的發(fā)震地點(diǎn)，只能指出發(fā)震區(qū)在突變分界線附近。存在干擾因素（山崩、暴雨、磁暴、洪水也可導(dǎo)致異常）導(dǎo)致誤報(bào)。國外地震也會造成國內(nèi)異常。異常日附近并不發(fā)震，80%地震發(fā)生在異常過后第27天或41天的前后4天，這使得預(yù)報(bào)的時間的準(zhǔn)確性打了折扣。解決思路新建臺站增大密度不現(xiàn)實(shí)?？砂亚嗪?、寧夏、陜西的地磁臺觀測數(shù)據(jù)整合起來，一是可以形成一個完整的閉合區(qū)域(矩形或圓形)利于繪制突變界線，二是可以湊出足夠多的臺站數(shù)量保證突變界線的精確度。加強(qiáng)臺站的維護(hù)管理，使得數(shù)據(jù)不要間斷。改進(jìn)地磁低點(diǎn)位移地震預(yù)報(bào)的程序(軟件)，使其功能更完善。clc;clear;vv=load('D:\MapSIS\H01結(jié)果_數(shù)據(jù)\62001_2_31230_分.TXT');vv=num2str(vv(end,1));v1=str2double(vv(1:4));v2=str2double(vv(5:6));v3=str2double(vv(7:8));vv=[v1v2v3];%******地磁數(shù)據(jù)使用的是格林尼治時間，加上7小時之后是蘭州當(dāng)?shù)貢r間path='D:\MapSIS\H01結(jié)果_數(shù)據(jù)\';%輸入路徑**ww=[2014,07,16];%輸入繪圖日日期**繪圖日期必須在工作日期之前%***********************************rr=datenum(ww)-datenum([2014,1,1])+1;[XX,nn]=inp(path,vv);gg=draw(XX,ww,nn,rr);%%地圖模塊---------------------nl={'蘭州1','蘭州2','蘭州3','山丹4','天水5','天水6','天水7','嘉峪關(guān)8','嘉峪關(guān)9',...'嘉峪關(guān)10','嘉峪關(guān)11','合作12','合作13','瑪曲14','瑪曲15','舟曲16','舟曲17','英鴿18',...'寺灘19','蘆陽20','松山21','古豐22','黃羊23','橫梁24','臨夏25','肅北26','肅北27'};%27個[ax,tl]=drap(XX,gg,nl,rr);%%------------------------noo=14;tt=dras(XX,noo,nl,ww,rr);5.地磁低點(diǎn)位移法的序代碼function[XX,nn,yy]=inp(path,vv)%UNTITLED2此處顯示有關(guān)此函數(shù)的摘要%此處顯示詳細(xì)說明FileName=dir(fullfile(path,'*.TXT'));FileName=struct2cell(FileName);Q=FileName(1,:);%提取文件名矩陣的第一行信息nn=length(Q);pp=datenum(vv)-datenum([2014,1,1])+1;%工作天數(shù)XX=zeros(nn,pp);%文件夾數(shù)(臺站數(shù))*工作天數(shù)fork=1:nn;%num(2)=27個文件夾=臺站數(shù)QQ=strcat(path,Q{k});SS=load(QQ);DD=SS(:,2:end);DD(DD==999999.0)=NaN;%把矩陣DAE中所有999999的數(shù)值用NaN來代替

[yy,XX(k,:)]=min(DD.');%輸出單個臺站從年初到工作日每日的地磁地點(diǎn)位移時間序列endXX(XX==1)=NaN;XX=XX+420;%格林尼治時間轉(zhuǎn)化為北京時間7*60=420子程序inpfunctiongg=draw(XX,ww,nn,rr)%UNTITLED2此處顯示有關(guān)此函數(shù)的摘要%此處顯示詳細(xì)說明figure;plot(1:nn,XX(:,rr),'*-b');set(gca,'xtick',1:1:nn,'xlim',[1nn],'ytick',(4:9)*120,'ylim',[13]*360,'yticklabel',...{'8點(diǎn)','10點(diǎn)','12點(diǎn)','14點(diǎn)','16點(diǎn)','18點(diǎn)'},'FontName','TimesNewRoman','FontSize',15);xlabel('儀器序號','FontName','TimesNewRoman','FontSize',16,'FontWeight','Bold');ylabel('地磁低點(diǎn)位移時間','FontName','TimesNewRoman','FontSize',16,'FontWeight','Bold');holdon;gridon;plot([1nn],[11]*max(XX(:,rr)),'r',[1nn],[11]*min(XX(:,rr)),'r');md=max(XX(:,rr))-min(XX(:,rr));uu=strcat('md=',num2str(md),'分鐘');text(nn/2,max(XX(:,rr)),uu,'FontName','TimesNewRoman','FontSize',16,'FontWeight','Bold');ifmd<120er='\color{red}md＜2HOUR';elseer='\color{red}md≥2HOUR';endgg=strcat('DATE:',num2str(ww(1)),'/',num2str(ww(2)),'/',num2str(ww(3)),'@',er);title(gg,'FontName','宋體','FontSize',19,'FontWeight','bold');%'ArialBlack'子程序drawfunction[ax,tl]=drap(XX,gg,nl,rr)figure('unit','normalized','position',[0,0,1,1]);ax=worldmap([32.242.9],[92.5109]);setm(gca,'MapParallels',0);%設(shè)置地圖軸線模式,'地圖(的軸線)平行'setm(gca,'MLineLocation',2,'MLabelLocation',2,...%設(shè)置經(jīng)度'PLineLocation',2,'PLabelLocation',2)%設(shè)置緯度間隔、緯度標(biāo)簽間隔geoshow('G:\20140306\map\中國縣界_譚平玉\bou2_4p.shp','FaceColor',[011]);rivers=shaperead('worldrivers','UseGeoCoords',true);%加載河流geoshow(rivers,'Color','blue')tightmap;set(gcf,'color','w');%圖形底色設(shè)為白色,默認(rèn)為灰色,f_figure;w_white;a_axisaxisoff;%關(guān)閉圖片邊框framemon;%添加地圖經(jīng)緯網(wǎng)的外邊框textm(42.7,97.5,gg,'FontName','ArialBlack','FontSize',16,'FontWeight','bold');xl=[103.8,101.04,104.92,98.22,102.91,102.04,104.37,103.23,103.88,...

104.15,103.49,102.83,103.06,103.33,103.26,94.9,];yl=[36.1,38.78,34.57,39.81,35.01,34.02,33.8,36.91,37.23,...

37.05,37.11,37.43,37.37,37.3,35.63,39.53,];tl={nl{1},nl{4},nl{5},nl{8},nl{12},nl{14},nl{16},nl{18},...nl{19},nl{20},nl{21},nl{22},nl{23},nl{24},nl{25},nl{26}};%共16個臺站XX=XX(:,rr)';XX=[XX(1),XX(4),XX(5),XX(8),XX(12),XX(14),XX(16),XX(18),...XX(19),XX(20),XX(21),XX(22),XX(23),XX(24),XX(25),XX(26)];ll=fix(XX./60);%求整數(shù)kl=rem(XX,60);%求余數(shù)fork=1:16;plotm(yl(k),xl(k),'o','MarkerSize',3,'MarkerFaceColor',[110]);textm(yl(k)*0.997,xl(k),tl{k},'FontSize',9);ifkl(k)<10ml=strcat(num2str(ll(k)),':0',num2str(kl(k)));elseml=strcat(num2str(ll(k)),':',num2str(kl(k)));endtextm(yl(k)*1.003,xl(k),ml,'FontName','ArialBlack','FontSize',12,'color','r');end子程序drapfunction[tt,tr,av]=dras(XX,noo,nl,ww,rr)figure('unit','normalized','position',[0,0,1,1]);tt=datenum(2014,1,1):datenum(ww);%set(gcf,'Position',[400,100,300,300],'color','w')tr=[tt(1),tt(end)];subplot(2,1,1,'position',[0.050.540.9380.4]);%[leftbottomwidthheight]plot(tt,XX(noo,1:rr))xlabel('日期(天)','FontName','TimesNewRoman','FontSize',15,'FontWeight','Bold');ylabel('日地磁低點(diǎn)(時)','FontName','TimesNewRoman','FontSize',15,'FontWeight','Bold');holdonav=mean(XX(noo,1:r