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

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