地球物理反演理論、算法及應(yīng)用一書(shū)首先闡述了地球物理反演理論,特別是從哲學(xué)的角度、信息論的角度、數(shù)學(xué)的角度和地質(zhì)的角度,闡述了正則化反演思想的本質(zhì),通俗簡(jiǎn)明地介紹了常見(jiàn)非線性最優(yōu)化問(wèn)題的基本思想和算法。接著,本書(shū)闡述了重磁反演、電法反演、地震層析成像等單一反演方法的算法和程序?qū)崿F(xiàn),闡述了重磁電交叉梯度聯(lián)合反演的算法和程序?qū)崿F(xiàn),反演有二維反演,也有三維反演。最后給出了一些反演應(yīng)用實(shí)例。有些算法和程序來(lái)源于國(guó)外開(kāi)源代碼,更多算法和代碼來(lái)自于本課題組的獨(dú)立研究成果。
更多科學(xué)出版社服務(wù),請(qǐng)掃碼獲取。
目錄
前言
基礎(chǔ)篇
第1章 最優(yōu)化理論與方法概述 3
1.1 數(shù)學(xué)基礎(chǔ)知識(shí) 3
1.1.1 向量與矩陣 3
1.1.2 多元函數(shù) 7
1.2 凸集與凸函數(shù) 11
1.2.1 凸集 11
1.2.2 凸集分離定理 16
1.2.3 凸函數(shù) 23
1.3 最優(yōu)化問(wèn)題 28
1.3.1 數(shù)學(xué)模型 28
1.3.2 最優(yōu)解 29
1.3.3 最優(yōu)性條件 30
1.4 最優(yōu)化算法結(jié)構(gòu) 39
1.4.1 迭代算法的基本格式 39
1.4.2 收斂性與收斂速度 40
1.4.3 終止準(zhǔn)則 42
1.5 線性搜索 44
1.5.1 精確線性搜索算法 46
1.5.2 非精確線性搜索算法 74
1.6 最優(yōu)化方法 79
1.6.1 最速下降法 79
1.6.2 牛頓法 96
1.6.3 共軛梯度法 105
1.6.4 擬牛頓法 135
1.6.5 信賴(lài)域法 164
1.7 非線性最小二乘問(wèn)題 174
1.7.1 高斯-牛頓(Gauss-Newton)法 175
1.7.2 利文貝格-馬夸特(Levenberg-Marquardt)法 178
1.7.3 擬牛頓法 185
第2章 線性正則化反演理論 187
2.1 非線性問(wèn)題的線性化 187
2.2 線性反演問(wèn)題的解 190
2.2.1 超定問(wèn)題的最小方差解 190
2.2.2 欠定問(wèn)題的簡(jiǎn)單模型解 191
2.2.3 混定問(wèn)題的解——馬夸特解 192
2.3 先驗(yàn)信息在模型構(gòu)制中的應(yīng)用 193
2.3.1 對(duì)模型參數(shù)的限制 193
2.3.2 對(duì)觀測(cè)數(shù)據(jù)的限制 193
2.4 觀測(cè)數(shù)據(jù)和模型參數(shù)估算值的方差 194
2.5 廣義逆與廣義逆反演法 195
2.5.1 廣義逆矩陣的概念 195
2.5.2 奇異值分解和自然逆 197
2.5.3 廣義逆反演法 200
2.6 廣義逆反演解的評(píng)價(jià) 201
2.6.1 數(shù)據(jù)分辨矩陣 201
2.6.2 參數(shù)分辨矩陣 203
2.6.3 數(shù)據(jù)重建、模型分辨以及解的誤差與特征值(奇異值)的關(guān)系 203
2.6.4 最佳折中解 205
2.7 奧卡姆(Occam)反演算法 206
2.7.1 Occam反演算法概述 206
2.7.2 Occam反演中求取拉格朗日乘子方法的改進(jìn) 207
參考文獻(xiàn) 217
應(yīng)用篇
第3章 2.5 維復(fù)電阻率法反演理論及應(yīng)用 221
3.1 復(fù)電阻率法反演理論 221
3.1.1 廣義最小二乘反演 222
3.1.2 電場(chǎng)振幅與電場(chǎng)相位聯(lián)合反演 223
3.1.3 多個(gè)排列數(shù)據(jù)的聯(lián)合反演 224
3.2 反演模型參數(shù)約束 224
3.2.1 模型光滑度約束 224
3.2.2 模型參數(shù)界限約束 226
3.3 靈敏度矩陣的推導(dǎo)與計(jì)算 227
3.3.1 電場(chǎng)振幅和電場(chǎng)相位靈敏度表達(dá)式 227
3.3.2 科爾-科爾(Cole-Cole)模型偏導(dǎo)數(shù)表達(dá)式 228
3.3.3 互換定理計(jì)算靈敏度 229
3.4 理論模型反演算例 231
3.4.1 山谷地形下低阻極化體反演 232
3.4.2 傾斜地形下高阻極化體反演 235
3.4.3 山脊地形下組合極化體反演 238
3.4.4 反演結(jié)果的分析與討論 241
3.5 廬樅沙溪斑巖銅礦的反演應(yīng)用 242
3.5.1 礦床地質(zhì)概況 242
3.5.2 野外數(shù)據(jù)采集 242
3.5.3 反演結(jié)果與分析 243
3.6 本章小結(jié) 247
參考文獻(xiàn) 247
第4章 帶地形三維大地電磁法反演理論及應(yīng)用 249
4.1 大地電磁測(cè)深法三維共軛梯度反演 249
4.1.1 目標(biāo)函數(shù) 249
4.1.2 反演流程 250
4.1.3 基于矢量有限元的雅可比偏導(dǎo)數(shù)計(jì)算 252
4.1.4 雅可比偏導(dǎo)數(shù)矩陣與一個(gè)向量乘積的“擬正演”問(wèn)題 255
4.2 理論模型反演算例 256
4.2.1 山谷地形下含低阻體模型 256
4.2.2 山峰地形下含高阻體模型 260
4.2.3 峰谷組合地形下含低阻體模型 264
4.2.4 數(shù)據(jù)含 10%高斯噪聲的三維反演 265
4.3 某礦區(qū)實(shí)測(cè)資料的三維反演應(yīng)用 266
4.3.1 關(guān)于某礦區(qū)資料 266
4.3.2 對(duì)反演所得三維電阻率模型初步分析 267
4.3.3 新疆某地區(qū)實(shí)測(cè)資料的三維反演 268
參考文獻(xiàn) 272
第5章 帶地形三維可控源電磁法反演理論 275
5.1 可控源電磁法三維共軛梯度反演 275
5.1.1 反演目標(biāo)函數(shù) 275
5.1.2 非線性共軛梯度法 276
5.1.3 目標(biāo)函數(shù)的梯度計(jì)算 276
5.1.4 反演步長(zhǎng)的線性搜索 278
5.2 理論模型反演算例 279
5.2.1 模型 1 280
5.2.2 模型 2 281
5.3 誤差分析 282
參考文獻(xiàn) 283
第6章 三維復(fù)電阻率反演理論及應(yīng)用 285
6.1 擬線性反演 285
6.1.1 反演物質(zhì)性質(zhì)系數(shù) 285
6.1.2 科爾-科爾(Cole-Cole)模型參數(shù)的確定 286
6.1.3 偏導(dǎo)數(shù)矩陣 287
6.2 理論模型反演算例 288
6.2.1 軸向偶極-偶極數(shù)據(jù)反演 288
6.2.2 中梯測(cè)深數(shù)據(jù)反演 292
6.3 實(shí)測(cè)數(shù)據(jù)反演 297
6.3.1 區(qū)域地質(zhì)背景及礦區(qū)地質(zhì)特征 297
6.3.2 花牛山礦區(qū)地質(zhì)概況 298
6.3.3 反演結(jié)果及分析 299
6.4 小結(jié) 302
參考文獻(xiàn) 302
第7章 遠(yuǎn)震層析成像方法理論及應(yīng)用 303
7.1 模型參數(shù)化 304
7.2 正演計(jì)算 306
7.3 反演計(jì)算 308
7.4 棋盤(pán)格測(cè)試與誤差分析 313
7.5 南海周邊東南亞地區(qū)實(shí)際數(shù)據(jù)處理 314
7.5.1 數(shù)據(jù) 314
7.5.2 棋盤(pán)測(cè)試與精度分析 314
7.5.3 層析成像結(jié)果與分析 315
參考文獻(xiàn) 318
第8章 聯(lián)合反演理論及應(yīng)用 322
8.1 雙物性單約束的交叉梯度聯(lián)合反演 322
8.2 雙物性聯(lián)合反演理論模型試驗(yàn) 325
8.2.1 大地電磁和地震走時(shí)聯(lián)合反演 325
8.2.2 大地電磁和磁法聯(lián)合反演 327
8.2.3 大地電磁和重力與大地電磁和重力梯度聯(lián)合反演對(duì)比分析 329
8.3 多物性多約束的交叉梯度聯(lián)合反演 334
8.4 多物性聯(lián)合反演理論模型試驗(yàn) 336
8.4.1 多物性多約束聯(lián)合反演的優(yōu)勢(shì)分析 336
8.4.2 歸一化算子的影響分析 338
8.5 RGB 合成圖 340
8.6 實(shí)測(cè)數(shù)據(jù)反演 341
8.7 小結(jié) 345
參考文獻(xiàn) 346
結(jié)束語(yǔ) 347