地磁低點(diǎn)位移

以甘肅16個(gè)地磁臺(tái)站的觀測(cè)結(jié)果為數(shù)據(jù)的地磁低點(diǎn)位移地震預(yù)報(bào)方法探究目錄我省地磁觀測(cè)系統(tǒng)的臺(tái)站目錄及數(shù)據(jù)代碼地球磁場(chǎng)的構(gòu)成及變化地磁低點(diǎn)位移預(yù)報(bào)地震方法地磁低點(diǎn)位移法預(yù)報(bào)地震存在的問(wèn)題地磁低點(diǎn)位移法的程序設(shè)計(jì)1.臺(tái)站目錄與數(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蘆陽(yáng)'62146_5_31230_分.TXT'11松山'62147_2_31230_分.TXT'紅砂無(wú)數(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'臺(tái)站代碼測(cè)點(diǎn)編號(hào)垂直分量采樣周期甘肅省地磁觀測(cè)臺(tái)站分布圖2.地球磁場(chǎng)的構(gòu)成及變化磁場(chǎng)穩(wěn)定磁場(chǎng)，內(nèi)源場(chǎng)占99%，其余是外源場(chǎng)變化磁場(chǎng)，外源場(chǎng)占66%，其余是內(nèi)源場(chǎng)按穩(wěn)定性劃分磁場(chǎng)內(nèi)源場(chǎng)外源磁場(chǎng)按來(lái)源劃分異常場(chǎng)基本磁場(chǎng)非偶極子場(chǎng)地心偶極子場(chǎng)局部磁異常，源于地殼淺處區(qū)域磁異常，源于地殼深部地球的地磁軸Sm-Nm與自轉(zhuǎn)軸S-N南北極垂直分量較大赤道附近水平分量較大南北極附近垂直分量較大甘肅，水平地磁北極36N0甘肅地區(qū)垂直分量大于水平分量磁場(chǎng)的變化-擾動(dòng)變化對(duì)觀測(cè)有影響平靜變化擾動(dòng)變化太陽(yáng)風(fēng)弓激波極尖區(qū)等離子體磁氣圈;

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

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

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

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

據(jù)丁鑒?！瘛瘛?5.00突變界線的性質(zhì)地磁低點(diǎn)位移法一般只預(yù)報(bào)6、7級(jí)強(qiáng)震。在少震區(qū)、地震平靜時(shí)期且臺(tái)網(wǎng)密度較大的地區(qū)，有時(shí)可預(yù)報(bào)5級(jí)地震。分界線也是地震、地磁、形變、流體短臨異常突出地區(qū)，即屬于短期地球物理場(chǎng)分界線。地震沿地磁低點(diǎn)位移分界線遷移。低點(diǎn)位移反映地質(zhì)構(gòu)(板塊塊體)造微動(dòng)態(tài)活動(dòng)。低點(diǎn)位移反映孕震內(nèi)外環(huán)境因素的變化，亦即內(nèi)外機(jī)理耦合作用。每天一個(gè)值，基本與觀測(cè)臺(tái)站地方是午時(shí)一致(倫敦時(shí)間5點(diǎn)=蘭州時(shí)間12點(diǎn))，由東向西按4分鐘/經(jīng)度延遲。以地磁垂直分量(Z)日變化極小時(shí)間為分析對(duì)象。變化磁場(chǎng)的感應(yīng)場(chǎng)的異常反映到地磁日變上來(lái),這種異常變化在地磁垂直分量上最明顯。當(dāng)日極大值當(dāng)日極小值出現(xiàn)的時(shí)間(倫敦時(shí)間計(jì)時(shí))采用資料第九頁(yè)圖+7+6+0背景場(chǎng)時(shí)間約為當(dāng)?shù)貢r(shí)間12點(diǎn)正常場(chǎng)在當(dāng)?shù)貢r(shí)間12點(diǎn)正常場(chǎng)在當(dāng)?shù)貢r(shí)間12點(diǎn)蘭州當(dāng)?shù)貢r(shí)間蘭州當(dāng)?shù)貢r(shí)間蘭州當(dāng)?shù)貢r(shí)間臺(tái)站太少，無(wú)法繪出低點(diǎn)高點(diǎn)的大區(qū)域。我省轄區(qū)是條狀分布的，繪出的突變界線不準(zhǔn)確格林尼治時(shí)間4.基于我省資料的地磁低點(diǎn)位移法預(yù)報(bào)地震存在的問(wèn)題我省管轄區(qū)域?yàn)殚L(zhǎng)條形，不利于突變界線的繪制。按照統(tǒng)計(jì)，臺(tái)站數(shù)量>50，才能畫出準(zhǔn)確的突變界線，而我省臺(tái)站只有16個(gè)，過(guò)于稀少。震例不足，導(dǎo)致研究缺乏數(shù)據(jù)。部分臺(tái)站維護(hù)不良導(dǎo)致缺數(shù)、無(wú)數(shù)據(jù)。地磁低點(diǎn)位移方法不能預(yù)報(bào)準(zhǔn)確的發(fā)震地點(diǎn)，只能指出發(fā)震區(qū)在突變分界線附近。存在干擾因素（山崩、暴雨、磁暴、洪水也可導(dǎo)致異常）導(dǎo)致誤報(bào)。國(guó)外地震也會(huì)造成國(guó)內(nèi)異常。異常日附近并不發(fā)震，80%地震發(fā)生在異常過(guò)后第27天或41天的前后4天，這使得預(yù)報(bào)的時(shí)間的準(zhǔn)確性打了折扣。解決思路新建臺(tái)站增大密度不現(xiàn)實(shí)。可把青海、寧夏、陜西的地磁臺(tái)觀測(cè)數(shù)據(jù)整合起來(lái)，一是可以形成一個(gè)完整的閉合區(qū)域(矩形或圓形)利于繪制突變界線，二是可以湊出足夠多的臺(tái)站數(shù)量保證突變界線的精確度。加強(qiáng)臺(tái)站的維護(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ù)使用的是格林尼治時(shí)間，加上7小時(shí)之后是蘭州當(dāng)?shù)貢r(shí)間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','蘆陽(yáng)20','松山21','古豐22','黃羊23','橫梁24','臨夏25','肅北26','肅北27'};%27個(gè)[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ì)說(shuō)明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ù)(臺(tái)站數(shù))*工作天數(shù)fork=1:nn;%num(2)=27個(gè)文件夾=臺(tái)站數(shù)QQ=strcat(path,Q{k});SS=load(QQ);DD=SS(:,2:end);DD(DD==999999.0)=NaN;%把矩陣DAE中所有999999的數(shù)值用NaN來(lái)代替

[yy,XX(k,:)]=min(DD.');%輸出單個(gè)臺(tái)站從年初到工作日每日的地磁地點(diǎn)位移時(shí)間序列endXX(XX==1)=NaN;XX=XX+420;%格林尼治時(shí)間轉(zhuǎn)化為北京時(shí)間7*60=420子程序inpfunctiongg=draw(XX,ww,nn,rr)%UNTITLED2此處顯示有關(guān)此函數(shù)的摘要%此處顯示詳細(xì)說(shuō)明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('儀器序號(hào)','FontName','TimesNewRoman','FontSize',16,'FontWeight','Bold');ylabel('地磁低點(diǎn)位移時(shí)間','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\中國(guó)縣界_譚平玉\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個(gè)臺(tái)站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)(時(shí))','FontName','TimesNewRoman','FontSize',15,'FontWeight','Bold');holdonav=mean(XX(noo,1:r