2023航空重力測量技術(shù)規(guī)范

航空重力測量技術(shù)規(guī)范目 次前言 IV1范圍..............................................................................12規(guī)范性引用文件....................................................................13術(shù)語和定義........................................................................14總則..............................................................................25技術(shù)設(shè)計..........................................................................25.1 25.2 25.3 25.4 35.5 35.6 35.7 55.8 55.9 55.10 55.11 66儀器設(shè)備..........................................................................76.1 76.2 76.3 76.4 7量行野工作 7測設(shè)安裝 7GNSS基布設(shè) 8地重基建立 8航重儀溫試 9導(dǎo)定儀試 9測測飛行 9輔飛測量 11航重原資錄 12原資現(xiàn)檢驗 12巖礦石本集與度數(shù)定 13據(jù)理圖編制 14數(shù)處理 14I基圖編制 20質(zhì)評價 21數(shù)處報編寫 24測成驗收 24空力料斷釋 24解分類 24推解基要求 24區(qū)重場釋 25局重異解釋 25航重推解件 26果告寫評審 26附錄A（料附） 航空力量作計寫提綱 29附錄B（范附） 航空力量外合度評方法 32附錄C（料附） 航空力檢結(jié)記表 34附錄D（料附） 航空力量量計表 35附錄E（料附） 航空力量果告寫提綱 36II航空重力測量技術(shù)規(guī)范范圍本文件規(guī)定了航空重力測量的技術(shù)設(shè)計、儀器設(shè)備、測量飛行與野外作業(yè)、數(shù)據(jù)處理與圖件編制、航空重力資料推斷解釋與成果報告編寫等技術(shù)要求。（GB/T18314全球定位系統(tǒng)(GPS)測量規(guī)范DZ/T0004重力調(diào)查技術(shù)規(guī)定（1：50000）DZ/T0069地球物理勘查圖圖式圖例及用色標準DZ/T0171-1997大比例尺重力勘查規(guī)范下列術(shù)語和定義適用于本文件。3.1航空重力測量airbornegravitysurvey利用航空重力儀、定位儀等組成的航空重力測量系統(tǒng)，借助于專業(yè)測量飛機進行的空中重力測量，通過測量地球重力場（異常）研究基礎(chǔ)地質(zhì)及尋找礦產(chǎn)和能源資源等的一種地球物理方法。3.2內(nèi)符合精度RMSofrepeatlinenoiselevel相同飛行條件和狀態(tài)下，各重復(fù)線重力測量數(shù)據(jù)與它們的平均重力數(shù)據(jù)差值的均方差。3.3外符合精度RMSofdifferentialsbetweenflightlineandgrounddata航空重力與相同位置上延到飛行高度的地面重力差值的均方差。3.4低通濾波截止頻率cut-offfrequencyoflow-passfilter濾除航空重力測量數(shù)據(jù)中的高頻干擾噪聲獲得滿足精度要求的重力信息的低通濾波頻率。1DZ/T2019—XXXX3.5航空重力空間分辨率spatialresolutionofairbornegravity表示成以航空重力系統(tǒng)的測量精度可區(qū)分重力異常的最短波長，通常定義為可區(qū)分重力異常半波長（c/2）縮略語下列縮略語適用于本文件。3.6.1GNSS全球?qū)Ш叫l(wèi)星系統(tǒng)GlobalNavigationSatelliteSystem總則航空重力測量主要用于基礎(chǔ)地質(zhì)研究、油氣資源和固體礦產(chǎn)勘查。航空重力標量Δg測量，測量精度應(yīng)優(yōu)于1×10-5m/s2。20002000收集和分析測區(qū)及鄰區(qū)地質(zhì)、地球物理、地球化學(xué)、遙感、鉆孔及測繪等資料，收集測區(qū)的地形、地理、水文和氣象資料，收集測區(qū)內(nèi)的各種重力資料，并列表給出區(qū)內(nèi)不同巖性密度的統(tǒng)計表。設(shè)計書編制前，宜前往測區(qū)現(xiàn)場踏勘，了解飛行作業(yè)條件，包括：管制測量區(qū)域空域的相關(guān)部門、空域繁忙程度、作業(yè)機場、夜航條件、空中與地面保障條件、交通條件、氣候變化和生活保障等。22～33km～5km。1:10000010（1:50000～1:2000001:50000。1。3飛行高度確定應(yīng)考慮以下因素：航空重力測量應(yīng)綜合考慮測量的地質(zhì)效果、測區(qū)的地形條件和飛行時的氣流大小，以較低的同一高度平飛。表1 測線距測比尺照表測線間距/m測量比例尺測點距/m5001:5000015～3010001:10000030～6020001:20000030～601:50000200m～400m。200m。1:100000200～400m1:200000400～600m100m。表2 不同形件測比尺、行度照表測量比例尺平均離地飛行高度m平原地區(qū)（高差＜100m）丘陵地區(qū)（高差＜200m）低山區(qū)（高差＜400m）山區(qū)（高差＜600m）高山區(qū)（高差＞600m）1:50000200～400300～500400～600500～8001:100000300～400400～600500～700600～1000800～12001:200000400～600500～700600～800700～1200800～1400注1：應(yīng)以不高于表中的平均離地飛行高度開展航空重力測量420m。3180～250km/h160km/h。400km/h。導(dǎo)航定位應(yīng)采用全球?qū)Ш叫l(wèi)星系統(tǒng)，或采用同等以上定位精度的導(dǎo)航定位系統(tǒng)。80m5km±80.0m±0.1m。102h10m。80m5km5km33455表3 航空由間力常精度考

單位：×10-5m/s2自由空間重力異常總精度外界加速度干擾改正均方誤差正常重力場改正均方誤差重力高度改正均方誤差0.60.30.490.10.150.80.30.720.10.151.00.30.940.10.15表4 外界擾速改均誤差考

單位：×10-5m/s2外界加速度干擾改正均方誤差垂向加速度改正均方誤差水平加速度改正均方誤差厄特渥斯改正均方誤差0.490.300.200.300.150.610.450.200.300.200.720.550.300.300.200.820.650.300.350.200.940.750.350.400.20表5 航空格力常精參考

單位：×10-5m/s2布格重力異?？偩茸杂煽臻g重力異常總精度地形改正均方誤差（200m以遠）地形高程均方誤差(m)0.610.600.07±5~±100.810.800.11±10~±201.031.000.24±20~±30注：地形改正均方差引用DZ/T0171-1997---表6，取其最大誤差A(yù)。6（），5×10-5m/s2/d，6.3.2.4 （）（1×10-5m/s2（m）(m)0.02+1×10-6×D，D為基線長度，單位m。71m（m）0.1m。3m。1mYant（前為正）、Zant（上為正），測量精度應(yīng)優(yōu)于0.1m。導(dǎo)航顯示器應(yīng)當安裝在飛機前部方便駕駛員監(jiān)控和操作的位置。4h。200km10200m），212h2GPSGB/T183142h±0.1m8DZ/T0004地面重力基點引點前，應(yīng)完成用于重力引點重力儀的調(diào)節(jié)和校驗、格值標定和性能試驗。收集測區(qū)附近的2000國家重力基本網(wǎng)系統(tǒng)的重力基點資料。20000.05×10-5m/s2。200048h3h）24h2h（）9（）GNSS6（6PDOP2.5，1m。＞5km（3GNSS40minGNSS815min。15°。1040min。20×10-5m/s240km。252B。。1%。/113~53%~5%（）100%GNSS（含）（RMS）0.05m/s（RMS）GNSSPDOPGNSS6PDOP2.5（GNSS）12（（GNSS、GNSS4（GNSS100%。巖（）巖（）（）30300g3cm×4cm×4cm4cm×4cm×4cm50cm×50cm×50cm（13DZ/T0004（（。（（）地（礦（））GNSSGNSSPDOPD14DZ/T2019—XXXXGaussUTM（1）計算

g g0) g????（）b z z U E H F K 0式中：g 10-5m/s2；bg 10-5m/s2；bzg ，10-5m/s2；zzg0 10-5m/s2；zUv ，10-5m/s2；UEa 10-5m/s2；EHa ，10-5m/s2；HFa 10-5m/s2；FKa 10-5m/s2；Kg0 （IAG）1980年國際正常重力公式或CGS2000的正常重力公式等模型計算得到，10-5m/s2。、載體向速度平速度正aH由解算軟計求得。、（2）計算v2 v2aEsvE E N ???????????（）Nh Mh式中：aE

----厄特渥斯改正值，10-5m/s2；15，Rad；Rad；E、N----分別是飛機東向和北向速度，m/s；N、M 分為酉和午圈曲半，m；h m。厄特渥斯改正誤差通常按公式（3）計算

RAAA?（3）式中：daE 10-5m/s2；A ），Rad；R m；V ，m/s。（4）計算

h)107

h2???4）式中：aF 10-5m/s2；，Rad；h m；h m。高度改正誤差通常按公式（5）計算

617sh???????????5）式中：daF----高度改正誤差，10m/s；-5 2，Rad；dh ，m。（IAG）1980（6）0 g0 式中：160g 10-5m/s2；0Rad。dg0（7）

DZ/T2019—XXXX0

???????????（7）式中：0dg ，10-5m/s2；0，Rad；dd m。2（8）為：1g0)t

???????????（8）K Z

t1t0（9）(g1g0)(g1g0)t

?????????9）tK Z t

10b b t10式中：aK

---- 10-5m/s2；g0g1 10-5m/s2；Z Zt0t s；g0，g1

----分別為前校、后校基點處的重力場值，10-5m/s2。b b在數(shù)據(jù)進行布格改正和水平調(diào)平之前，通常需要對數(shù)據(jù)進行兩方面的整理。15km17DZ/T2019—XXXX（如1）1(x,yz0),,)是具有密度圖1 地形正意圖Z該點的引力位沿ZZR

R

x)y)

(x)2(y)2(z)2]00y)x)x)2y)2z)200(x)(y)

xy2z2

???（1）0z)arctan 02 2 2

2x1y(z)(x)(y)(z0)

1z1式中：

----重力地形改正值，m/s2；G（6.67×10-11m3/(kg×s2)；地形直棱柱體巖石剩余和水體虧空密度，kg/m3；（水準面之上陸地地形直棱柱體巖石剩余密度，一般選用2.67×103kg/m3；水準面之下水底地形直棱柱體水體虧空密度，海水虧空密度一般選用(1.02~2.67)×103kg/m3，淡水虧空密度一般選用(1.00~2.67)×103kg/m3）Rm；x1、x2、y1、y2、z1、z2分別是每個直棱柱體對應(yīng)的角點坐標。1:100000()Radar20～166.7km。局部地形網(wǎng)格數(shù)據(jù)要求更為精細，網(wǎng)格間距較18DZ/T2019—XXXX）（10）

????????????（1）式中：bg 10-5m/s2；bg 10-5m/s2；aT

----地形改正值，10-5m/s2。（a）計算控制線法調(diào)平后航空重力異常數(shù)據(jù)沿垂直測線方向的網(wǎng)格數(shù)據(jù)G1；b）G1G2G3G1–G2；G3c）將每條測線的噪聲值從控制線法調(diào)平后航空重力異常數(shù)據(jù)中去除就獲得了第一次微調(diào)平處理后的航空重力數(shù)據(jù)。19(a)-(e調(diào)平后的航空重力異常數(shù)據(jù)需進行數(shù)據(jù)噪聲處理。宜采用的主要步驟如下：a）計算微調(diào)平后航空重力異常數(shù)據(jù)沿垂直測線方向的網(wǎng)格數(shù)據(jù)G1；b）計算微調(diào)平后航空重力異常奇數(shù)和偶數(shù)號測線數(shù)據(jù)沿垂直測線方向的網(wǎng)格數(shù)據(jù)G2、G3；c）計算G4=G1-G2、G5=G1-G3分別獲得偶數(shù)和奇數(shù)號測線的噪聲網(wǎng)格數(shù)據(jù)；d）合并和平均奇數(shù)、偶數(shù)號測線的噪聲網(wǎng)格獲得所有測線的噪聲網(wǎng)格數(shù)據(jù)；將所有測線的噪聲網(wǎng)格數(shù)據(jù)從G1中去除后獲得噪聲去除后的航空重力異常數(shù)據(jù)。航空重力數(shù)據(jù)處理結(jié)束后，應(yīng)編制以下基礎(chǔ)圖件：航空自由空間重力異常剖面平面圖，航空自由空間重力異常等值線平面圖，航空布格重力異常剖面平面圖，航空布格重力異常等值線平面圖，航空重力測網(wǎng)分布圖（航跡圖），飛行高度剖面平面圖等。航空重力剖面平面圖、航空重力等值線平面圖的成圖比例尺應(yīng)等于或小于測量比例尺；當在一個測區(qū)內(nèi)采用兩種測量比例尺時，將小比例尺作為成圖比例尺。其它成圖比例尺應(yīng)小于測量比例尺。1mm1/4。DZ／T0069基礎(chǔ)圖件技術(shù)說明部分應(yīng)分別包含以下內(nèi)容：a）測量飛行日期、使用飛機型號、機場、測量比例尺和色標等；b）導(dǎo)航定位方法及儀器設(shè)備型號、實際達到的導(dǎo)航和定位精度；c）航空重力儀系統(tǒng)型號、靈敏度（或噪聲水平）、數(shù)據(jù)收錄方式及采樣點間距；d）航空重力數(shù)據(jù)主要校正內(nèi)容；e）高度計型號、測量飛行高度及其分布情況；f）成圖坐標及投影。20DZ/T2019—XXXXd）(10-5m/s2)（12）為：

i1 2nn2i??????i1 2nn2i式中：ii（i=1,2,3，?n）110-5m/s2；in 685表6 測線架測質(zhì)評等級表內(nèi) 容整架次單條測線評價等級均方差（10-5m/s2）σ≤0.8σ≤1.0一級0.8＜σ≤1.01.0＜σ≤1.2二級1.0＜σ≤1.21.2＜σ≤1.5三級σ＞1.2σ＞1.5不合格注：交叉點處測線與控制線飛行高度差大于30m的不參加統(tǒng)計）符合精度計算詳見附錄。GNSS（21飛行高度評價由平均高度

表示，也可統(tǒng)計全測區(qū)高度分組出現(xiàn)的頻數(shù)及所占分數(shù)。平均高度h（m）（13）為：式中：

NN1N i???????????????13）1Ni1△hi im；N W（14）：W式中：

n/N%

????????????????14）n N 按設(shè)計的偏航距要求，計算全測區(qū)測點偏離預(yù)定測線的測線連續(xù)長度（大于設(shè)計偏航距的線段長度），用以評價導(dǎo)航質(zhì)量。沿垂直于測線方向，計算飛行測量獲得的實際測線間距與設(shè)計書預(yù)定的測線間距的差值，計算這些差值的平均值即為全區(qū)的測網(wǎng)疏密度。測網(wǎng)疏密度是用來評價飛行測線平面上實際均勻分布情況的指標。測網(wǎng)疏密度計算方法（15）如下：1NND diNi1

????????????(15)式中：di----第i個位置處實際測線間距，m；N ----設(shè)計書預(yù)定的測線間距，m。按照一定的間隔沿測線方向進行取樣，獲得測線分布的均勻程度。GNSSGNSS22GNSSGNSS）（不Sp(m161N11N1Ni1(XiX)2式中：

?????????????(16)Xi----表示某固定點上第i次定位系統(tǒng)的坐標觀測值，m；X----諸Xi的算術(shù)平均值，m；N i =?N。）利用每架次（日）航空重力前后校的數(shù)據(jù)，計算出架次零漂量，并進行等級評定。2Zg(10-5m/s2（17）為：ZgZgZgZ

1g

???????????????（1）ZZg(10-5m/s2（18）ZZg g1g0)g1g0???????????（18）Z Z Z b b式中：g0g1 10-5m/s2；Z Zg0，g1

----分別為前校、后?；c處的重力場值，10-5m/s2。b b航空重力零漂gZ

(10-5m/s2)評價分級：一級二級 1.0三級 2.0不合格

gZgZgZgZ

≤1.0≤2.0≤3.0＞3.0Z如果架次零漂g大于3×10-5m/s2，為不合格架次，須進行重飛。Z23i1 i1 2nn2i

????????????????（1）式中：，10-5m/s2。i 第（i=1,2,3，…n）10-5m/s2；n a）項目簡介，數(shù)據(jù)處理使用的軟件程序及相應(yīng)版本；b）用于數(shù)據(jù)處理(包括預(yù)處理)的資料清單；c）獲取原始資料的測量飛行日期，使用的飛機型號、機場、測量比例尺和測線方向；d）航空重力數(shù)據(jù)編輯、各項改正、采用的濾波方法及參數(shù)、處理流程、圖件編制等；e）圖件繪制時使用的坐標、投影和轉(zhuǎn)換方法等；）8.3航空重力資料推斷解釋根據(jù)項目任務(wù)目標可分為區(qū)域性基礎(chǔ)地質(zhì)調(diào)查和油氣資源調(diào)查2種類型。24（、ΔgΔg（Δg）ΔgΔg梯度G—4（G—2009—0001）（25大致估計密度地質(zhì)體形態(tài)、產(chǎn)狀、空間位置和物性參數(shù)。EulerParkerΔgΔgAB、C、DEＡ類：基巖凸起異常，指重力基底的局部構(gòu)造高點。B類：中－基性侵入巖異常，一般指基性侵入巖的范圍。C類：中酸性侵入巖異常，指中酸侵入巖的范圍或局部高點。D類：火山巖異常，推斷與火山巖有關(guān)的局部異常。E類：與沉積層有關(guān)的異常。根據(jù)目標任務(wù)要求，編制航空重力推斷解釋圖件。26（E）a）b）d）e）結(jié)論與建議。a）c）a）c）a）航空自由空間和布格重力異常剖面平面圖；b）航空自由空間和布格重力異常等值線平面圖。27a）b）推斷基底巖相構(gòu)造圖；c）推斷區(qū)域地質(zhì)構(gòu)造圖；d）局部構(gòu)造異常及油氣遠景評價圖。（（）28附 錄 A（資料性附錄）航空重力測量工作設(shè)計編寫提綱前言簡述項目來源、任務(wù)、選區(qū)依據(jù)、工作時間、提交成果等，以及其他有關(guān)要求。簡述野外作業(yè)、室內(nèi)編圖、基礎(chǔ)圖件草圖驗收和成果報告編寫時間等安排計劃。（。（礦（簡述出開工前系統(tǒng)和各項測試準備工作，包括測量系統(tǒng)在室內(nèi)、室外的測試方法和結(jié)果要求。29簡述測區(qū)氣象、交通、地形等自然概況（附測區(qū)交通圖、地形圖、降雨分布圖等）。/簡述航空重力測量使用的導(dǎo)航定位系統(tǒng)、方法及其質(zhì)量要求。簡述航空重力動態(tài)精度和質(zhì)量檢查測試方法，控制線飛行時航空重力系統(tǒng)的儀器狀況及氣象要求，重復(fù)線測量方法與要求，不同高度或加密測量方法與要求。簡述資料編錄整理方法，各類測量數(shù)據(jù)文件命名方法等。簡述測量總精度使用的統(tǒng)計方法和要求達到的測量精度。簡述數(shù)據(jù)預(yù)處理方法、流程、質(zhì)量控制方法，重測、補測線段確定。簡述處理成圖使用的軟件系統(tǒng)、各項改正、基礎(chǔ)圖件編制、提交的圖件。30簡述航空物探解釋的基本思路，擬采用的數(shù)據(jù)轉(zhuǎn)換處理方法，解釋的技術(shù)路線等。提交階段性和最終成果的名稱、圖件種類及驗收時間。簡述項目組人員組成、職責及分工；進度安排與效率、技術(shù)裝備等。按有關(guān)預(yù)算規(guī)定執(zhí)行。完成任務(wù)計劃的預(yù)算說明和相應(yīng)的預(yù)算表格。組織措施、質(zhì)量保證、生產(chǎn)安全措施等。31附 錄 B（規(guī)范性附錄）航空重力測量內(nèi)外符合精度評價方法i1 ni1 nijn2

， （j,,,m）???????B）

mmj1m

???????????????B.）ij——水平調(diào)整后第j條重復(fù)線公共段各點重力觀測值Fij與該點各重復(fù)線重力觀測的平均值Fi之差，10-5m/s2；ijFjF)，（in；jm）Fi

mmj1 ，（in）mn每條重復(fù)線重力數(shù)據(jù)的水平均值（10-5m/s2）：n Fi1 ，（jm）j nn所有重復(fù)線重力平均場數(shù)據(jù)的水平均值（10-5m/s2）：n Fi1 n32式中：m——重復(fù)線的數(shù)目；n——重復(fù)線公共段數(shù)據(jù)點數(shù)。（10-5m/s2）m n( m n( )2ijj1i1mn

?????????????（B3）（ij（或）（Pi）10-5/s2 ij jm n [(FFP)P] ij jm n [(FFP)P]2ij1i1mn外符合標準場的重力水平均值P（10-5m/s2）：nniPi1 ?????????????（B5）n式中：ijF——ji10-5m/s2；ijiP——在第i個測點的外符合實際標準重力值，10-5m/s2。im——參加統(tǒng)計的測線數(shù)；n——參加統(tǒng)計的測點個數(shù)。33附 錄 C（資料性附錄）航空重力儀檢測結(jié)果記錄表航空重力儀檢測結(jié)果記錄表，見表C.1。表C.1航空重力儀檢測結(jié)果記錄表測區(qū)稱： 儀型： 測地：靜態(tài)定位精度（m）參