Bagaimana cara mengubah Distance, Azimuth, Dip ke XYZ?

Saya memiliki spreadsheet excel dengan data bor dan survei. Data header berisi ID lubang, dan koordinat lokasi, dan data survei berisi survei downhole terkait dengan nilai Distance, Azimuth, dan Dip.

Karena saya tahu lokasi lubang dan ketinggian permukaan, saya ingin dapat mengubah tabel survei menjadi koordinat XYZ juga. Apakah ada yang punya fungsi / prosedur / contoh? (VB & ArcObjects)

Data Header:

Data Header

Data Survei:

masukkan deskripsi gambar di sini

— Jakub Sisak GeoGraphics
sumber

Tidak bisakah Anda menggunakannya sebagai koordinat xyz? Saya berasumsi bahwa dip memiliki asal 0, bukan?

— Emily

@ Emily - Ya XYZ pertama diberikan (X: 425990, Y: 5409010, Z: 350). Nilai dip adalah 0 dan jaraknya 0. Bagaimana cara menghitung XYZ untuk Jarak: 41, Azimuth: 359, Dip: -71? (Pada setiap titik survei arah dan kemiringan akan berbeda, menghasilkan lubang yang melenceng dan berputar) Mungkin ada formula sederhana ...

— Jakub Sisak GeoGraphics

Jawaban:

Pertanyaannya meminta konversi antara koordinat bola dan koordinat kartesius . Spreadsheet ini menjabarkan formula:

Tangkapan layar spreadsheet

Garis biru adalah input, hitam adalah perhitungan antara, dan merah adalah output. Di dalam rumus, nilai-nilai dirujuk oleh nama-nama di kolom [Parameter] (ditugaskan melalui operasi Sisipkan | Nama | Buat).

Mereka berbeda dari yang ada di sebagian besar referensi matematika / fisika karena dalam geografi, azimuth biasanya diambil timur dari utara daripada utara timur. Ini membuat azimuth geografis menjadi komplemen dari yang matematika (jumlahnya mencapai 90 derajat). Mengganti sudut dengan komplemennya dalam fungsi trigonometri apa saja menukarnya dengan mitra "co" -nya: sinus dan cosinus dipertukarkan, bersinggungan dan cotangent, garis potong dan kosecant. Juga, dalam banyak sistem matematika, "dip" diekspresikan sebagai sudut dari true vertical (co-latitude) daripada sebagai angle dari horizontal (a latitude), lagi-lagi menyebabkan pertukaran sinus dan kosinus.

Edit 9/20/13

Untuk jarak downhole Anda mungkin ingin meniadakan dZ.

— whuber
sumber

Terima kasih. Kelihatan bagus! Harus mengkonfirmasi jenis azimuth. Bagaimana cara menemukan radian dalam formula?

— Jakub Sisak GeoGraphics

Radian = Derajat / 180 * Pi

— whuber

Tepat. Terima kasih. Aku menatap formula di tautan di atas dan menggaruk kepalaku. Jelas seperti bel sekarang.

— Jakub Sisak GeoGraphics

@ Jakub Anda sadar ini adalah spreadsheet yang berfungsi, valid, bukan? Anda bisa mengetik rumus ini ke Excel, beri nama sel di kolom [Nilai] seperti yang ditunjukkan di kolom [Parameter] di sebelah kirinya, dan itu akan berjalan. Teks biru adalah input; hitam adalah perhitungan antara; dan merah adalah output. Setelah merasa nyaman dengan itu, Anda dapat memodifikasi spreadsheet kedua untuk melakukan perhitungan untuk setiap entri. Satu-satunya trik adalah bergabung dengan koordinat (X0, Y0, Z0) dari yang pertama: lakukan dengan basis data atau melalui VLOOKUP ().

— whuber

Saya tidak sadar. Bahkan lebih baik! Jadi RADIANS jelas merupakan fungsi Excel yang saya lewati baik Azimuth atau dip.

— Jakub Sisak GeoGraphics

Meskipun ini adalah pertanyaan lama, jawaban lainnya tidak tepat. Konversi Jarak (Kedalaman Terukur), Dip (Inklinasi), Koordinat Azimuth ke 3D tergantung pada bagaimana Anda menginterpretasikan apa yang terjadi antara lokasi di mana pengukuran dilakukan (stasiun survei). Praktik standar saat ini adalah "Kurva Minimum" di mana asumsinya adalah busur lingkaran menghubungkan setiap lokasi survei.

http://www.drillingformulas.com/minimum-curvature-method/ memberikan detail lengkap tentang cara menghitung lokasi X, Y, dan Z. Bagian yang relevan adalah:

dMD = Distance2 - Distance1
B = acos(cos(I2 - I1) - (sin(I1)*sin(I2)*(1-cos(A2-A1))))
RF = 2 / B * tan(B / 2)
dX = dMD/2 * (sin(I1)*sin(A1) + sin(I2)*sin(A2))*RF
dY = dMD/2 * (sin(I1)*cos(A1) + sin(I2)*cos(A2))*RF
dZ = dMD/2 * (cos(I1) + cos(I2))*RF

X2 = X1 + dX
Y2 = Y1 + dX
Z2 = Z1 + dX

— Eric
sumber

Saya menemukan bahwa jika ada segmen lurus (I1 == I2 dan A1 == A2), maka B berakhir 0 dan sehingga RF menghasilkan kesenjangan dengan nol kesalahan. Dalam hal B adalah nol, maka saya mengatur RF ke 1 (karena batas RF sebagai B pergi ke 0 adalah 1). Kalau tidak, ini bekerja dengan baik, terima kasih!

— Saxon Druce