區(qū)域重力調(diào)查規(guī)范

中華人民共和國地質(zhì)礦產(chǎn)行業(yè)標準代替DZ/T0082—93IDZ/T0082—2021前言 Ⅲ 12規(guī)范性引用文件 1 13.1術語和定義 1 2 24.1應用條件 24.2應用范圍 34.3基本要求 35技術設計 35.1資料收集和分析 35.2野外踏勘 4 45.4總精度及誤差分配 45.5基點網(wǎng)建立 55.6測地工作 6 75.8物性工作 7 85.10設計書編寫 85.11設計書審查與變更 86重力儀 8 86.2重力儀格值標定 8 7.3重力點布設與觀測 7.5控制點、測點平面位置和高程的測量 7.10原始資料質(zhì)量評價及資料驗收 ⅡDZ/T0082—2021 8.1一般要求 20 24 8.10自由空間重力異常及精度計算 8.11布格重力異常及精度計算 8.12均衡重力異常計算 8.13密度統(tǒng)計與精度計算 9基礎圖件編制 9.1圖件種類及一般要求 10數(shù)據(jù)處理和解釋推斷 10.1一般要求 10.4定量反演 10.5綜合解釋推斷 10.6解釋推斷成果的表達 11.1一般要求 11.2成果報告編寫 11.3成果提交 附錄A(資料性)重力基本網(wǎng)和重力儀格值標定場 附錄B(資料性)重力儀檢查與調(diào)節(jié) 41附錄C(規(guī)范性)重力儀觀測值計算方法 附錄D(規(guī)范性)密度測定方法 附錄E(資料性)地形改正和均衡改正 附錄F(規(guī)范性)表格樣式 附錄G(規(guī)范性)兩點間多臺重力儀多次觀測結果的精度統(tǒng)計方法 附錄H(規(guī)范性)重力固體潮改正 附錄I(規(guī)范性)重力基點網(wǎng)平差 附錄K(資料性)主要地質(zhì)目標的識別標志 附錄L(資料性)區(qū)域重力調(diào)查成果報告編寫提綱 Ⅲ本文件按照GB/T1.1—2020《標準化工作導則第1部分：標準化文件的結構和起草規(guī)則》的規(guī)定———1:200000測量比例尺改為1:250000比例尺(見3.1.1);——增加了兩個術語(見3.1.9、3.1.10)和一個縮略語(——修改了五統(tǒng)一基本要求(見4.3);——增加了資料收集和分析(見5.1)及野外踏勘(見5.2)—提高了布格重力異常總精度(見5.4);——修改了基點網(wǎng)建立的技術指標(見5.5.2); ——修改了近區(qū)地形改正方法(見5.7.4); ——增加了利用DEM高程數(shù)據(jù)進行地形改正的技術精度評價(見8.6.2);——修改了重力正常場改正公式(見8.7.1); 本文件由中國地質(zhì)科學院地球物理地球化學勘查研究所起草?！?993年首次發(fā)布為DZ/T0082—93;1本文件適用于基礎地質(zhì)、礦產(chǎn)地質(zhì)工作中的陸地及淺海、灘涂1:100000和1:250000的區(qū)域重力GB/T14499—1993地球物理勘查技術符號DZ/T0004—2015重力調(diào)查技術規(guī)范(1:50000)DZ/T0069—1993地球物理勘查圖圖式圖例及用色標準DZ/T0153—2014物化探工程測量規(guī)范以研究區(qū)域重力場區(qū)域地質(zhì)為主要目的的重力勘查，常用比例尺為1:50000、1:100000、1:2相對于1號點的增量值和1號點相對于2號點的增量值。三程循環(huán)觀測法中一個閉合段內(nèi)的兩個非獨立增量的平均值或雙程往返觀測法求得的重力增4總則31:100000、1:250000區(qū)域重力調(diào)查，主要用于：4.3.1統(tǒng)一采用2000國家重力基本網(wǎng)系統(tǒng)。4.3.2統(tǒng)一采用2000國家大地坐標系和大地高程基準。4.3.3統(tǒng)一采用2000橢球面正常重力公式計算正常重力值。5.1.1.1.1應收集國家測繪部門最新出版的1:50000和1:25000的地形圖、三角點、GNSS控制點、5.1.1.1.2宜采用1:50000的地形圖作為野外數(shù)據(jù)采集手圖。5.1.1.2.4收集調(diào)查區(qū)內(nèi)已有1:250000、1:200000、1:50000區(qū)域地球化學調(diào)查基礎數(shù)據(jù)和成果45.1.1.2.5對可利用的成果圖件進行統(tǒng)一標準處理后數(shù)字化，配準統(tǒng)一坐標系統(tǒng)，套合到地理底5.3.1本文件涉及區(qū)域重力調(diào)查的基本比例尺為1:100000和1:250000。1:250000的調(diào)查以基區(qū)面積的10%。比例尺自由網(wǎng)測點密度(每個測點所控制的面積)5比例尺常均方誤差10-?m/s測點重力均方誤差點位影響士0.11士0.001士0.15士0.28士0.15士0.28士0.001單位為10?m/s2(建二級基點網(wǎng)時)測點重力觀測士0.11士0.05士0.056一級基點網(wǎng)二級基點網(wǎng)1.傳遞重力值2.控制二級基點網(wǎng)的重力聯(lián)測1.傳遞重力值1.按省(市、區(qū))統(tǒng)一設計，一般一次建成點網(wǎng)的需要3.至少包含(或聯(lián)測于)一個國家重力控制不超過12,并使起算用的國家重力控制點位于網(wǎng)的中部工安排情況分區(qū)建立，充分利用一級基點網(wǎng)進行控制1.用于聯(lián)測的重力儀按三程循環(huán)觀測或雙程往返觀測法時，均方誤差E一般不大于士50×10-?m/s2多一倍的堅強邊相連接1.用于聯(lián)測的重力儀按三程循環(huán)觀測或雙程往返觀測法時，均方誤差eo一般不大于士60×10-8m/s21.網(wǎng)的精度ewa不低于士50×10-m/s22.聯(lián)測精度eb不低于±30×10-?m/s21.網(wǎng)的精度ewa不低于±100×10-?m/s2.聯(lián)測精度e不低于±40×10-?m/s2比例尺mm士0.55.6.2各級基點的平面位置和高程采用GNSS測定或從1:50000、1:100000地形圖上量取。7比例尺近區(qū)中區(qū)遠一區(qū)遠二區(qū)士0.070士0.28士0.213a)中區(qū)地形改正宜使用國家測繪地理信息局1:50000DEM高程數(shù)據(jù)或用數(shù)字化1:50000地896.2.2.1建立省(市、區(qū))級重力儀格值標定場時，所用重力儀格值測定的相對均方誤差應小于級基點時，兩段格值的相對變化應不大于1/3600;進行測點觀測時DZ/T0082—20216.3重力儀性能試驗6.3.1重力儀性能試驗內(nèi)容6.3.2靜態(tài)試驗6.3.3動態(tài)試驗6.3.3.1試驗時間應不少于12h,選擇兩個或兩個以上的試驗點，試驗點間重力差不小于3×10-?m/s2,兩點間單程觀測時間間隔一般不大于20min。6.3.3.2對重力儀的動態(tài)觀測結果應進行固體潮改正。經(jīng)改正后繪制動態(tài)零點位移曲線，統(tǒng)計儀器的動態(tài)觀測精度，其觀測均方誤差應小于設計的測點重力觀測均方誤差(eg)的1/2或設計的基點聯(lián)測均方誤差ε,資料整理方法詳見附錄C。動態(tài)觀測均方誤差計算公式為e——動態(tài)觀測均方誤差，單位為10-?m/s2;m——增量的總個數(shù)；6.3.3.3重力測點觀測閉合于基點的時間應在動態(tài)零點位移曲線的線性時段內(nèi)。動態(tài)零點位移曲線與n——獨立增量數(shù)。6.3.4一致性試驗6.3.4.1一致性試驗點一般不少于15個，點間重力值的變化應足夠大(一般不小于3×10-?m/s2),采6.3.4.2各臺儀器間一致性均方誤差計算公式為a;——第j臺儀器在i點上的觀測值，單位為10-?m/s2;m——參加一致性試驗觀測的儀器臺數(shù)；6.3.5調(diào)節(jié)測程后讀數(shù)穩(wěn)定時間試驗7野外工作7.1儀器準備40×10-?m/s2,讀數(shù)精度小于或等于10×10-?m/s2的要求。7.1.1.2重力儀工作前應進行儀器調(diào)節(jié)、格值標定、靜態(tài)試驗、動態(tài)試驗以及多臺儀器間的一致性7.1.1.3標稱精度低于20×10-?m格值的變化情況。標稱精度高于20×10-?m/s2的相對重力儀一般每三年檢查一次儀器的格值或格值7.2.3.3閉合差不超過設計的平均增量的均方誤差eb的2√N倍河北02,內(nèi)蒙古03,山西04,遼寧05,吉林06,黑龍江07,上海08,江蘇09,浙江10,安徽11,江西12,福建13,山東14,廣東15,廣西16,湖北17,湖南18,河南19,四川20,云南21,貴州22,西藏23,陜西24,甘肅25,青海26,新疆27,寧夏28,海南29,重慶30。體位于縱坐標為4541000和4542000、橫坐標為20544000和20545000的四條線所限定的1km2的顆個m高程收斂精度mCORS服務無限制注：收斂精度為固定解的要求，浮點解可放寬到兩倍執(zhí)3%～5%,檢查點應不少于30個。點數(shù)不得超過檢查點數(shù)的1%,否則應擴大檢查工作量7.5.3.1對于野外實測取得的測點高程，應與比例尺為1:50000或1:100000的地形圖上高程進行7.8.3剖面質(zhì)量檢查工作量不少于5%,且檢查總點數(shù)不少于30個。c)測點重力觀測的記錄(自動記錄須打印輸出);7.9.7記錄的數(shù)字均采用阿拉伯數(shù)字。記錄的有效數(shù)字應和所用儀器的讀數(shù)精度以及測量的方法相適7.10.1.2分類標準如下。1)野外原始圖件；c)室內(nèi)計算的全部內(nèi)容經(jīng)過100%復算，正確無誤。8.1.2計算精度b)格值及校正系數(shù)計算準確至10-6。8.1.2.2基點網(wǎng)聯(lián)測中，重力增量的單項計算和各項改正均取準至0.001×10-?m/s2?；c重力值取準至0.001×10-?m/s2。8.1.2.3測點重力值觀測的單項計算和各項改正均取準至0.001×10-?m/s2,測點重力值取準至8.1.2.4布格重力異常、自由空間重力異常計算中的各項取準至0.001×10-?m/s2,最后結果取準至8.2.1小測程重力儀格值計算見公式(5),格值測定的相對均方誤差計算見公式(6),不同測段格值的相對變化計算見公式(7),重力儀格值的相對變化(相對于使用格值)k重力儀格值；ms——重力儀讀格差平均值的均方誤差(為第i個讀格差與平均讀格差的差),ms=用平均讀格差的相對均方誤近似代表格值測定的相對均方誤為k的均方誤差)。當不忽略格值標定場兩點間重力差值測量相對均方誤差時，格值測定的相對均方誤差式中：C——重力儀格值表校正系數(shù)；△g;——第i臺儀器的第j個獨立增量，單位為10-?m/s2。mAg——標定場均方誤差，單位為10-?m/s2,計算詳見附錄G;m?——偶然誤差，單位為10-?m/s2;m?——第一類半系統(tǒng)誤差，單位為10-?m/s2;m?—第二類半系統(tǒng)誤差，單位為10-?m/s2。8.4基點網(wǎng)平差及精度計算8.4.1固體潮改正和零點位移改正基點重力聯(lián)測結果計算時應進行固體潮改正和零點漂移改正。計算方法詳見附錄H和附錄C。對8.4.2基點網(wǎng)平差8.4.2.1工作應分級進行。一次建成的一級基點網(wǎng)，應在國家重力基本網(wǎng)的控制下在全省(區(qū))范圍內(nèi)8.4.2.2基點網(wǎng)平差可采用條件平差法，以各邊段的獨立增量數(shù)為權進行。平差后應計算各基點重力8.4.3聯(lián)測精度N——基點網(wǎng)的聯(lián)測邊段數(shù)； (13) (14)δ;——基點網(wǎng)第j邊段上各第i個獨立增量與該邊段平均增量之差，單位為10-?m/s2。8.4.4基點重力值均方誤差8.4.4.1根據(jù)單位權均方誤差(μ)計算基點重力值均方誤差(eG),計算公式見式(15)和式(μ——單位權均方誤差，單位為10-?m/s2;EG——重力值均方誤差，單位為10-?m/s2;r——多余條件數(shù)(閉合圈數(shù));P;——第j邊段的權系數(shù)(其數(shù)值等于該邊段的獨立增量個數(shù));8.4.4.2根據(jù)支基點的重力聯(lián)測均方誤差(cm)計算支基點重力值均方誤差(e?),計算公式見式(17)和 (18)n——檢查點總數(shù)。8.6.1.2中區(qū)地形改正采用計算機進行，宜使用國家測繪地理信息局1:50000DEM高程數(shù)據(jù)或用數(shù)字化儀將1:50000地形圖數(shù)字化形成50m×50m或100m×100m節(jié)點高程，使用共用點法(或直接算方法參見附錄E。8.6.2精度估算8.6.2.1近區(qū)地形改正精度計算公式如下：δ;——檢查值與原計算值之差，單位為10?m/s2;n——參與統(tǒng)計的檢查點數(shù)。 (22) (23)8.6.2.2中區(qū)地形改正采用節(jié)點高程計算時，采用同精度評價方法進行精度評價?？刹捎靡韵路椒ㄓ媋)當采用人工矢量化形成的DEM計算中區(qū)地形改正值時，用重新數(shù)字化后的DEM進行精度評b)當采用從測繪部門收集的DEM計算中區(qū)地形改正值時，用收集DEM數(shù)據(jù)附帶的高程中誤差進行中區(qū)地形改正誤差評價。檢查率不低于3%。誤差計算見公式(22)。8.6.2.3遠一區(qū)(2km～20km)地形改正誤差可采用分不同景觀區(qū)，通過試驗的方法確定。當采用原地質(zhì)礦產(chǎn)部區(qū)域重力數(shù)據(jù)庫進行地形改正時，可將0.114×10-?m/s2作為遠區(qū)地形改正誤差。參見附錄E。8.6.2.4遠二區(qū)(20km～166.7km)地形改正方法參見附錄E。當?shù)匦胃恼嬎悴捎迷刭|(zhì)礦產(chǎn)部區(qū)8.6.2.5誤差可采用分不同景觀區(qū)，通過試驗的方法確定。當采用原地質(zhì)礦產(chǎn)部區(qū)域重力數(shù)據(jù)庫進行遠區(qū)地形改正時，可采用表2中列出的遠區(qū)地形改正精度作為遠區(qū)地形改正誤差。參見附錄E。εgr——地形改正總均方誤差，單位為10-?m/s2;Egr?——0m~50m地形改正均方誤差，單位為10-?m/s2;egr?——50m~2km地形改正均方誤差，單位為10-?m/s2;egrg——2km~20km地形改正均方誤差，單位為10-?m/s2;Egr——20km～166.7km地形改正均方誤差，單位為10-?m/s2。8.7正常重力值改正及誤差計算8.7.1將CGCS2000橢球參數(shù)代go=978032.53349(1+0.00530244sin2φ—0.000005δgh=[0.3086(1+0.0007cos2φ)—0.72×10-h]h其他符號含義見8.9.1.3。δgh=—0.3086(1+0.0007cos2φ)(h?-ho)其他符號含義見8.9.1.3。δgh=0.3086(1+0.0007cos2p)(h?-h?+h?)其他符號含義見8.9.1.3。8.8.2高度改正均方誤差計算egh=±[0.3086(1+0.0007cos2φ)-1.44×10-7h]eh其他符號含義見8.9.2.3。Egh=±0.3086(1+0.0007cos2φ)ehO其他符號含義見8.9.2.3。8.9.1中間層改正其他符號含義見8.9.1.3。其他符號含義見8.9.1.3。h——重力點大地高程(測點低于高程基準面時為負值),單位為米(m);ho——重力觀測時從大地高起算的瞬時潮高(其數(shù)值一般在0m～5m范圍內(nèi)變化，瞬時海面低于h?——海面觀測時重力儀至海底的高度(始終為正值),單位為米(m);h?——重力觀測時的瞬時水深(始終為正值),單位為米(m);其他符號含義見8.9.2.3。其他符號含義見8.9.2.3。a在淺海域重力儀至海底的高度(h?)和觀測時刻的瞬時水深(h?)的測量誤差很小，其影響可忽略不8.9.3布格改正及精度計算△gF=g-go+δgh△gp=g-g?+δg+0.0419p?(2h?-h?)g——測點重力值，單位為10-?m/s2;po——海水密度(1.03g/cm3);EA=±√e3+e3+eE△F=±√e3c+e2+e+(0.0419poEo)2Eg——點位誤差引起的正常重力值均方誤差，單位為10-?m/s2。8.11.1布格重力異常計算公式為△gs——布格重力異常值，單位為10-?go——正常重力值，單位為10-?m/s2;g——測點重力值，單位為10-?m/s2;δgh——高度改正值，單位為10-?m/s2;ggr——地形改正值，單位為10-?m/s2。8.11.2布格重力異常均方誤差計算公式為eg——測點重力值均方誤差，單位為10-?m/s2;Egr——地形改正均方誤差，單位為10-?m/s2;8.12均衡重力異常計算8.12.1一般采用艾里一海斯卡寧均衡模式進行均衡改正。改正方法參見附錄E。8.12.2均衡重力異常計算公式為△g?=△gs+δgiδgi——均衡改正值，單位為10-?m/s2。8.13.1密度統(tǒng)計8.13.1.3當同類巖石密度標本的數(shù)量超過30塊時：8.13.1.4在按巖性統(tǒng)計密度的基礎上，以厚度為權計算地層單元(一般為系或統(tǒng)或同比例尺區(qū)域地質(zhì)8.13.2精度計算δ;——各點上第i次測定值(包括該點參與計算平均值的原始測定值和所有檢查測定值)與該9基礎圖件編制9.1.2一般要求9.1.2.2繪圖儀的分辨率要求在600dpi以上；繪圖儀精度要求最大累積誤差不超過士0.a)點位數(shù)據(jù)圖的主要內(nèi)容有重力測點點位、點號、布格或號。比例尺為1:250000、1:100000的點位數(shù)據(jù)圖上表示國家重力基本點、國家重力一等點b)各級物探重力基點的編號，應按7.2.4.3要求執(zhí)行。d)各級重力基點用不同的符號表示見表8。當基點和測點重合時，保留基點符號及編號，省略測符號符號大小★國家重力基本點紅[RGB(255,0,0),CMYK(0,99,100,0)]國家重力一等點紅[RGB(255,0,0),CMYK(0,99,100,0)]黑[RGB(35,31,32),CMYK(0,0,0,100)]△物探重力Ⅱ級基點黑[RGB(35,31,32),CMYK(0,0,0,100)]9.1.3.2布格重力異常和自由空間重力異常平面圖。搜索半徑5km;比例尺為1:250000的圖件，網(wǎng)格距1000m×1000m,搜索半徑10km)。距要求見表9。比例尺布格重力自由空間重力異常注記大小(46,49,146),CMYK(100,100,0,0)]。零等值線和注記及重力高、重力低符號比例尺為1:250000的圖框采用高斯投影。詳見附錄J。——比例尺1:100000,1:2500009.2.2.1比例尺為1:100000、比例尺為1:250000的重力異常平面圖應繪內(nèi)外圖框線，分度帶標注9.2.2.2比例尺為1:100000、比例尺為1:250000的圖框計算出理論經(jīng)度差、緯度差的格網(wǎng)點和圖比例尺理論經(jīng)度差理論緯度差網(wǎng)點經(jīng)緯線寬5黑色[RGB(35,31,32)CMYK(0,0,0,100)]黑色[RGB(35,31,32)CMYK(0,0,0,100)]57網(wǎng)是我國測繪部門1957年在蘇聯(lián)航空重力隊協(xié)助下建立的。該網(wǎng)由27個基本點及80個I等點差分別為±0.3×10-?m/s2和士0.4×10-?m/s2,無絕對重力點控制，精度較低，不少點的標記已被使用。85網(wǎng)至57網(wǎng)之間，在我國存在系統(tǒng)差，其值為13.5×10-?m/s2左右(57網(wǎng)至85網(wǎng))。震局協(xié)作完成。該網(wǎng)由21個基準點、126個基本點和112個簡稱引點組成?；鶞庶c10-?m/s2;基本點平均中誤差為±6.6×10-?m/s2;引點平均中誤差為±8.7×10-?m/s2;網(wǎng)的平均中誤差為±7.4×10-?m/s2。2000網(wǎng)與85網(wǎng)之間無系統(tǒng)差，僅僅是2000網(wǎng)精度較85網(wǎng)略高。A.3.2標定場的重力聯(lián)測聯(lián)測工作至少應采用三臺標稱精度高于±20×10-?m/s2的相對重力儀進行。并按6.2.1、6.2.2、6.2.3的有關要求在國家重力長基線或國家級重力儀格值標定場上標定儀器的格值或格值表校正系數(shù)。臺儀器取得的獨立增量數(shù)不少于10個。(資料性)B.1不同類型重力儀檢查與調(diào)節(jié)ZSM型、CG-3型重力儀的光線位移靈敏度應調(diào)節(jié)至1.6×10-?m/s2~2.0×10-?m/s2之間，Worden重力儀的光線位移靈敏度應調(diào)節(jié)至2.0×10-?m/s2~5.0×10-?m/s2之間。B.1.1.3水準器的檢查與調(diào)節(jié)采用測水泡曲線的方法來檢查水準器是否調(diào)節(jié)正確。要求重力儀水泡曲線的極值點偏離正確位置(即水泡居中時的位置)不超過一小格(即圓周的B.1.2.1光線位移靈敏度的測定與調(diào)節(jié)位移靈敏度應調(diào)節(jié)至：改變讀數(shù)一格(約相當于重力變化1×10-?m/s2)目鏡筒內(nèi)黑影線移動9個~11個目鏡分劃。B.1.2.2正確讀數(shù)線的檢查與調(diào)節(jié)B.1.2.3橫水準器的檢查與調(diào)節(jié)B.1.2.4電子讀數(shù)零位和檢流計零位的檢查與調(diào)節(jié)B.1.2.5電子靈敏度的調(diào)節(jié)d)訪問SERVICE菜單，選擇XYOFFSET改正。按F5<STf)再次訪問LEVELING界面：若不記錄數(shù)據(jù)，按下MEASURE鍵兩次；若記錄數(shù)據(jù)，按<FINALDATA>以在AUTOGRAVFINALDATA界面下查看數(shù)據(jù)。系統(tǒng)會彈出一個b)訪問SERVICE菜單。選擇UserCalibration選項。按下ENTERE鍵，選中XYOFFSETc)按下<MEASURE/CLR>鍵。按下F5<LF5<FINALDATA>可在AUTOGRAVFINALDATA界面下查看數(shù)據(jù)。系統(tǒng)會彈出一個F5<FINALDATA>可在AUTOGRAVFINALDATA界面下查看數(shù)據(jù)單并讀數(shù)。在橫梁到達最低位置時，在后面的附表(CALIBRATETHELEVELCORRECTION數(shù)的影響是由它的余弦分量引起。在刻度值0位置處，一個格值對重力的影響值大約是0.0004×方向。FBK-Corr中兩組讀數(shù)之差為0.017×10-?m/s2表示要做一個檢流計格值的校正。當差為果機械參考位置已被適當調(diào)節(jié)到了檢流計的0刻度位置，在檢流計0位置左右各10個格值(也即0.1rad)所測得的兩組觀測重力值之差為(FBK-Corr)-0.042×10-?m/s2。與0點或中心位置的刻度差9時，與0點或中心位置的重力差-0.034×10-?m/s2;d)進入UltraGrav軟件。h)當進行連續(xù)測量時，在主界面選擇“CONTINUOUS”按鍵后進入做上角為“EnterDial節(jié)前、后零點漂移均隨時間按線性規(guī)律變化。在比例尺為1:100000(或1:250000)的重力測量工作B.2.3.4各次測程調(diào)節(jié)重復B.2.3.2式中：型號測程范圍金屬425半自動自動15自動(規(guī)范性)測路線為重力儀從基點開始，經(jīng)過一系列重力測點觀測，最后閉合于基點。本方法主要用于LCR型、C.2單程觀測法C.2.1石英重力儀和不附有格值表的LCR型重力儀1)停頓之前重力值計算公式為2)停頓之后重力值計算公式為1)停頓之前重力值計算公式為DZ/T0082—20212)停頓之后重力值計算公式為K——重力儀格值，單位為10-?m/s2;值(單位為10-?m/s2);值(單位為10-?m/s2);(單位為10-?m/s2);位為10-?m/s2)、停頓時間[單位為秒(s)]、停頓結束時的固體潮與停頓開始時的固體潮改正值之差(單位為10-?m/s2)。C.2.2附有格值表的重力儀b)有較長時間的停頓，需進行靜掉格改正時，停頓前、后測點重力值計算公式見式(C.8)、式1)停頓之前重力值計算公式為2)停頓之后重力值計算公式為g觀=G?+C(g'-gí)+(R;-R?)b)有較長時間的停頓，需進行靜掉格改正時，停頓前、后測點重力值計算公式見式(C.11)、式1)停頓之前重力值計算公式為g觀=G?+C(g;-gí)+(R,—R?)一2)停頓之后重力值計算公式為gí、T?、R?——起始基點上的重力觀測值(單位為10-?m/s2)、觀測時間[單位為秒g2、T?、R?——閉合基點上的重力觀測值(單位為10-?m/s2)、觀測時間[單位為秒g'、T、R某測點上的重力觀測值(單位為10-?m/s2)、觀測時間[單位為秒(s)]、固體潮改正值(單位為10-?m/s2);位為10-?m/s2)、停頓時間[單位為秒(s)]、停頓結束時的固體潮與停頓開始時的固體潮改正值之差(單位為10-?m/s2)。C.3雙程往返觀測法K——重力儀格值，單位為10-?m/s2。點相對于i點的重力增量值，單位為10-?m/s2;.C.4三程循環(huán)觀測法…S?、T?、R?——在1號點的第一次重力觀測值(讀格)、觀測時間[單位為秒(s)]、固體潮改正值(單位為10-?m/s2);Si、Ti、R{——在1號點的第二次重力觀測值(讀格)、觀測時間[單位為秒(s)]、固體潮改正值(單位為10-?m/s2);S?、T?、R?——在2號點的第一次重力觀測值(讀格)、觀測時間[單位為秒(s)]、固體潮改正值(單位為10-?m/s2);S2、T2、R2——在2號點的第二次重力觀測值(讀格)、觀測時間[單位為秒(s)]、固體潮改正值位為10-?m/s2)。 (規(guī)范性)D.1密度測定方法種類使秤臂與水平線成α?角度);精確度可達±(0.01～0.02)g/cm3。設用天平稱得標本在空氣中的質(zhì)量為p?,稱得標本在水中的質(zhì)量為p?,則標本的密度可用公式D.4大樣法公式見式(D.3):D.5密度測定結果的整理D.5.1密度測定結果的整理原則同類巖石的密度測定值，通常服從算術正態(tài)分布規(guī)律。對同一類巖(礦)石標本(有相當數(shù)量和代表D.5.2密度的算術平均值N——標本總塊數(shù)；D.5.3密度的標準離差當同一類巖石標本數(shù)目大于30塊時，可用計算法來統(tǒng)計平均值σ和標準離差D。方法是先將密度值按相等間隔△o分組。分組數(shù)n與標本總塊數(shù)之間的關系，在對數(shù)坐標軸上按線性比例變化?？蓞⒖急鞤.1標本密度分組表標本數(shù)(N)塊分組數(shù)(n)組456789DZ/T0082—2021n——分組數(shù)；N——標本總塊數(shù)。(資料性)E.1.1錐形公式n——方位數(shù)；G——萬有引力常數(shù)[6.67×10-?cm3/(g·s2)];i——地形傾角，單位為度()。E.1.2斜面公式式中各字母含義見E.1.1。E.1.3臺階公式E.1.3.1單層臺階階公式計算臺階地形改正值。量取h和rm后可采用公式長堤內(nèi)外高度一致時E.1.3.4傾斜臺階E.1.4扇形公式△gr—采用扇形公式計算的地形改正值，單位為10-?m/s2;G萬有引力常數(shù)[6.67×10-8cm3/(g·s2)],計算時輸入6.67×E.2.1中區(qū)和遠一區(qū)地形改正方法E.2.1.1中區(qū)和遠一區(qū)地形改正的計E.2.1.3遠一區(qū)地形改正采用平面公式計算地形改正值，使用1km×1km高程節(jié)點網(wǎng)按狹義地形改E.2.2共用點法地形改正E.2.2.1計算方法E.2.2.2計算公式G引力常數(shù)[6.67×10??cm3/(g·s2)];l——積分格距，選用1000m(或所用節(jié)點高程的網(wǎng)格距);R——計算點到補角的距離(圖E.3中PB),單位為米(m),R=1.105r;DZ/T0082—2021Z;——節(jié)點高程，單位為米(m);W——權系數(shù)；用20km畫圓，凡是節(jié)點落在圓外的就不參加計算，相當于外接口為鋸齒狀，誤差在允許范圍E.2.2.6精度計算由公式(E.8)對h;求導得出精度計算公式(E.13):其他字母含義見E.2.2.2。E.2.3遠一區(qū)地形改正結果對比以錐頂為中心計算1km×1km節(jié)點網(wǎng)高程公式見式(E.15):(△g計-△g理)/以湖南省地礦局物探大隊提供的60個點手算地形改正值做對比。手算時環(huán)帶劃分為：1.0～1.5;1.5～2.0;2.0～3.0;3.0～5.0;5.0～7.0;7.0～10.0;對這60個點統(tǒng)計了2km～20km的對比計算結果，見表E.2。按本文件規(guī)定遠一區(qū)(2km～20km)地形改正均方誤差為士0.114×10-?m/s2。共用點法計算結果精度滿足要求。系統(tǒng)誤差E.3.1球面圓域地形改正采用球面圓域模型計算遠二區(qū)地形改正值。見球體理論模型示意圖E.4,以地球球心為坐標原點，OZ軸通過測點P,XOZ面沿著點P所在的地理子午面，OX軸與地球旋轉(zhuǎn)軸ON分別在0Z軸的兩側。圖E.4球體理論模型示意圖R——地球平均半徑，取6371.025km;對式(E.16)積分得式(E.17):(R+Z)(cos20—cosθ)×1n{(2y-2(R+Z)cosθ+2[y塊的影響值為一個環(huán)帶的1/m,再令E(θ,r)=(2-r2-rcosθ-3cos20)(1+r2-2r3(cosθ—cos30)1n[r—cosE.3.1.2球殼質(zhì)塊的高程E.5球面坐標關系圖λ——某一球殼質(zhì)塊中心點的地理緯度(縱坐標),單位為度()。平均值作為該扇形塊的平均高程H。16273849858遠二區(qū)地形改正(20km～166.7km)%手算資料采用西安地區(qū)五點。手算方法是先把5'×5'高程數(shù)據(jù)展在1:1000000地形圖上，按表點12345XYE.4均衡改正E.4.1均衡改正模式實的地殼調(diào)整成為具有平均密度p=2.67g/cm3和正常地殼厚度T=30km的均勻地殼。E.4.2球面均衡改正地殼柱體體積V在球極坐標下積分表達式見式(E.24):即 表格樣式標定場名最大重力差點號經(jīng)度(或Y)(或X)點間重力增量123建場前所用標定場名年月日省(區(qū)、市)縣(區(qū))標石平面圖及照片點位略圖(應說明周圍100m建筑物分布情況)繪制日期填表日期表F.2重力基點檔案點號等級點名經(jīng)度(或Y)1:10萬圖幅緯度(或X)高程(Z)高程來源基點聯(lián)測時間受托管單位及點位所在地省(區(qū))縣(市)鄉(xiāng)村基點聯(lián)測路線標石平面圖及照片點位略圖繪制日期填表日期表F.3重力觀測點觀測記錄表時分123平均零點改正值表F.4GNSS觀測記錄表計算者：日驗者：重力圖幅號時間時間點位描述表F.5重力測區(qū)索引表12工作比例尺和測區(qū)面積3測區(qū)角點坐標角點順序經(jīng)度測區(qū)范圍示意圖1234567894中間層校正及地形改正所用密度5點位中誤差m6高程中誤差m78m9各級基點網(wǎng)均方誤差工作起止年月表F.6統(tǒng)一點號(或Y)(Z)m實測重力值(g)自由空間布格重力近區(qū)中區(qū)遠一區(qū)遠二區(qū)DZ/T0082—2nk123jk回觀測平均值123△n回觀測平均值平均值與總平均值之差…………nk12jk和1020δ0nδo00表G.3vi;(單位為10-?m/s2)計算表nkl2jk和10nk12jk和20V0n0和00000按式(G.12)計算出平均重力增量△go的均方誤差mA?0:(規(guī)范性)δ=—4.83+15.73sin2g′-1.59sin?g′F(φ)=0.998327+0.00167cos2φH.2.1求出計算時刻的儒略世紀數(shù)TT?-2415020.0——自1899年12月31日(世界時)12時起算至某一天實際天數(shù)；T——計算日的儒略世紀數(shù)(自1899年12月31日12時起算至某一時間t的儒略H.2.2計算六個天文引數(shù)s=270.43659°+481267.89057T+0.00198h=279.69668°+36000.76892T+0.0003p=334.32956°+4069.03403°T—0.01032°T2—0.0000N=259.18328°—1934.14201T+0.00208T2+0.000002T3ps=281.22083°+1.71902°T+0.00045°T2+0.00000e=23.45229°—0.01301T—0.000002°T2s——月亮的平黃經(jīng)，單位為度();h——太陽的平黃經(jīng)，單位為度();H.3求月亮的Cm/rm及cosZm0.0082cos2(s—h)+0.0006cos(2s—3h+ps)+0.0009coλ=s—0.0032sin(h—ps)一0.001sin(2h—2p)+0.001sin(s—3h+p+ps)0.0222sin(s—2h+p)+0.0007sin(s—h—p+ps)—0.0006sin(s—h)十0.1098sin(s—p)一0.0005sin(s+h—p-ps)+0.0008sin(2s—3h+pg)十0.0115sin(2s-2h)+0.0037sin(2s-2p)—0.0020sin(2s-2N)十0.0009sin(3s—2h—p)β=—0.0048sin(P-N)—0.0008sin(2h—p-N)+0.00.0895sin(s一N)+0.001sin(2s—2h+p-N)+0.0049sin(2s—p-N)十cosZm=sinp'sino+cosφ'cosocosHθ=(t—δt)×15°+h+L-180°(8?=8h,為世界時和北京時之差)φ′=φ-0°.193296sin2φλ——月亮的黃經(jīng)，單位為度();0——地方恒星時，單位為度();H.4求太陽λg=h+0.0335sin(h—ps)+0.0004sin(2h—2ps)………………(H.16)cosZ?=sing'sinesinλ,+cosg'(cosλscosθ+sinθcoseH.5計算固體潮改正值例如：B=31°20′,L=93°,2003年5月6日，19:45分。T?-2415020.0=37745.5為自1899年12月31日(世界時)12時起算至計算日(2003年5月6日)(規(guī)范性)重力基點網(wǎng)平差I.1平差計算物探重力基點網(wǎng)聯(lián)測結果的平差，采用條件平差方法較為適宜。下面結合實例介紹平差計算過程。I.2基點重力值計算I.2.1繪制基點網(wǎng)分布示意圖某重力一級基點網(wǎng)由八條邊組成三個閉(附)合圈，圖I.1為該基點網(wǎng)(一級)平差示意圖。圖1.1基點網(wǎng)(一級)平差示意圖在各邊段上標出平均重力增量值及變化方向(箭頭指向重力值減小的方向。其中G?~G?是重力一級基點，兩個畫五角星的點是國家重力基本點或I等點)。求出各閉(附)合圈的閉合差(圈內(nèi)箭頭為順時針方向各邊平均重力增量值之和減去箭頭為逆時針方向各邊平均重力增量值之和),并標在閉(附)合圈中央。把權(P等于各邊獨立增量數(shù))標在各邊段上。I.2.2列出改正數(shù)條件方程式-1,其余的a;=0;…解本例方程式(I.4),求得聯(lián)系數(shù)為abCSP1112113114l005100610071008總和232I.2.4計算改正數(shù)v;的值v?=7×10-?m/s2;G;=G+f?x?+f?x?+…+f,xnG?=G+x?+x?+x?=0×10-?+36.994×10-?+86.292×10-?+13.592×10-?=136.G?=G+x?+x?=0×10-?+86.292×10-?+13.592×10-?=99.8G?=G+x?=0×10-?+13.592×10-?=13.59G?=G+x?=0×10-?+32.970×10-?=32.97G?=G+x?+x?=0×10-?+32.970×10-?+47.252×10-?=80.2求得該基點的轉(zhuǎn)換系數(shù)qa=-0.8537,qb=-0.3659,q