2.1 .CƠ SỞ LỰA CHỌN MƠ HÌNH
2.2. CÁCH THỨC TÍNH MƠ HÌNH
2.2.2 Mơ hình sai phân hữu hạn (Finite Difference (FD))
Cơng cụ mơ hình sai phân hữu hạn (FD) cho phép mơ phịng sự lan truyền sóng
địa chấn hoặc sóng điện từ tương ứng bằng phương tiện của phương pháp FD-method
cho các nguồn khác nhau (mặt sóng, nguồn điểm cũng như “Exploding-Reflector”- source). Khi một kết quả một tuyến đơn hoặc trường sóng phức tạp được lưu lại và hiển thị sau đó. Trong phần sau đây chúng ta sẽ mô tả việc mô phỏng GPR cho mặt cắt 2D cửa sổ 0 (standard GPR-data acquisition).
B1. Trước hết phải tạo ra một mơ hình mới hoặc phải tải một mơ hình đã có bằng cách sử dụng tùy chọn file/load model.
B2. Kích hoạt tùy chọn FD
B3. Hộp FD-GroupBox mở ra (xem hình dưới). Trong nhóm hộp này, có thể
đưa vào các tham số FD cần thiết.
B4. Trước hết, phải nhập vào tần số chình frequency cho sự mô phỏng theo đơn vị Hz (địa chấn) hoặc MHz (GPR).
B5. Nhập vào loại nguồn cần thiết source type – với ví dụ này cho ranh giới
phản xạ. Loại ranh giới phản xạ cho phép mô phỏng mặt cắt 2D (cửa sổ 0) trong một tính tốn.
B6. Nhập vào giá trị số gia DeltaX trong không gian theo phương x và theo phương z. Đối với việc tính tốn FD, mơ phỏng lớp phải được raster hóa với một
gia số đã cho theo phương x và phương z (tùy chọn DeltaX). Cũng như vậy, gia số thời gian cũng phải được gán (tùy chọn DeltaT). Kích thước của gia số khơng gian và thời gian phụ thuộc vào độ dài sóng cực tiểu cũng như vận tốc truyền sóng. Chương trình sẽ tự động xác định giá trị tới hạn của gia số khơng gian (1/8 độ dài sóng cực tiểu cho FD-scheme 4-space, 1/12 độ dài sóng cực tiểu cho FD-scheme 2 space, một cách tương ứng) và trình bày giá trị này trong hộp giá trị tới hạn kéo xuống bên phải, giá trị này không thể bỏ qua. Việc chọn DeltaX quá lớn thì làm tăng kết quả trong sự phân tán của sóng, vì vậy nếu kết quả chứng tỏ một sự phân tán như vậy, cần phải giảm gia số X.
B7. Nhập vào gia số DeltaT cần thiết. Gia số thời gian cực đại phụ thuộc vào vận tốc cực đại, cũng như gia số không gian đã cho DeltaX (một cách xấp xỉ
xi<=1./%2V cho sự lan truyền sóng điện từ và t<=1./(Vp+Vs) cho sự lan truyền sóng
đàn hồi). Gia số thời gian tới hạn của tập hợp dữ liệu hiện hành được trình bày trong
hộp thoại cá giá trị tới hạn tính tốn ở đáy. Việc gia tăng DeltaT quá lớn, làm gia tăng tính khơng ổn định của kết quả, điều này có ý nghĩa là sự gia tăng biên độ tỉ lệ với thời gian theo quy luật hàm mũ. Vì vậy, nếu các kết quả chứng tỏ sự gia tăng biên độ như vật, cần phải giảm gia số DeltaT.
B8.Nhập vào thời gian tổng cho sự mô phỏng Tmax.
B9.Nhập các điều kiện biên boundary condition, ví dụ giới hạn hấp thu cho sự mô phỏng GPR).