Số liệu đo qi(zj) được đưa vào theo phương phỏp đo, số điểm đo, kớch thước a theo trỡnh tự đo, sau đú chương trỡnh chạy theo thuật toỏn sau:
Theo phương trỡnh (3.6) để giảm sai số εi cú thể sử dụng phương
trỡnh Gauxơ - Newton xỏc định sự thay đổi tham số mụ hỡnh để giảm tổng
sai số bỡnh phương của phương trỡnh (3.7): 1
k i
p J
Với sự thay đổi tham số này, tớnh toỏn lại vectơ thay đổi tham số
nhận được mụ hỡnh mới:
1 1
n n n
g g p (3.8)
Tuy nhiờn trong thực tế bài toỏn ngược Địa vật lý ma trận J-1J cú thể
là kỳ dị nờn phương trỡnh bỡnh phương tối thiểu (3.6) khụng cú lời giải.
Nếu kết quả đo qi(zj) cú sai số lớn, ma trận sản phẩm J-1J gần kỳ dị,
kết quả giải bài toỏn theo phương trỡnh (3.8) sẽ sai với thực tế; trong trường
hợp này phải dựng phương phỏp hồi quy thực hiện theo phương phỏp bỡnh
phương tối thiểu cưỡng bức với phương trỡnh mới:
1 1 (J .J f) p Ji (3.9) Ở đõy: - 1 1 . . x x z z
- fx là bộ lọc theo phương ngang
- fz là bộ lọc theo phương thẳng đứng
- λ là hệ số suy giảm hay trọng số của cỏc bộ lọc làm trơn
Sau khi εi=εmin đạt kết quả, mụ hỡnh gn được chấp nhận là mụ hỡnh kết quả, chương trỡnh giải bài toỏn ngược sẽ dừng lại, trờn màn hỡnh mỏy hiện lờn:
- Lỏt cắt điện trở suất biểu kiến qi = ρk(x, z) của tài liệu đo.
- Lỏt cắt điện trở suất được tớnh toỏn ở lần lặp thứ i: gi(x, z)
- Chấp nhận εmin thỡ cho mụ hỡnh ρ(x, z)=gn(x, z) ở lần tớnh lặp thứ n
là mụ hỡnh kết quả.