GemPy \u5730\u8d28\u5efa\u6a21\u6307\u5357 | GemPy Geological Modeling Guide<\/title>\n <style>\n :root {\n --bg: #f8f9fa;\n --surface: #ffffff;\n --ink: #1a1a2e;\n --muted: #5a6578;\n --rule: #dee2e6;\n --accent: #2563eb;\n --accent2: #059669;\n --code-bg: #1e293b;\n --code-text: #e2e8f0;\n --warn-bg: #fff7ed;\n --warn-border: #f97316;\n --info-bg: #eff6ff;\n --info-border: #3b82f6;\n }\n\n [data-theme=\"dark\"] {\n --bg: #0f172a;\n --surface: #1e293b;\n --ink: #f1f5f9;\n --muted: #94a3b8;\n --rule: #334155;\n --accent: #60a5fa;\n --accent2: #34d399;\n --code-bg: #020617;\n --code-text: #cbd5e1;\n --warn-bg: #2a1b0f;\n --warn-border: #f97316;\n --info-bg: #172554;\n --info-border: #3b82f6;\n }\n\n * { margin: 0; padding: 0; box-sizing: border-box; }\n\n body {\n font-family: -apple-system, BlinkMacSystemFont, \"Segoe UI\", Roboto, \"Helvetica Neue\", Arial, \"Noto Sans SC\", \"PingFang SC\", \"Microsoft YaHei\", sans-serif;\n background: var(--bg);\n color: var(--ink);\n line-height: 1.7;\n transition: background 0.3s, color 0.3s;\n }\n\n .container {\n max-width: 960px;\n margin: 0 auto;\n padding: 0 20px;\n }\n\n \/* Header *\/\n header {\n background: var(--surface);\n border-bottom: 1px solid var(--rule);\n padding: 1rem 0;\n position: sticky;\n top: 0;\n z-index: 100;\n }\n\n .header-inner {\n display: flex;\n justify-content: space-between;\n align-items: center;\n }\n\n .logo {\n font-size: 1.25rem;\n font-weight: 700;\n color: var(--accent);\n display: flex;\n align-items: center;\n gap: 0.5rem;\n }\n\n .logo svg {\n width: 28px;\n height: 28px;\n }\n\n .controls {\n display: flex;\n gap: 0.75rem;\n align-items: center;\n }\n\n .lang-switch {\n display: flex;\n background: var(--bg);\n border: 1px solid var(--rule);\n border-radius: 6px;\n overflow: hidden;\n }\n\n .lang-switch button {\n border: none;\n background: transparent;\n color: var(--muted);\n padding: 0.4rem 0.9rem;\n cursor: pointer;\n font-size: 0.85rem;\n transition: all 0.2s;\n }\n\n .lang-switch button.active {\n background: var(--accent);\n color: white;\n }\n\n .theme-toggle {\n background: var(--bg);\n border: 1px solid var(--rule);\n border-radius: 6px;\n padding: 0.4rem 0.7rem;\n cursor: pointer;\n color: var(--muted);\n font-size: 1rem;\n line-height: 1;\n }\n\n \/* Hero *\/\n .hero {\n padding: 3rem 0 2rem;\n text-align: center;\n }\n\n .hero h1 {\n font-size: 2rem;\n margin-bottom: 0.75rem;\n line-height: 1.3;\n }\n\n .hero p {\n color: var(--muted);\n font-size: 1.05rem;\n max-width: 640px;\n margin: 0 auto;\n }\n\n .badge {\n display: inline-block;\n background: var(--accent);\n color: white;\n padding: 0.25rem 0.75rem;\n border-radius: 999px;\n font-size: 0.75rem;\n font-weight: 600;\n margin-bottom: 1rem;\n text-transform: uppercase;\n letter-spacing: 0.05em;\n }\n\n \/* Sections *\/\n section {\n background: var(--surface);\n border: 1px solid var(--rule);\n border-radius: 10px;\n padding: 1.75rem;\n margin-bottom: 1.5rem;\n }\n\n section h2 {\n font-size: 1.35rem;\n margin-bottom: 1rem;\n padding-bottom: 0.5rem;\n border-bottom: 2px solid var(--accent);\n display: inline-block;\n }\n\n section h3 {\n font-size: 1.1rem;\n margin: 1.5rem 0 0.75rem;\n color: var(--accent);\n }\n\n section p, section li {\n margin-bottom: 0.75rem;\n color: var(--muted);\n }\n\n ul, ol {\n padding-left: 1.5rem;\n margin-bottom: 1rem;\n }\n\n li { margin-bottom: 0.4rem; }\n\n \/* Tables *\/\n .table-wrap {\n overflow-x: auto;\n margin: 1rem 0;\n }\n\n table {\n width: 100%;\n border-collapse: collapse;\n font-size: 0.9rem;\n min-width: 500px;\n }\n\n th, td {\n padding: 0.6rem 0.8rem;\n text-align: left;\n border-bottom: 1px solid var(--rule);\n }\n\n th {\n background: var(--bg);\n font-weight: 600;\n color: var(--ink);\n }\n\n tr:hover td { background: var(--bg); }\n\n \/* Code blocks *\/\n pre {\n background: var(--code-bg);\n color: var(--code-text);\n padding: 1.25rem;\n border-radius: 8px;\n overflow-x: auto;\n font-family: \"SF Mono\", Monaco, \"Cascadia Code\", Consolas, monospace;\n font-size: 0.85rem;\n line-height: 1.6;\n margin: 1rem 0;\n }\n\n code {\n background: var(--bg);\n color: var(--accent);\n padding: 0.15rem 0.4rem;\n border-radius: 4px;\n font-family: \"SF Mono\", Monaco, \"Cascadia Code\", Consolas, monospace;\n font-size: 0.88em;\n }\n\n pre code {\n background: transparent;\n color: inherit;\n padding: 0;\n }\n\n \/* Callouts *\/\n .callout {\n padding: 1rem 1.25rem;\n border-radius: 8px;\n margin: 1rem 0;\n border-left: 4px solid;\n }\n\n .callout.info {\n background: var(--info-bg);\n border-color: var(--info-border);\n }\n\n .callout.warn {\n background: var(--warn-bg);\n border-color: var(--warn-border);\n }\n\n .callout strong {\n display: block;\n margin-bottom: 0.3rem;\n color: var(--ink);\n }\n\n \/* Links *\/\n a {\n color: var(--accent);\n text-decoration: none;\n }\n\n a:hover { text-decoration: underline; }\n\n .link-list a {\n display: block;\n padding: 0.6rem 0;\n border-bottom: 1px solid var(--rule);\n }\n\n .link-list a:last-child { border-bottom: none; }\n\n \/* Tags *\/\n .tag {\n display: inline-block;\n background: var(--bg);\n border: 1px solid var(--rule);\n padding: 0.2rem 0.6rem;\n border-radius: 4px;\n font-size: 0.8rem;\n color: var(--muted);\n margin-right: 0.4rem;\n margin-bottom: 0.4rem;\n }\n\n .tag.required {\n background: var(--accent);\n color: white;\n border-color: var(--accent);\n }\n\n .tag.optional {\n background: var(--accent2);\n color: white;\n border-color: var(--accent2);\n }\n\n \/* Footer *\/\n footer {\n text-align: center;\n padding: 2rem 0;\n color: var(--muted);\n font-size: 0.85rem;\n }\n\n \/* Responsive *\/\n @media (max-width: 600px) {\n .hero h1 { font-size: 1.5rem; }\n section { padding: 1.25rem; }\n .header-inner { flex-wrap: wrap; gap: 0.5rem; }\n }\n\n \/* Hidden language sections *\/\n .lang-section { display: none; }\n .lang-section.active { display: block; }\n <\/style>\n<\/head>\n<body>\n <header>\n <div class=\"container header-inner\">\n <div class=\"logo\">\n <svg viewBox=\"0 0 24 24\" fill=\"none\" stroke=\"currentColor\" stroke-width=\"2\" stroke-linecap=\"round\" stroke-linejoin=\"round\">\n <path d=\"M12 2L2 7l10 5 10-5-10-5z\"\/>\n <path d=\"M2 17l10 5 10-5\"\/>\n <path d=\"M2 12l10 5 10-5\"\/>\n <\/svg>\n GemPy Guide\n <\/div>\n <div class=\"controls\">\n <div class=\"lang-switch\">\n <button class=\"active\" onclick=\"setLang('zh')\" id=\"btn-zh\">\u4e2d\u6587<\/button>\n <button onclick=\"setLang('en')\" id=\"btn-en\">English<\/button>\n <\/div>\n <button class=\"theme-toggle\" onclick=\"toggleTheme()\" title=\"Toggle theme\">☼<\/button>\n <\/div>\n <\/div>\n <\/header>\n\n <div class=\"container\">\n \n <div class=\"lang-section active\" id=\"lang-zh\">\n <div class=\"hero\">\n <span class=\"badge\">GemPy v3 (2024)<\/span>\n <h1>GemPy \u5730\u8d28\u5efa\u6a21\u6570\u636e\u51c6\u5907\u6307\u5357<\/h1>\n <p>\u57fa\u4e8e\u5f00\u6e90\u9690\u5f0f\u4e09\u7ef4\u5730\u8d28\u5efa\u6a21\u5e93 GemPy\uff0c\u4ecb\u7ecd\u4ece\u94bb\u5b54\u67f1\u72b6\u56fe\u6570\u636e\u5230 3D \u5730\u8d28\u6a21\u578b\u7684\u5b8c\u6574\u6570\u636e\u51c6\u5907\u6d41\u7a0b\u4e0e\u57fa\u672c\u5efa\u6a21\u793a\u4f8b\u3002<\/p>\n <\/div>\n\n <section>\n <h2>\u4e00\u3001GemPy \u7b80\u4ecb<\/h2>\n <p><a href=\"https:\/\/www.gempy.org\/\" target=\"_blank\">GemPy<\/a> \u662f\u4e00\u4e2a\u57fa\u4e8e Python \u7684\u5f00\u6e90\u4e09\u7ef4\u7ed3\u6784\u5730\u8d28\u5efa\u6a21\u5e93\uff0c\u91c7\u7528<strong>\u9690\u5f0f\u5efa\u6a21\uff08Implicit Modeling\uff09<\/strong>\u65b9\u6cd5\uff0c\u901a\u8fc7\u901a\u7528\u534f\u514b\u91cc\u91d1\u63d2\u503c\u6784\u5efa 3D \u6807\u91cf\u573a\u3002\u5b83\u80fd\u591f\u9ad8\u6548\u5730\u521b\u5efa\u5305\u542b\u8936\u76b1\u3001\u65ad\u5c42\u3001\u4e0d\u6574\u5408\u7b49\u590d\u6742\u6784\u9020\u7279\u5f81\u7684\u5730\u8d28\u6a21\u578b\uff0c\u5e76\u539f\u751f\u652f\u6301\u6982\u7387\u5efa\u6a21\u4e0e\u4e0d\u786e\u5b9a\u6027\u5206\u6790\u3002<\/p>\n <div class=\"callout info\">\n <strong>\u6838\u5fc3\u7279\u70b9<\/strong>\n GPU \u52a0\u901f\u8ba1\u7b97 | \u652f\u6301 PyTorch\/NumPy \u540e\u7aef | 2D\/3D \u53ef\u89c6\u5316\uff08PyVista\uff09| \u53ef\u76f4\u63a5\u5bfc\u51fa\u81f3 Blender \/ LiquidEarth\n <\/div>\n <\/section>\n\n <section>\n <h2>\u4e8c\u3001\u5efa\u6a21\u6240\u9700\u7684\u6700\u5c0f\u6570\u636e\u96c6<\/h2>\n <p>\u5728 GemPy \u4e2d\u6784\u5efa\u5730\u8d28\u6a21\u578b\uff0c\u9700\u8981\u4e24\u7c7b\u6838\u5fc3\u8f93\u5165\u6570\u636e\u3002\u8fd9\u4e9b\u6570\u636e\u53ef\u4ee5\u4ece\u94bb\u5b54\u67f1\u72b6\u56fe\u3001\u5730\u8d28\u5256\u9762\u56fe\u3001\u91ce\u5916\u9732\u5934\u6d4b\u91cf\u7b49\u8d44\u6599\u4e2d\u63d0\u53d6\u3002<\/p>\n\n <h3>1. \u754c\u9762\u70b9\uff08Surface Points\uff09<\/h3>\n <p>\u754c\u9762\u70b9\u662f\u4e09\u7ef4\u7a7a\u95f4\u4e2d\u6807\u8bb0\u5730\u8d28\u754c\u9762\uff08\u5982\u5730\u5c42\u5e95\u754c\u9762\u3001\u65ad\u5c42\u9762\u3001\u4e0d\u6574\u5408\u9762\uff09\u4f4d\u7f6e\u7684\u5750\u6807\u70b9\u3002<strong>\u5728 GemPy \u4e2d\uff0c\u6bcf\u4e2a\u754c\u9762\u4ee3\u8868\u5176\u6240\u8d4b\u5355\u5143\u7684\u5e95\u754c\u9762\u3002<\/strong><\/p>\n <div class=\"table-wrap\">\n <table>\n <tr>\n <th>\u5b57\u6bb5\u540d<\/th>\n <th>\u7c7b\u578b<\/th>\n <th>\u8bf4\u660e<\/th>\n <th>\u5fc5\u586b<\/th>\n <\/tr>\n <tr>\n <td><code>X<\/code><\/td>\n <td>\u6570\u503c<\/td>\n <td>\u4e1c\u5411\u5750\u6807\uff08\u4e0e\u5efa\u6a21\u8303\u56f4\u5355\u4f4d\u4e00\u81f4\uff09<\/td>\n <td><span class=\"tag required\">\u662f<\/span><\/td>\n <\/tr>\n <tr>\n <td><code>Y<\/code><\/td>\n <td>\u6570\u503c<\/td>\n <td>\u5317\u5411\u5750\u6807<\/td>\n <td><span class=\"tag required\">\u662f<\/span><\/td>\n <\/tr>\n <tr>\n <td><code>Z<\/code><\/td>\n <td>\u6570\u503c<\/td>\n <td>\u9ad8\u7a0b\/\u6df1\u5ea6\uff08\u901a\u5e38\u5730\u4e0b\u4e3a\u8d1f\u503c\uff09<\/td>\n <td><span class=\"tag required\">\u662f<\/span><\/td>\n <\/tr>\n <tr>\n <td><code>formation<\/code><\/td>\n <td>\u5b57\u7b26\u4e32<\/td>\n <td>\u5730\u5c42\/\u6784\u9020\u5355\u5143\u540d\u79f0\uff08\u5982 “Sandstone_1″\uff09<\/td>\n <td><span class=\"tag required\">\u662f<\/span><\/td>\n <\/tr>\n <\/table>\n <\/div>\n\n <h3>2. \u4ea7\u72b6\u6d4b\u91cf\uff08Orientations\uff09<\/h3>\n <p>\u4ea7\u72b6\u6570\u636e\u63cf\u8ff0\u5730\u8d28\u754c\u9762\u7684\u7a7a\u95f4\u503e\u659c\u65b9\u5411\uff0c\u7528\u4e8e\u7ea6\u675f\u63d2\u503c\u7684\u51e0\u4f55\u5f62\u6001\u3002\u6bcf\u4e2a\u6784\u9020\u7ec4\u81f3\u5c11\u9700\u8981\u4e00\u6761\u4ea7\u72b6\u6570\u636e\u3002<\/p>\n <div class=\"table-wrap\">\n <table>\n <tr>\n <th>\u5b57\u6bb5\u540d<\/th>\n <th>\u7c7b\u578b<\/th>\n <th>\u8bf4\u660e<\/th>\n <th>\u5fc5\u586b<\/th>\n <\/tr>\n <tr>\n <td><code>X<\/code>, <code>Y<\/code>, <code>Z<\/code><\/td>\n <td>\u6570\u503c<\/td>\n <td>\u6d4b\u91cf\u70b9\u7684\u4e09\u7ef4\u5750\u6807<\/td>\n <td><span class=\"tag required\">\u662f<\/span><\/td>\n <\/tr>\n <tr>\n <td><code>azimuth<\/code><\/td>\n <td>\u6570\u503c<\/td>\n <td>\u503e\u5411\uff080\u2013360\u5ea6\uff09<\/td>\n <td><span class=\"tag required\">\u662f<\/span><\/td>\n <\/tr>\n <tr>\n <td><code>dip<\/code><\/td>\n <td>\u6570\u503c<\/td>\n <td>\u503e\u89d2\uff080\u201390\u5ea6\uff09<\/td>\n <td><span class=\"tag required\">\u662f<\/span><\/td>\n <\/tr>\n <tr>\n <td><code>polarity<\/code><\/td>\n <td>\u6570\u503c<\/td>\n <td>\u6781\u6027\uff08\u901a\u5e38\u53d6 1\uff0c\u5012\u7f6e\u6784\u9020\u53ef\u53d6 -1\uff09<\/td>\n <td><span class=\"tag required\">\u662f<\/span><\/td>\n <\/tr>\n <tr>\n <td><code>formation<\/code><\/td>\n <td>\u5b57\u7b26\u4e32<\/td>\n <td>\u6240\u5c5e\u5730\u5c42\/\u6784\u9020\u5355\u5143\u540d\u79f0<\/td>\n <td><span class=\"tag required\">\u662f<\/span><\/td>\n <\/tr>\n <\/table>\n <\/div>\n\n <div class=\"callout warn\">\n <strong>\u6700\u5c0f\u6570\u636e\u89c4\u5219<\/strong>\n \u6bcf\u4e2a\u6784\u9020\u7ec4\uff08Series\/Stack\uff09\u4e2d\uff0c\u81f3\u5c11\u9700\u8981\u4e00\u4e2a\u754c\u9762\u7684<strong>\u4e24\u4e2a\u754c\u9762\u70b9<\/strong> + <strong>\u4e00\u6761\u4ea7\u72b6<\/strong>\u3002\u6dfb\u52a0\u66f4\u591a\u540c\u7ec4\u754c\u9762\u65f6\uff0c\u6bcf\u4e2a\u65b0\u754c\u9762\u53ea\u9700\u81f3\u5c11\u4e00\u4e2a\u754c\u9762\u70b9\u5373\u53ef\u5171\u4eab\u8be5\u7ec4\u7684\u4ea7\u72b6\u4fe1\u606f\u3002\n <\/div>\n <\/section>\n\n <section>\n <h2>\u4e09\u3001\u4ece\u94bb\u5b54\u67f1\u72b6\u56fe\u63d0\u53d6\u6570\u636e<\/h2>\n <p>\u94bb\u5b54\u67f1\u72b6\u56fe\u662f\u5730\u8d28\u5efa\u6a21\u6700\u5e38\u89c1\u7684\u539f\u59cb\u8d44\u6599\u3002\u4ee5\u4e0b\u662f\u4ece\u94bb\u5b54\u67f1\u72b6\u56fe\u8f6c\u6362\u4e3a GemPy \u8f93\u5165\u6570\u636e\u7684\u5b8c\u6574\u6d41\u7a0b\uff1a<\/p>\n\n <h3>\u6b65\u9aa4 1\uff1a\u6536\u96c6\u539f\u59cb\u94bb\u5b54\u8d44\u6599<\/h3>\n <p>\u7528\u6237\u9700\u8981\u63d0\u4f9b\u4ee5\u4e0b\u6750\u6599\uff1a<\/p>\n <ol>\n <li><strong>\u94bb\u5b54\u57fa\u672c\u4fe1\u606f\u8868<\/strong>\uff1a\u5b54\u53f7\u3001\u5b54\u53e3\u5750\u6807\uff08X, Y\uff09\u3001\u5b54\u53e3\u9ad8\u7a0b\uff08Z\uff09\u3001\u94bb\u5b54\u6df1\u5ea6\u3001\u94bb\u5b54\u503e\u89d2\u4e0e\u65b9\u4f4d\u89d2\uff08\u82e5\u4e3a\u659c\u5b54\uff09<\/li>\n <li><strong>\u94bb\u5b54\u5206\u5c42\u6570\u636e\u8868<\/strong>\uff1a\u6bcf\u5c42\u7684\u5c42\u5e95\u6df1\u5ea6\u3001\u5c42\u5e95\u9ad8\u7a0b\u3001\u5730\u5c42\u4ee3\u53f7\u3001\u5730\u5c42\u540d\u79f0\u3001\u5ca9\u6027\u63cf\u8ff0<\/li>\n <li><strong>\u5730\u5c42\u4ea7\u72b6\u6570\u636e<\/strong>\uff1a\u94bb\u5b54\u5185\u6216\u94bb\u5b54\u9644\u8fd1\u6d4b\u5f97\u7684\u5730\u5c42\u503e\u5411\u4e0e\u503e\u89d2<\/li>\n <li><strong>\u7814\u7a76\u533a\u8303\u56f4<\/strong>\uff1a\u5efa\u6a21\u7684\u5e73\u9762\u8303\u56f4\uff08X_min, X_max, Y_min, Y_max\uff09\u4e0e\u9ad8\u7a0b\u8303\u56f4\uff08Z_min, Z_max\uff09<\/li>\n <li><strong>\u5730\u5c42\u5e8f\u5217\u5173\u7cfb<\/strong>\uff1a\u5404\u5730\u5c42\u7684\u6c89\u79ef\u987a\u5e8f\u3001\u662f\u5426\u5b58\u5728\u4e0d\u6574\u5408\/\u4fb5\u5165\u5173\u7cfb\u3001\u65ad\u5c42\u4fe1\u606f<\/li>\n <\/ol>\n\n <h3>\u6b65\u9aa4 2\uff1a\u8ba1\u7b97\u754c\u9762\u70b9\u5750\u6807<\/h3>\n <p>\u5bf9\u6bcf\u53e3\u94bb\u5b54\u7684\u6bcf\u4e2a\u5730\u5c42\u5e95\u754c\u9762\uff0c\u8ba1\u7b97\u5176\u4e09\u7ef4\u5750\u6807\uff1a<\/p>\n <pre><code># \u76f4\u5b54\u793a\u4f8b\uff1a\u5b54\u53e3\u5750\u6807 (X0, Y0, Z0)\uff0c\u5c42\u5e95\u6df1\u5ea6 D\nX = X0\nY = Y0\nZ = Z0 - D # \u82e5 Z \u5411\u4e0b\u4e3a\u8d1f\n\n# \u659c\u5b54\u9700\u6839\u636e\u5929\u9876\u89d2 theta \u548c\u65b9\u4f4d\u89d2 alpha \u8fdb\u884c\u6362\u7b97\n# X = X0 + D * sin(theta) * sin(alpha)\n# Y = Y0 + D * sin(theta) * cos(alpha)\n# Z = Z0 - D * cos(theta)<\/code><\/pre>\n\n <h3>\u6b65\u9aa4 3\uff1a\u51c6\u5907 CSV \u6587\u4ef6<\/h3>\n <p>\u5c06\u63d0\u53d6\u7684\u754c\u9762\u70b9\u4fdd\u5b58\u4e3a <code>surface_points.csv<\/code>\uff1a<\/p>\n <pre><code>X,Y,Z,formation\n0,0,800,rock3\n0,0,450,rock2\n0,0,250,rock1\n500,0,800,rock3\n...<\/code><\/pre>\n\n <p>\u5c06\u4ea7\u72b6\u6570\u636e\u4fdd\u5b58\u4e3a <code>orientations.csv<\/code>\uff1a<\/p>\n <pre><code>X,Y,Z,azimuth,dip,polarity,formation\n100,500,650,90,315,1,rock2\n100,500,450,90,315,1,rock1\n...<\/code><\/pre>\n\n <div class=\"callout info\">\n <strong>\u5750\u6807\u7cfb\u63d0\u793a<\/strong>\n GemPy \u4f7f\u7528\u53f3\u624b\u5750\u6807\u7cfb\u3002Z \u8f74\u53ef\u5411\u4e0a\u4e3a\u6b63\uff08\u5730\u8868\u9ad8\u7a0b\uff09\u6216\u5411\u4e0b\u4e3a\u8d1f\uff08\u6df1\u5ea6\uff09\uff0c\u53ea\u9700\u5728\u6574\u4e2a\u6a21\u578b\u4e2d\u4fdd\u6301\u4e00\u81f4\u3002\u5efa\u8bae\u4e0e\u7814\u7a76\u533a\u5e38\u7528\u5750\u6807\u7cfb\u5bf9\u9f50\u3002\n <\/div>\n <\/section>\n\n <section>\n <h2>\u56db\u3001GemPy \u5efa\u6a21\u57fa\u672c\u6d41\u7a0b\uff08Demo\uff09<\/h2>\n\n <h3>1. \u5b89\u88c5 GemPy<\/h3>\n <pre><code># \u6838\u5fc3\u5e93\npip install gempy\n\n# \u53ef\u89c6\u5316\u6a21\u5757\npip install gempy-viewer\n\n# 3D \u53ef\u89c6\u5316\u4f9d\u8d56\uff08\u63a8\u8350\u5b89\u88c5\uff09\npip install pyvista\n\n# \u82e5\u9700 GPU \u52a0\u901f\uff08PyTorch \u540e\u7aef\uff09\npip install torch<\/code><\/pre>\n\n <h3>2. \u4ece CSV \u5bfc\u5165\u6570\u636e\u5e76\u521b\u5efa\u6a21\u578b<\/h3>\n <pre><code>import gempy as gp\nimport gempy_viewer as gpv\n\n# \u6570\u636e\u8def\u5f84\ndata_path = \"https:\/\/raw.githubusercontent.com\/cgre-aachen\/gempy_data\/master\/\"\npath_to_data = data_path + \"\/data\/input_data\/video_tutorials_v3\/\"\n\n# \u521b\u5efa\u5730\u8d28\u6a21\u578b\ngeo_model = gp.create_geomodel(\n project_name='tutorial_model',\n extent=[0, 2500, 0, 1000, 0, 1000],\n resolution=[100, 40, 40],\n importer_helper=gp.data.ImporterHelper(\n path_to_orientations=path_to_data + \"tutorial_model_orientations.csv\",\n path_to_surface_points=path_to_data + \"tutorial_model_surface_points.csv\"\n )\n)<\/code><\/pre>\n\n <h3>3. \u5b9a\u4e49\u6784\u9020\u5e8f\u5217\u4e0e\u6620\u5c04\u5173\u7cfb<\/h3>\n <pre><code># \u6620\u5c04\u5730\u5c42\u5230\u6784\u9020\u7ec4\uff08\u4ece\u5e74\u8f7b\u5230\u5e74\u8001\u6392\u5e8f\uff09\ngp.map_stack_to_surfaces(\n gempy_model=geo_model,\n mapping_object={\n \"Strat_Series2\": (\"rock3\"),\n \"Strat_Series1\": (\"rock2\", \"rock1\")\n }\n)<\/code><\/pre>\n\n <h3>4. \u8ba1\u7b97\u6a21\u578b<\/h3>\n <pre><code># \u8ba1\u7b97\u5730\u8d28\u6a21\u578b\nimport os\nos.environ[\"DEFAULT_BACKEND\"] = \"PYTORCH\" # \u6216 NUMPY\n\ngp.compute_model(geo_model)<\/code><\/pre>\n\n <h3>5. \u53ef\u89c6\u5316\u7ed3\u679c<\/h3>\n <pre><code># 2D \u5256\u9762\u53ef\u89c6\u5316\ngpv.plot_2d(geo_model, cell_number=20, show_data=True)\n\n# 3D \u4ea4\u4e92\u5f0f\u53ef\u89c6\u5316\ngpv.plot_3d(geo_model, show_lith=True, show_boundaries=True)<\/code><\/pre>\n\n <h3>6. \u4ece\u96f6\u624b\u52a8\u521b\u5efa\uff08\u4e0d\u4f9d\u8d56 CSV\uff09<\/h3>\n <pre><code># \u521b\u5efa\u7a7a\u6a21\u578b\ngeo_model = gp.create_geomodel(\n project_name='Model1',\n extent=[0, 780, -200, 200, -582, 0],\n resolution=(50, 50, 50),\n structural_frame=gp.data.StructuralFrame.initialize_default_structure()\n)\n\n# \u624b\u52a8\u6dfb\u52a0\u754c\u9762\u70b9\ngp.add_surface_points(\n geo_model=geo_model,\n x=[225, 460, 617],\n y=[0, 0, 0],\n z=[-95, -100, -10],\n elements_names=['Limestone', 'Limestone', 'Limestone']\n)\n\n# \u624b\u52a8\u6dfb\u52a0\u4ea7\u72b6\ngp.add_orientations(\n geo_model=geo_model,\n x=[350], y=[0], z=[-120],\n elements_names=['Limestone'],\n pole_vector=[[0, 0, 1]] # \u6c34\u5e73\u5c42\u7406\uff0c\u6cd5\u5411\u91cf\u671d\u4e0a\n)\n\n# \u8ba1\u7b97\u5e76\u53ef\u89c6\u5316\ngp.compute_model(geo_model)\ngpv.plot_2d(geo_model)<\/code><\/pre>\n <\/section>\n\n <section>\n <h2>\u4e94\u3001\u5b8c\u6574\u6570\u636e\u9700\u6c42\u6e05\u5355<\/h2>\n <p>\u4e3a\u4f7f\u7528 GemPy \u6210\u529f\u5efa\u7acb\u4e09\u7ef4\u5730\u8d28\u6a21\u578b\uff0c\u8bf7\u786e\u4fdd\u63d0\u4f9b\u4ee5\u4e0b\u8d44\u6599\uff1a<\/p>\n <div class=\"table-wrap\">\n <table>\n <tr>\n <th>\u6750\u6599<\/th>\n <th>\u683c\u5f0f<\/th>\n <th>\u5fc5\u8981\u6027<\/th>\n <th>\u8bf4\u660e<\/th>\n <\/tr>\n <tr>\n <td>\u94bb\u5b54\u5750\u6807\u4e0e\u5b54\u53e3\u9ad8\u7a0b<\/td>\n <td>Excel \/ CSV<\/td>\n <td><span class=\"tag required\">\u5fc5\u9700<\/span><\/td>\n <td>\u6bcf\u53e3\u5b54\u7684 X, Y, Z \u5750\u6807<\/td>\n <\/tr>\n <tr>\n <td>\u94bb\u5b54\u5206\u5c42\u6570\u636e<\/td>\n <td>Excel \/ CSV<\/td>\n <td><span class=\"tag required\">\u5fc5\u9700<\/span><\/td>\n <td>\u6bcf\u5c42\u5e95\u6df1\u3001\u5730\u5c42\u540d\u79f0\u3001\u5ca9\u6027<\/td>\n <\/tr>\n <tr>\n <td>\u5730\u5c42\u4ea7\u72b6\u6d4b\u91cf<\/td>\n <td>Excel \/ CSV<\/td>\n <td><span class=\"tag required\">\u5fc5\u9700<\/span><\/td>\n <td>\u503e\u5411\u3001\u503e\u89d2\u53ca\u6d4b\u91cf\u4f4d\u7f6e\u5750\u6807<\/td>\n <\/tr>\n <tr>\n <td>\u7814\u7a76\u533a\u8303\u56f4<\/td>\n <td>\u6587\u672c\u8bf4\u660e<\/td>\n <td><span class=\"tag required\">\u5fc5\u9700<\/span><\/td>\n <td>X\/Y\/Z \u7684\u6700\u5c0f\u6700\u5927\u503c<\/td>\n <\/tr>\n <tr>\n <td>\u5730\u5c42\u5e8f\u5217\u5173\u7cfb<\/td>\n <td>\u6587\u672c\/\u56fe\u8868<\/td>\n <td><span class=\"tag required\">\u5fc5\u9700<\/span><\/td>\n <td>\u65b0\u8001\u5173\u7cfb\u3001\u4e0d\u6574\u5408\u3001\u4fb5\u5165\u4f53<\/td>\n <\/tr>\n <tr>\n <td>\u65ad\u5c42\u6570\u636e<\/td>\n <td>Excel \/ CSV<\/td>\n <td><span class=\"tag optional\">\u53ef\u9009<\/span><\/td>\n <td>\u65ad\u5c42\u9762\u4ea7\u72b6\u3001\u65ad\u8ddd\u3001\u5ef6\u4f38\u8303\u56f4<\/td>\n <\/tr>\n <tr>\n <td>\u5730\u5f62\u6570\u636e<\/td>\n <td>DEM \/ GeoTIFF<\/td>\n <td><span class=\"tag optional\">\u53ef\u9009<\/span><\/td>\n <td>\u7814\u7a76\u533a\u6570\u5b57\u9ad8\u7a0b\u6a21\u578b<\/td>\n <\/tr>\n <tr>\n <td>\u5730\u8d28\u5256\u9762\u56fe<\/td>\n <td>\u56fe\u7247 \/ CAD<\/td>\n <td><span class=\"tag optional\">\u53ef\u9009<\/span><\/td>\n <td>\u8f85\u52a9\u9a8c\u8bc1\u6a21\u578b\u5408\u7406\u6027<\/td>\n <\/tr>\n <\/table>\n <\/div>\n <\/section>\n\n <section>\n <h2>\u516d\u3001\u5b98\u65b9\u6587\u6863\u4e0e\u5e2e\u52a9\u94fe\u63a5<\/h2>\n <div class=\"link-list\">\n <a href=\"https:\/\/www.gempy.org\/\" target=\"_blank\" rel=\"noopener\">GemPy \u5b98\u65b9\u7f51\u7ad9 — \u9879\u76ee\u4e3b\u9875\u4e0e\u529f\u80fd\u4ecb\u7ecd<\/a>\n <a href=\"https:\/\/docs.gempy.org\/\" target=\"_blank\" rel=\"noopener\">GemPy \u5b98\u65b9\u6587\u6863 — \u5b8c\u6574 API \u4e0e\u6559\u7a0b<\/a>\n <a href=\"https:\/\/docs.gempy.org\/installation.html\" target=\"_blank\" rel=\"noopener\">\u5b89\u88c5\u6307\u5357 — pip \/ conda \/ \u624b\u52a8\u5b89\u88c5<\/a>\n <a href=\"https:\/\/docs.gempy.org\/tutorials\/a_getting_started\/get_started.html\" target=\"_blank\" rel=\"noopener\">\u5165\u95e8\u6559\u7a0b — \u4ece\u96f6\u5f00\u59cb\u6784\u5efa\u7b2c\u4e00\u4e2a\u6a21\u578b<\/a>\n <a href=\"https:\/\/docs.gempy.org\/tutorials\/ch1_fundamentals\/ch1_1_basics.html\" target=\"_blank\" rel=\"noopener\">\u57fa\u7840\u5efa\u6a21\u6559\u7a0b — CSV \u6570\u636e\u5bfc\u5165\u4e0e\u6784\u9020\u5e8f\u5217<\/a>\n <a href=\"https:\/\/github.com\/cgre-aachen\/gempy\" target=\"_blank\" rel=\"noopener\">GitHub \u4ed3\u5e93 — \u6e90\u7801\u3001Issues \u4e0e\u793a\u4f8b\u6570\u636e<\/a>\n <a href=\"https:\/\/www.youtube.com\/@GemPy3D\" target=\"_blank\" rel=\"noopener\">YouTube \u89c6\u9891\u6559\u7a0b — \u5206\u6b65\u64cd\u4f5c\u6f14\u793a<\/a>\n <a href=\"https:\/\/www.geosci-model-dev.net\/12\/1\/2019\/gmd-12-1-2019.pdf\" target=\"_blank\" rel=\"noopener\">GemPy 1.0 \u5b66\u672f\u8bba\u6587 — \u65b9\u6cd5\u8bba\u4e0e\u79d1\u5b66\u80cc\u666f<\/a>\n <\/div>\n <\/section>\n <\/div>\n\n \n <div class=\"lang-section\" id=\"lang-en\">\n <div class=\"hero\">\n <span class=\"badge\">GemPy v3 (2024)<\/span>\n <h1>GemPy Geological Modeling Data Preparation Guide<\/h1>\n <p>A comprehensive guide to preparing borehole log data for 3D structural geological modeling using the open-source implicit modeling library GemPy.<\/p>\n <\/div>\n\n <section>\n <h2>1. Introduction to GemPy<\/h2>\n <p><a href=\"https:\/\/www.gempy.org\/\" target=\"_blank\">GemPy<\/a> is an open-source Python library for 3D structural geological modeling. It uses an <strong>implicit modeling<\/strong> approach based on universal cokriging to construct 3D scalar fields, enabling fast creation of complex models including folds, faults, and unconformities. It also natively supports probabilistic modeling and uncertainty analysis.<\/p>\n <div class=\"callout info\">\n <strong>Key Features<\/strong>\n GPU-accelerated computation | PyTorch\/NumPy backend | 2D\/3D visualization via PyVista | Export to Blender \/ LiquidEarth\n <\/div>\n <\/section>\n\n <section>\n <h2>2. Minimum Required Datasets<\/h2>\n <p>To build a geological model in GemPy, two core input data types are required. These can be extracted from borehole logs, geological cross-sections, and outcrop measurements.<\/p>\n\n <h3>2.1 Surface Points<\/h3>\n <p>Surface points are 3D coordinates marking geological interfaces (e.g., layer boundaries, fault planes, unconformities). <strong>In GemPy, each surface represents the bottom of the unit it is assigned to.<\/strong><\/p>\n <div class=\"table-wrap\">\n <table>\n <tr>\n <th>Field<\/th>\n <th>Type<\/th>\n <th>Description<\/th>\n <th>Required<\/th>\n <\/tr>\n <tr>\n <td><code>X<\/code><\/td>\n <td>Numeric<\/td>\n <td>Easting coordinate<\/td>\n <td><span class=\"tag required\">Yes<\/span><\/td>\n <\/tr>\n <tr>\n <td><code>Y<\/code><\/td>\n <td>Numeric<\/td>\n <td>Northing coordinate<\/td>\n <td><span class=\"tag required\">Yes<\/span><\/td>\n <\/tr>\n <tr>\n <td><code>Z<\/code><\/td>\n <td>Numeric<\/td>\n <td>Elevation \/ depth (negative for subsurface)<\/td>\n <td><span class=\"tag required\">Yes<\/span><\/td>\n <\/tr>\n <tr>\n <td><code>formation<\/code><\/td>\n <td>String<\/td>\n <td>Formation \/ structural element name<\/td>\n <td><span class=\"tag required\">Yes<\/span><\/td>\n <\/tr>\n <\/table>\n <\/div>\n\n <h3>2.2 Orientation Measurements<\/h3>\n <p>Orientation data describes the spatial dip of geological interfaces, constraining the geometry of interpolation. At least one orientation per structural group is required.<\/p>\n <div class=\"table-wrap\">\n <table>\n <tr>\n <th>Field<\/th>\n <th>Type<\/th>\n <th>Description<\/th>\n <th>Required<\/th>\n <\/tr>\n <tr>\n <td><code>X<\/code>, <code>Y<\/code>, <code>Z<\/code><\/td>\n <td>Numeric<\/td>\n <td>Measurement location coordinates<\/td>\n <td><span class=\"tag required\">Yes<\/span><\/td>\n <\/tr>\n <tr>\n <td><code>azimuth<\/code><\/td>\n <td>Numeric<\/td>\n <td>Dip direction (0\u2013360 degrees)<\/td>\n <td><span class=\"tag required\">Yes<\/span><\/td>\n <\/tr>\n <tr>\n <td><code>dip<\/code><\/td>\n <td>Numeric<\/td>\n <td>Dip angle (0\u201390 degrees)<\/td>\n <td><span class=\"tag required\">Yes<\/span><\/td>\n <\/tr>\n <tr>\n <td><code>polarity<\/code><\/td>\n <td>Numeric<\/td>\n <td>Polarity (usually 1; -1 for inverted structures)<\/td>\n <td><span class=\"tag required\">Yes<\/span><\/td>\n <\/tr>\n <tr>\n <td><code>formation<\/code><\/td>\n <td>String<\/td>\n <td>Associated formation \/ element name<\/td>\n <td><span class=\"tag required\">Yes<\/span><\/td>\n <\/tr>\n <\/table>\n <\/div>\n\n <div class=\"callout warn\">\n <strong>Minimum Data Rule<\/strong>\n Each structural group (Series\/Stack) requires at least <strong>two surface points<\/strong> for one surface + <strong>one orientation<\/strong>. Additional surfaces in the same group only need one surface point each, as they share orientation information.\n <\/div>\n <\/section>\n\n <section>\n <h2>3. Extracting Data from Borehole Logs<\/h2>\n <p>Borehole logs are the most common source of data for geological modeling. Below is the complete workflow to convert borehole data into GemPy inputs:<\/p>\n\n <h3>Step 1: Collect Raw Borehole Materials<\/h3>\n <p>The user should provide the following:<\/p>\n <ol>\n <li><strong>Borehole location table<\/strong>: Borehole ID, collar coordinates (X, Y), collar elevation (Z), total depth, inclination and azimuth (for deviated holes)<\/li>\n <li><strong>Stratigraphic interval table<\/strong>: Bottom depth of each layer, bottom elevation, formation code, formation name, lithology description<\/li>\n <li><strong>Structural orientation data<\/strong>: Dip direction and dip angle measured within or near the borehole<\/li>\n <li><strong>Study area extent<\/strong>: X\/Y\/Z minimum and maximum values defining the modeling domain<\/li>\n <li><strong>Stratigraphic sequence<\/strong>: Depositional order, presence of unconformities or intrusions, fault information<\/li>\n <\/ol>\n\n <h3>Step 2: Compute Surface Point Coordinates<\/h3>\n <p>For each layer bottom in each borehole, calculate the 3D coordinates:<\/p>\n <pre><code># Vertical borehole example: collar (X0, Y0, Z0), layer bottom depth D\nX = X0\nY = Y0\nZ = Z0 - D # if Z is negative downward\n\n# Deviated borehole conversion using zenith angle theta and azimuth alpha:\n# X = X0 + D * sin(theta) * sin(alpha)\n# Y = Y0 + D * sin(theta) * cos(alpha)\n# Z = Z0 - D * cos(theta)<\/code><\/pre>\n\n <h3>Step 3: Prepare CSV Files<\/h3>\n <p>Save extracted surface points as <code>surface_points.csv<\/code>:<\/p>\n <pre><code>X,Y,Z,formation\n0,0,800,rock3\n0,0,450,rock2\n0,0,250,rock1\n500,0,800,rock3\n...<\/code><\/pre>\n\n <p>Save orientation data as <code>orientations.csv<\/code>:<\/p>\n <pre><code>X,Y,Z,azimuth,dip,polarity,formation\n100,500,650,90,315,1,rock2\n100,500,450,90,315,1,rock1\n...<\/code><\/pre>\n\n <div class=\"callout info\">\n <strong>Coordinate System Note<\/strong>\n GemPy uses a right-handed coordinate system. The Z axis can be positive upward (elevation) or negative downward (depth); consistency across the model is what matters. Align with your study area’s standard coordinate system.\n <\/div>\n <\/section>\n\n <section>\n <h2>4. Basic GemPy Modeling Workflow (Demo)<\/h2>\n\n <h3>4.1 Install GemPy<\/h3>\n <pre><code># Core library\npip install gempy\n\n# Visualization module\npip install gempy-viewer\n\n# 3D visualization dependency (recommended)\npip install pyvista\n\n# For GPU acceleration (PyTorch backend)\npip install torch<\/code><\/pre>\n\n <h3>4.2 Import Data from CSV and Create Model<\/h3>\n <pre><code>import gempy as gp\nimport gempy_viewer as gpv\n\n# Data paths\ndata_path = \"https:\/\/raw.githubusercontent.com\/cgre-aachen\/gempy_data\/master\/\"\npath_to_data = data_path + \"\/data\/input_data\/video_tutorials_v3\/\"\n\n# Create geological model\ngeo_model = gp.create_geomodel(\n project_name='tutorial_model',\n extent=[0, 2500, 0, 1000, 0, 1000],\n resolution=[100, 40, 40],\n importer_helper=gp.data.ImporterHelper(\n path_to_orientations=path_to_data + \"tutorial_model_orientations.csv\",\n path_to_surface_points=path_to_data + \"tutorial_model_surface_points.csv\"\n )\n)<\/code><\/pre>\n\n <h3>4.3 Define Structural Sequence<\/h3>\n <pre><code># Map formations to structural groups (youngest to oldest)\ngp.map_stack_to_surfaces(\n gempy_model=geo_model,\n mapping_object={\n \"Strat_Series2\": (\"rock3\"),\n \"Strat_Series1\": (\"rock2\", \"rock1\")\n }\n)<\/code><\/pre>\n\n <h3>4.4 Compute the Model<\/h3>\n <pre><code>import os\nos.environ[\"DEFAULT_BACKEND\"] = \"PYTORCH\" # or \"NUMPY\"\n\ngp.compute_model(geo_model)<\/code><\/pre>\n\n <h3>4.5 Visualize Results<\/h3>\n <pre><code># 2D cross-section\ngpv.plot_2d(geo_model, cell_number=20, show_data=True)\n\n# 3D interactive visualization\ngpv.plot_3d(geo_model, show_lith=True, show_boundaries=True)<\/code><\/pre>\n\n <h3>4.6 Create from Scratch (Without CSV)<\/h3>\n <pre><code># Create empty model\ngeo_model = gp.create_geomodel(\n project_name='Model1',\n extent=[0, 780, -200, 200, -582, 0],\n resolution=(50, 50, 50),\n structural_frame=gp.data.StructuralFrame.initialize_default_structure()\n)\n\n# Add surface points manually\ngp.add_surface_points(\n geo_model=geo_model,\n x=[225, 460, 617],\n y=[0, 0, 0],\n z=[-95, -100, -10],\n elements_names=['Limestone', 'Limestone', 'Limestone']\n)\n\n# Add orientation manually\ngp.add_orientations(\n geo_model=geo_model,\n x=[350], y=[0], z=[-120],\n elements_names=['Limestone'],\n pole_vector=[[0, 0, 1]] # Horizontal bedding, normal vector up\n)\n\n# Compute and visualize\ngp.compute_model(geo_model)\ngpv.plot_2d(geo_model)<\/code><\/pre>\n <\/section>\n\n <section>\n <h2>5. Complete Data Requirements Checklist<\/h2>\n <p>To successfully build a 3D geological model with GemPy, please ensure the following materials are provided:<\/p>\n <div class=\"table-wrap\">\n <table>\n <tr>\n <th>Material<\/th>\n <th>Format<\/th>\n <th>Required<\/th>\n <th>Notes<\/th>\n <\/tr>\n <tr>\n <td>Borehole coordinates & collar elevation<\/td>\n <td>Excel \/ CSV<\/td>\n <td><span class=\"tag required\">Required<\/span><\/td>\n <td>X, Y, Z for each borehole<\/td>\n <\/tr>\n <tr>\n <td>Stratigraphic interval data<\/td>\n <td>Excel \/ CSV<\/td>\n <td><span class=\"tag required\">Required<\/span><\/td>\n <td>Bottom depth, formation name, lithology<\/td>\n <\/tr>\n <tr>\n <td>Structural orientation measurements<\/td>\n <td>Excel \/ CSV<\/td>\n <td><span class=\"tag required\">Required<\/span><\/td>\n <td>Dip direction, dip angle, measurement location<\/td>\n <\/tr>\n <tr>\n <td>Study area extent<\/td>\n <td>Text description<\/td>\n <td><span class=\"tag required\">Required<\/span><\/td>\n <td>X\/Y\/Z min and max values<\/td>\n <\/tr>\n <tr>\n <td>Stratigraphic sequence<\/td>\n <td>Text \/ Diagram<\/td>\n <td><span class=\"tag required\">Required<\/span><\/td>\n <td>Age relations, unconformities, intrusions<\/td>\n <\/tr>\n <tr>\n <td>Fault data<\/td>\n <td>Excel \/ CSV<\/td>\n <td><span class=\"tag optional\">Optional<\/span><\/td>\n <td>Fault plane orientation, displacement, extent<\/td>\n <\/tr>\n <tr>\n <td>Topography data<\/td>\n <td>DEM \/ GeoTIFF<\/td>\n <td><span class=\"tag optional\">Optional<\/span><\/td>\n <td>Digital elevation model of the study area<\/td>\n <\/tr>\n <tr>\n <td>Geological cross-sections<\/td>\n <td>Image \/ CAD<\/td>\n <td><span class=\"tag optional\">Optional<\/span><\/td>\n <td>For model validation<\/td>\n <\/tr>\n <\/table>\n <\/div>\n <\/section>\n\n <section>\n <h2>6. Official Documentation & Help Links<\/h2>\n <div class=\"link-list\">\n <a href=\"https:\/\/www.gempy.org\/\" target=\"_blank\" rel=\"noopener\">GemPy Official Website — Project homepage and features<\/a>\n <a href=\"https:\/\/docs.gempy.org\/\" target=\"_blank\" rel=\"noopener\">GemPy Documentation — Full API reference and tutorials<\/a>\n <a href=\"https:\/\/docs.gempy.org\/installation.html\" target=\"_blank\" rel=\"noopener\">Installation Guide — pip \/ conda \/ manual installation<\/a>\n <a href=\"https:\/\/docs.gempy.org\/tutorials\/a_getting_started\/get_started.html\" target=\"_blank\" rel=\"noopener\">Getting Started — Build your first model from scratch<\/a>\n <a href=\"https:\/\/docs.gempy.org\/tutorials\/ch1_fundamentals\/ch1_1_basics.html\" target=\"_blank\" rel=\"noopener\">Fundamentals Tutorial — CSV import and structural sequences<\/a>\n <a href=\"https:\/\/github.com\/cgre-aachen\/gempy\" target=\"_blank\" rel=\"noopener\">GitHub Repository — Source code, Issues, and example data<\/a>\n <a href=\"https:\/\/www.youtube.com\/@GemPy3D\" target=\"_blank\" rel=\"noopener\">YouTube Video Tutorials — Step-by-step demonstrations<\/a>\n <a href=\"https:\/\/www.geosci-model-dev.net\/12\/1\/2019\/gmd-12-1-2019.pdf\" target=\"_blank\" rel=\"noopener\">GemPy 1.0 Scientific Paper — Methodology and background<\/a>\n <\/div>\n <\/section>\n <\/div>\n\n <footer>\n <p>Guide prepared based on GemPy v3 official documentation (2024). For the latest updates, visit <a href=\"https:\/\/www.gempy.org\/\" target=\"_blank\">gempy.org<\/a>.<\/p>\n <\/footer>\n <\/div>\n\n <script>\n function setLang(lang) {\n document.querySelectorAll('.lang-section').forEach(el => el.classList.remove('active'));\n document.getElementById('lang-' + lang).classList.add('active');\n document.querySelectorAll('.lang-switch button').forEach(btn => btn.classList.remove('active'));\n document.getElementById('btn-' + lang).classList.add('active');\n document.documentElement.lang = lang === 'zh' ? 'zh-CN' : 'en';\n }\n\n function toggleTheme() {\n const html = document.documentElement;\n if (html.hasAttribute('data-theme')) {\n html.removeAttribute('data-theme');\n localStorage.removeItem('theme');\n } else {\n html.setAttribute('data-theme', 'dark');\n localStorage.setItem('theme', 'dark');\n }\n }\n\n \/\/ Restore theme preference\n if (localStorage.getItem('theme') === 'dark') {\n document.documentElement.setAttribute('data-theme', 'dark');\n }\n <\/script>\n<\/body>\n<\/html>\n\n","protected":false},"excerpt":{"rendered":"<p>GemPy \u5730\u8d28\u5efa\u6a21\u6307\u5357 | GemPy Geological Modeling Guide GemPy Gu […]<\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"open","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[11,9],"tags":[],"class_list":["post-297","post","type-post","status-publish","format-standard","hentry","category-cad","category-ai"],"blocksy_meta":[],"_links":{"self":[{"href":"https:\/\/numsimlab.com\/index.php?rest_route=\/wp\/v2\/posts\/297","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/numsimlab.com\/index.php?rest_route=\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/numsimlab.com\/index.php?rest_route=\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/numsimlab.com\/index.php?rest_route=\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/numsimlab.com\/index.php?rest_route=%2Fwp%2Fv2%2Fcomments&post=297"}],"version-history":[{"count":1,"href":"https:\/\/numsimlab.com\/index.php?rest_route=\/wp\/v2\/posts\/297\/revisions"}],"predecessor-version":[{"id":298,"href":"https:\/\/numsimlab.com\/index.php?rest_route=\/wp\/v2\/posts\/297\/revisions\/298"}],"wp:attachment":[{"href":"https:\/\/numsimlab.com\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=297"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/numsimlab.com\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=297"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/numsimlab.com\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=297"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}

{"id":297,"date":"2026-07-03T17:54:41","date_gmt":"2026-07-03T09:54:41","guid":{"rendered":"https:\/\/numsimlab.com\/?p=297"},"modified":"2026-07-03T17:54:41","modified_gmt":"2026-07-03T09:54:41","slug":"gempy-%e5%9c%b0%e8%b4%a8%e5%bb%ba%e6%a8%a1%e6%95%b0%e6%8d%ae%e5%87%86%e5%a4%87%e6%8c%87%e5%8d%97","status":"publish","type":"post","link":"https:\/\/numsimlab.com\/?p=297","title":{"rendered":"GemPy \u5730\u8d28\u5efa\u6a21\u6570\u636e\u51c6\u5907\u6307\u5357"},"content":{"rendered":"\n\n\n\n\n \n \n GemPy \u5730\u8d28\u5efa\u6a21\u6307\u5357 | GemPy Geological Modeling Guide<\/title>\n <style>\n :root {\n --bg: #f8f9fa;\n --surface: #ffffff;\n --ink: #1a1a2e;\n --muted: #5a6578;\n --rule: #dee2e6;\n --accent: #2563eb;\n --accent2: #059669;\n --code-bg: #1e293b;\n --code-text: #e2e8f0;\n --warn-bg: #fff7ed;\n --warn-border: #f97316;\n --info-bg: #eff6ff;\n --info-border: #3b82f6;\n }\n\n [data-theme=\"dark\"] {\n --bg: #0f172a;\n --surface: #1e293b;\n --ink: #f1f5f9;\n --muted: #94a3b8;\n --rule: #334155;\n --accent: #60a5fa;\n --accent2: #34d399;\n --code-bg: #020617;\n --code-text: #cbd5e1;\n --warn-bg: #2a1b0f;\n --warn-border: #f97316;\n --info-bg: #172554;\n --info-border: #3b82f6;\n }\n\n * { margin: 0; padding: 0; box-sizing: border-box; }\n\n body {\n font-family: -apple-system, BlinkMacSystemFont, \"Segoe UI\", Roboto, \"Helvetica Neue\", Arial, \"Noto Sans SC\", \"PingFang SC\", \"Microsoft YaHei\", sans-serif;\n background: var(--bg);\n color: var(--ink);\n line-height: 1.7;\n transition: background 0.3s, color 0.3s;\n }\n\n .container {\n max-width: 960px;\n margin: 0 auto;\n padding: 0 20px;\n }\n\n \/* Header *\/\n header {\n background: var(--surface);\n border-bottom: 1px solid var(--rule);\n padding: 1rem 0;\n position: sticky;\n top: 0;\n z-index: 100;\n }\n\n .header-inner {\n display: flex;\n justify-content: space-between;\n align-items: center;\n }\n\n .logo {\n font-size: 1.25rem;\n font-weight: 700;\n color: var(--accent);\n display: flex;\n align-items: center;\n gap: 0.5rem;\n }\n\n .logo svg {\n width: 28px;\n height: 28px;\n }\n\n .controls {\n display: flex;\n gap: 0.75rem;\n align-items: center;\n }\n\n .lang-switch {\n display: flex;\n background: var(--bg);\n border: 1px solid var(--rule);\n border-radius: 6px;\n overflow: hidden;\n }\n\n .lang-switch button {\n border: none;\n background: transparent;\n color: var(--muted);\n padding: 0.4rem 0.9rem;\n cursor: pointer;\n font-size: 0.85rem;\n transition: all 0.2s;\n }\n\n .lang-switch button.active {\n background: var(--accent);\n color: white;\n }\n\n .theme-toggle {\n background: var(--bg);\n border: 1px solid var(--rule);\n border-radius: 6px;\n padding: 0.4rem 0.7rem;\n cursor: pointer;\n color: var(--muted);\n font-size: 1rem;\n line-height: 1;\n }\n\n \/* Hero *\/\n .hero {\n padding: 3rem 0 2rem;\n text-align: center;\n }\n\n .hero h1 {\n font-size: 2rem;\n margin-bottom: 0.75rem;\n line-height: 1.3;\n }\n\n .hero p {\n color: var(--muted);\n font-size: 1.05rem;\n max-width: 640px;\n margin: 0 auto;\n }\n\n .badge {\n display: inline-block;\n background: var(--accent);\n color: white;\n padding: 0.25rem 0.75rem;\n border-radius: 999px;\n font-size: 0.75rem;\n font-weight: 600;\n margin-bottom: 1rem;\n text-transform: uppercase;\n letter-spacing: 0.05em;\n }\n\n \/* Sections *\/\n section {\n background: var(--surface);\n border: 1px solid var(--rule);\n border-radius: 10px;\n padding: 1.75rem;\n margin-bottom: 1.5rem;\n }\n\n section h2 {\n font-size: 1.35rem;\n margin-bottom: 1rem;\n padding-bottom: 0.5rem;\n border-bottom: 2px solid var(--accent);\n display: inline-block;\n }\n\n section h3 {\n font-size: 1.1rem;\n margin: 1.5rem 0 0.75rem;\n color: var(--accent);\n }\n\n section p, section li {\n margin-bottom: 0.75rem;\n color: var(--muted);\n }\n\n ul, ol {\n padding-left: 1.5rem;\n margin-bottom: 1rem;\n }\n\n li { margin-bottom: 0.4rem; }\n\n \/* Tables *\/\n .table-wrap {\n overflow-x: auto;\n margin: 1rem 0;\n }\n\n table {\n width: 100%;\n border-collapse: collapse;\n font-size: 0.9rem;\n min-width: 500px;\n }\n\n th, td {\n padding: 0.6rem 0.8rem;\n text-align: left;\n border-bottom: 1px solid var(--rule);\n }\n\n th {\n background: var(--bg);\n font-weight: 600;\n color: var(--ink);\n }\n\n tr:hover td { background: var(--bg); }\n\n \/* Code blocks *\/\n pre {\n background: var(--code-bg);\n color: var(--code-text);\n padding: 1.25rem;\n border-radius: 8px;\n overflow-x: auto;\n font-family: \"SF Mono\", Monaco, \"Cascadia Code\", Consolas, monospace;\n font-size: 0.85rem;\n line-height: 1.6;\n margin: 1rem 0;\n }\n\n code {\n background: var(--bg);\n color: var(--accent);\n padding: 0.15rem 0.4rem;\n border-radius: 4px;\n font-family: \"SF Mono\", Monaco, \"Cascadia Code\", Consolas, monospace;\n font-size: 0.88em;\n }\n\n pre code {\n background: transparent;\n color: inherit;\n padding: 0;\n }\n\n \/* Callouts *\/\n .callout {\n padding: 1rem 1.25rem;\n border-radius: 8px;\n margin: 1rem 0;\n border-left: 4px solid;\n }\n\n .callout.info {\n background: var(--info-bg);\n border-color: var(--info-border);\n }\n\n .callout.warn {\n background: var(--warn-bg);\n border-color: var(--warn-border);\n }\n\n .callout strong {\n display: block;\n margin-bottom: 0.3rem;\n color: var(--ink);\n }\n\n \/* Links *\/\n a {\n color: var(--accent);\n text-decoration: none;\n }\n\n a:hover { text-decoration: underline; }\n\n .link-list a {\n display: block;\n padding: 0.6rem 0;\n border-bottom: 1px solid var(--rule);\n }\n\n .link-list a:last-child { border-bottom: none; }\n\n \/* Tags *\/\n .tag {\n display: inline-block;\n background: var(--bg);\n border: 1px solid var(--rule);\n padding: 0.2rem 0.6rem;\n border-radius: 4px;\n font-size: 0.8rem;\n color: var(--muted);\n margin-right: 0.4rem;\n margin-bottom: 0.4rem;\n }\n\n .tag.required {\n background: var(--accent);\n color: white;\n border-color: var(--accent);\n }\n\n .tag.optional {\n background: var(--accent2);\n color: white;\n border-color: var(--accent2);\n }\n\n \/* Footer *\/\n footer {\n text-align: center;\n padding: 2rem 0;\n color: var(--muted);\n font-size: 0.85rem;\n }\n\n \/* Responsive *\/\n @media (max-width: 600px) {\n .hero h1 { font-size: 1.5rem; }\n section { padding: 1.25rem; }\n .header-inner { flex-wrap: wrap; gap: 0.5rem; }\n }\n\n \/* Hidden language sections *\/\n .lang-section { display: none; }\n .lang-section.active { display: block; }\n <\/style>\n<\/head>\n<body>\n <header>\n <div class=\"container header-inner\">\n <div class=\"logo\">\n <svg viewBox=\"0 0 24 24\" fill=\"none\" stroke=\"currentColor\" stroke-width=\"2\" stroke-linecap=\"round\" stroke-linejoin=\"round\">\n <path d=\"M12 2L2 7l10 5 10-5-10-5z\"\/>\n <path d=\"M2 17l10 5 10-5\"\/>\n <path d=\"M2 12l10 5 10-5\"\/>\n <\/svg>\n GemPy Guide\n <\/div>\n <div class=\"controls\">\n <div class=\"lang-switch\">\n <button class=\"active\" onclick=\"setLang('zh')\" id=\"btn-zh\">\u4e2d\u6587<\/button>\n <button onclick=\"setLang('en')\" id=\"btn-en\">English<\/button>\n <\/div>\n <button class=\"theme-toggle\" onclick=\"toggleTheme()\" title=\"Toggle theme\">☼<\/button>\n <\/div>\n <\/div>\n <\/header>\n\n <div class=\"container\">\n \n <div class=\"lang-section active\" id=\"lang-zh\">\n <div class=\"hero\">\n <span class=\"badge\">GemPy v3 (2024)<\/span>\n <h1>GemPy \u5730\u8d28\u5efa\u6a21\u6570\u636e\u51c6\u5907\u6307\u5357<\/h1>\n <p>\u57fa\u4e8e\u5f00\u6e90\u9690\u5f0f\u4e09\u7ef4\u5730\u8d28\u5efa\u6a21\u5e93 GemPy\uff0c\u4ecb\u7ecd\u4ece\u94bb\u5b54\u67f1\u72b6\u56fe\u6570\u636e\u5230 3D \u5730\u8d28\u6a21\u578b\u7684\u5b8c\u6574\u6570\u636e\u51c6\u5907\u6d41\u7a0b\u4e0e\u57fa\u672c\u5efa\u6a21\u793a\u4f8b\u3002<\/p>\n <\/div>\n\n <section>\n <h2>\u4e00\u3001GemPy \u7b80\u4ecb<\/h2>\n <p><a href=\"https:\/\/www.gempy.org\/\" target=\"_blank\">GemPy<\/a> \u662f\u4e00\u4e2a\u57fa\u4e8e Python \u7684\u5f00\u6e90\u4e09\u7ef4\u7ed3\u6784\u5730\u8d28\u5efa\u6a21\u5e93\uff0c\u91c7\u7528<strong>\u9690\u5f0f\u5efa\u6a21\uff08Implicit Modeling\uff09<\/strong>\u65b9\u6cd5\uff0c\u901a\u8fc7\u901a\u7528\u534f\u514b\u91cc\u91d1\u63d2\u503c\u6784\u5efa 3D \u6807\u91cf\u573a\u3002\u5b83\u80fd\u591f\u9ad8\u6548\u5730\u521b\u5efa\u5305\u542b\u8936\u76b1\u3001\u65ad\u5c42\u3001\u4e0d\u6574\u5408\u7b49\u590d\u6742\u6784\u9020\u7279\u5f81\u7684\u5730\u8d28\u6a21\u578b\uff0c\u5e76\u539f\u751f\u652f\u6301\u6982\u7387\u5efa\u6a21\u4e0e\u4e0d\u786e\u5b9a\u6027\u5206\u6790\u3002<\/p>\n <div class=\"callout info\">\n <strong>\u6838\u5fc3\u7279\u70b9<\/strong>\n GPU \u52a0\u901f\u8ba1\u7b97 | \u652f\u6301 PyTorch\/NumPy \u540e\u7aef | 2D\/3D \u53ef\u89c6\u5316\uff08PyVista\uff09| \u53ef\u76f4\u63a5\u5bfc\u51fa\u81f3 Blender \/ LiquidEarth\n <\/div>\n <\/section>\n\n <section>\n <h2>\u4e8c\u3001\u5efa\u6a21\u6240\u9700\u7684\u6700\u5c0f\u6570\u636e\u96c6<\/h2>\n <p>\u5728 GemPy \u4e2d\u6784\u5efa\u5730\u8d28\u6a21\u578b\uff0c\u9700\u8981\u4e24\u7c7b\u6838\u5fc3\u8f93\u5165\u6570\u636e\u3002\u8fd9\u4e9b\u6570\u636e\u53ef\u4ee5\u4ece\u94bb\u5b54\u67f1\u72b6\u56fe\u3001\u5730\u8d28\u5256\u9762\u56fe\u3001\u91ce\u5916\u9732\u5934\u6d4b\u91cf\u7b49\u8d44\u6599\u4e2d\u63d0\u53d6\u3002<\/p>\n\n <h3>1. \u754c\u9762\u70b9\uff08Surface Points\uff09<\/h3>\n <p>\u754c\u9762\u70b9\u662f\u4e09\u7ef4\u7a7a\u95f4\u4e2d\u6807\u8bb0\u5730\u8d28\u754c\u9762\uff08\u5982\u5730\u5c42\u5e95\u754c\u9762\u3001\u65ad\u5c42\u9762\u3001\u4e0d\u6574\u5408\u9762\uff09\u4f4d\u7f6e\u7684\u5750\u6807\u70b9\u3002<strong>\u5728 GemPy \u4e2d\uff0c\u6bcf\u4e2a\u754c\u9762\u4ee3\u8868\u5176\u6240\u8d4b\u5355\u5143\u7684\u5e95\u754c\u9762\u3002<\/strong><\/p>\n <div class=\"table-wrap\">\n <table>\n <tr>\n <th>\u5b57\u6bb5\u540d<\/th>\n <th>\u7c7b\u578b<\/th>\n <th>\u8bf4\u660e<\/th>\n <th>\u5fc5\u586b<\/th>\n <\/tr>\n <tr>\n <td><code>X<\/code><\/td>\n <td>\u6570\u503c<\/td>\n <td>\u4e1c\u5411\u5750\u6807\uff08\u4e0e\u5efa\u6a21\u8303\u56f4\u5355\u4f4d\u4e00\u81f4\uff09<\/td>\n <td><span class=\"tag required\">\u662f<\/span><\/td>\n <\/tr>\n <tr>\n <td><code>Y<\/code><\/td>\n <td>\u6570\u503c<\/td>\n <td>\u5317\u5411\u5750\u6807<\/td>\n <td><span class=\"tag required\">\u662f<\/span><\/td>\n <\/tr>\n <tr>\n <td><code>Z<\/code><\/td>\n <td>\u6570\u503c<\/td>\n <td>\u9ad8\u7a0b\/\u6df1\u5ea6\uff08\u901a\u5e38\u5730\u4e0b\u4e3a\u8d1f\u503c\uff09<\/td>\n <td><span class=\"tag required\">\u662f<\/span><\/td>\n <\/tr>\n <tr>\n <td><code>formation<\/code><\/td>\n <td>\u5b57\u7b26\u4e32<\/td>\n <td>\u5730\u5c42\/\u6784\u9020\u5355\u5143\u540d\u79f0\uff08\u5982 “Sandstone_1″\uff09<\/td>\n <td><span class=\"tag required\">\u662f<\/span><\/td>\n <\/tr>\n <\/table>\n <\/div>\n\n <h3>2. \u4ea7\u72b6\u6d4b\u91cf\uff08Orientations\uff09<\/h3>\n <p>\u4ea7\u72b6\u6570\u636e\u63cf\u8ff0\u5730\u8d28\u754c\u9762\u7684\u7a7a\u95f4\u503e\u659c\u65b9\u5411\uff0c\u7528\u4e8e\u7ea6\u675f\u63d2\u503c\u7684\u51e0\u4f55\u5f62\u6001\u3002\u6bcf\u4e2a\u6784\u9020\u7ec4\u81f3\u5c11\u9700\u8981\u4e00\u6761\u4ea7\u72b6\u6570\u636e\u3002<\/p>\n <div class=\"table-wrap\">\n <table>\n <tr>\n <th>\u5b57\u6bb5\u540d<\/th>\n <th>\u7c7b\u578b<\/th>\n <th>\u8bf4\u660e<\/th>\n <th>\u5fc5\u586b<\/th>\n <\/tr>\n <tr>\n <td><code>X<\/code>, <code>Y<\/code>, <code>Z<\/code><\/td>\n <td>\u6570\u503c<\/td>\n <td>\u6d4b\u91cf\u70b9\u7684\u4e09\u7ef4\u5750\u6807<\/td>\n <td><span class=\"tag required\">\u662f<\/span><\/td>\n <\/tr>\n <tr>\n <td><code>azimuth<\/code><\/td>\n <td>\u6570\u503c<\/td>\n <td>\u503e\u5411\uff080\u2013360\u5ea6\uff09<\/td>\n <td><span class=\"tag required\">\u662f<\/span><\/td>\n <\/tr>\n <tr>\n <td><code>dip<\/code><\/td>\n <td>\u6570\u503c<\/td>\n <td>\u503e\u89d2\uff080\u201390\u5ea6\uff09<\/td>\n <td><span class=\"tag required\">\u662f<\/span><\/td>\n <\/tr>\n <tr>\n <td><code>polarity<\/code><\/td>\n <td>\u6570\u503c<\/td>\n <td>\u6781\u6027\uff08\u901a\u5e38\u53d6 1\uff0c\u5012\u7f6e\u6784\u9020\u53ef\u53d6 -1\uff09<\/td>\n <td><span class=\"tag required\">\u662f<\/span><\/td>\n <\/tr>\n <tr>\n <td><code>formation<\/code><\/td>\n <td>\u5b57\u7b26\u4e32<\/td>\n <td>\u6240\u5c5e\u5730\u5c42\/\u6784\u9020\u5355\u5143\u540d\u79f0<\/td>\n <td><span class=\"tag required\">\u662f<\/span><\/td>\n <\/tr>\n <\/table>\n <\/div>\n\n <div class=\"callout warn\">\n <strong>\u6700\u5c0f\u6570\u636e\u89c4\u5219<\/strong>\n \u6bcf\u4e2a\u6784\u9020\u7ec4\uff08Series\/Stack\uff09\u4e2d\uff0c\u81f3\u5c11\u9700\u8981\u4e00\u4e2a\u754c\u9762\u7684<strong>\u4e24\u4e2a\u754c\u9762\u70b9<\/strong> + <strong>\u4e00\u6761\u4ea7\u72b6<\/strong>\u3002\u6dfb\u52a0\u66f4\u591a\u540c\u7ec4\u754c\u9762\u65f6\uff0c\u6bcf\u4e2a\u65b0\u754c\u9762\u53ea\u9700\u81f3\u5c11\u4e00\u4e2a\u754c\u9762\u70b9\u5373\u53ef\u5171\u4eab\u8be5\u7ec4\u7684\u4ea7\u72b6\u4fe1\u606f\u3002\n <\/div>\n <\/section>\n\n <section>\n <h2>\u4e09\u3001\u4ece\u94bb\u5b54\u67f1\u72b6\u56fe\u63d0\u53d6\u6570\u636e<\/h2>\n <p>\u94bb\u5b54\u67f1\u72b6\u56fe\u662f\u5730\u8d28\u5efa\u6a21\u6700\u5e38\u89c1\u7684\u539f\u59cb\u8d44\u6599\u3002\u4ee5\u4e0b\u662f\u4ece\u94bb\u5b54\u67f1\u72b6\u56fe\u8f6c\u6362\u4e3a GemPy \u8f93\u5165\u6570\u636e\u7684\u5b8c\u6574\u6d41\u7a0b\uff1a<\/p>\n\n <h3>\u6b65\u9aa4 1\uff1a\u6536\u96c6\u539f\u59cb\u94bb\u5b54\u8d44\u6599<\/h3>\n <p>\u7528\u6237\u9700\u8981\u63d0\u4f9b\u4ee5\u4e0b\u6750\u6599\uff1a<\/p>\n <ol>\n <li><strong>\u94bb\u5b54\u57fa\u672c\u4fe1\u606f\u8868<\/strong>\uff1a\u5b54\u53f7\u3001\u5b54\u53e3\u5750\u6807\uff08X, Y\uff09\u3001\u5b54\u53e3\u9ad8\u7a0b\uff08Z\uff09\u3001\u94bb\u5b54\u6df1\u5ea6\u3001\u94bb\u5b54\u503e\u89d2\u4e0e\u65b9\u4f4d\u89d2\uff08\u82e5\u4e3a\u659c\u5b54\uff09<\/li>\n <li><strong>\u94bb\u5b54\u5206\u5c42\u6570\u636e\u8868<\/strong>\uff1a\u6bcf\u5c42\u7684\u5c42\u5e95\u6df1\u5ea6\u3001\u5c42\u5e95\u9ad8\u7a0b\u3001\u5730\u5c42\u4ee3\u53f7\u3001\u5730\u5c42\u540d\u79f0\u3001\u5ca9\u6027\u63cf\u8ff0<\/li>\n <li><strong>\u5730\u5c42\u4ea7\u72b6\u6570\u636e<\/strong>\uff1a\u94bb\u5b54\u5185\u6216\u94bb\u5b54\u9644\u8fd1\u6d4b\u5f97\u7684\u5730\u5c42\u503e\u5411\u4e0e\u503e\u89d2<\/li>\n <li><strong>\u7814\u7a76\u533a\u8303\u56f4<\/strong>\uff1a\u5efa\u6a21\u7684\u5e73\u9762\u8303\u56f4\uff08X_min, X_max, Y_min, Y_max\uff09\u4e0e\u9ad8\u7a0b\u8303\u56f4\uff08Z_min, Z_max\uff09<\/li>\n <li><strong>\u5730\u5c42\u5e8f\u5217\u5173\u7cfb<\/strong>\uff1a\u5404\u5730\u5c42\u7684\u6c89\u79ef\u987a\u5e8f\u3001\u662f\u5426\u5b58\u5728\u4e0d\u6574\u5408\/\u4fb5\u5165\u5173\u7cfb\u3001\u65ad\u5c42\u4fe1\u606f<\/li>\n <\/ol>\n\n <h3>\u6b65\u9aa4 2\uff1a\u8ba1\u7b97\u754c\u9762\u70b9\u5750\u6807<\/h3>\n <p>\u5bf9\u6bcf\u53e3\u94bb\u5b54\u7684\u6bcf\u4e2a\u5730\u5c42\u5e95\u754c\u9762\uff0c\u8ba1\u7b97\u5176\u4e09\u7ef4\u5750\u6807\uff1a<\/p>\n <pre><code># \u76f4\u5b54\u793a\u4f8b\uff1a\u5b54\u53e3\u5750\u6807 (X0, Y0, Z0)\uff0c\u5c42\u5e95\u6df1\u5ea6 D\nX = X0\nY = Y0\nZ = Z0 - D # \u82e5 Z \u5411\u4e0b\u4e3a\u8d1f\n\n# \u659c\u5b54\u9700\u6839\u636e\u5929\u9876\u89d2 theta \u548c\u65b9\u4f4d\u89d2 alpha \u8fdb\u884c\u6362\u7b97\n# X = X0 + D * sin(theta) * sin(alpha)\n# Y = Y0 + D * sin(theta) * cos(alpha)\n# Z = Z0 - D * cos(theta)<\/code><\/pre>\n\n <h3>\u6b65\u9aa4 3\uff1a\u51c6\u5907 CSV \u6587\u4ef6<\/h3>\n <p>\u5c06\u63d0\u53d6\u7684\u754c\u9762\u70b9\u4fdd\u5b58\u4e3a <code>surface_points.csv<\/code>\uff1a<\/p>\n <pre><code>X,Y,Z,formation\n0,0,800,rock3\n0,0,450,rock2\n0,0,250,rock1\n500,0,800,rock3\n...<\/code><\/pre>\n\n <p>\u5c06\u4ea7\u72b6\u6570\u636e\u4fdd\u5b58\u4e3a <code>orientations.csv<\/code>\uff1a<\/p>\n <pre><code>X,Y,Z,azimuth,dip,polarity,formation\n100,500,650,90,315,1,rock2\n100,500,450,90,315,1,rock1\n...<\/code><\/pre>\n\n <div class=\"callout info\">\n <strong>\u5750\u6807\u7cfb\u63d0\u793a<\/strong>\n GemPy \u4f7f\u7528\u53f3\u624b\u5750\u6807\u7cfb\u3002Z \u8f74\u53ef\u5411\u4e0a\u4e3a\u6b63\uff08\u5730\u8868\u9ad8\u7a0b\uff09\u6216\u5411\u4e0b\u4e3a\u8d1f\uff08\u6df1\u5ea6\uff09\uff0c\u53ea\u9700\u5728\u6574\u4e2a\u6a21\u578b\u4e2d\u4fdd\u6301\u4e00\u81f4\u3002\u5efa\u8bae\u4e0e\u7814\u7a76\u533a\u5e38\u7528\u5750\u6807\u7cfb\u5bf9\u9f50\u3002\n <\/div>\n <\/section>\n\n <section>\n <h2>\u56db\u3001GemPy \u5efa\u6a21\u57fa\u672c\u6d41\u7a0b\uff08Demo\uff09<\/h2>\n\n <h3>1. \u5b89\u88c5 GemPy<\/h3>\n <pre><code># \u6838\u5fc3\u5e93\npip install gempy\n\n# \u53ef\u89c6\u5316\u6a21\u5757\npip install gempy-viewer\n\n# 3D \u53ef\u89c6\u5316\u4f9d\u8d56\uff08\u63a8\u8350\u5b89\u88c5\uff09\npip install pyvista\n\n# \u82e5\u9700 GPU \u52a0\u901f\uff08PyTorch \u540e\u7aef\uff09\npip install torch<\/code><\/pre>\n\n <h3>2. \u4ece CSV \u5bfc\u5165\u6570\u636e\u5e76\u521b\u5efa\u6a21\u578b<\/h3>\n <pre><code>import gempy as gp\nimport gempy_viewer as gpv\n\n# \u6570\u636e\u8def\u5f84\ndata_path = \"https:\/\/raw.githubusercontent.com\/cgre-aachen\/gempy_data\/master\/\"\npath_to_data = data_path + \"\/data\/input_data\/video_tutorials_v3\/\"\n\n# \u521b\u5efa\u5730\u8d28\u6a21\u578b\ngeo_model = gp.create_geomodel(\n project_name='tutorial_model',\n extent=[0, 2500, 0, 1000, 0, 1000],\n resolution=[100, 40, 40],\n importer_helper=gp.data.ImporterHelper(\n path_to_orientations=path_to_data + \"tutorial_model_orientations.csv\",\n path_to_surface_points=path_to_data + \"tutorial_model_surface_points.csv\"\n )\n)<\/code><\/pre>\n\n <h3>3. \u5b9a\u4e49\u6784\u9020\u5e8f\u5217\u4e0e\u6620\u5c04\u5173\u7cfb<\/h3>\n <pre><code># \u6620\u5c04\u5730\u5c42\u5230\u6784\u9020\u7ec4\uff08\u4ece\u5e74\u8f7b\u5230\u5e74\u8001\u6392\u5e8f\uff09\ngp.map_stack_to_surfaces(\n gempy_model=geo_model,\n mapping_object={\n \"Strat_Series2\": (\"rock3\"),\n \"Strat_Series1\": (\"rock2\", \"rock1\")\n }\n)<\/code><\/pre>\n\n <h3>4. \u8ba1\u7b97\u6a21\u578b<\/h3>\n <pre><code># \u8ba1\u7b97\u5730\u8d28\u6a21\u578b\nimport os\nos.environ[\"DEFAULT_BACKEND\"] = \"PYTORCH\" # \u6216 NUMPY\n\ngp.compute_model(geo_model)<\/code><\/pre>\n\n <h3>5. \u53ef\u89c6\u5316\u7ed3\u679c<\/h3>\n <pre><code># 2D \u5256\u9762\u53ef\u89c6\u5316\ngpv.plot_2d(geo_model, cell_number=20, show_data=True)\n\n# 3D \u4ea4\u4e92\u5f0f\u53ef\u89c6\u5316\ngpv.plot_3d(geo_model, show_lith=True, show_boundaries=True)<\/code><\/pre>\n\n <h3>6. \u4ece\u96f6\u624b\u52a8\u521b\u5efa\uff08\u4e0d\u4f9d\u8d56 CSV\uff09<\/h3>\n <pre><code># \u521b\u5efa\u7a7a\u6a21\u578b\ngeo_model = gp.create_geomodel(\n project_name='Model1',\n extent=[0, 780, -200, 200, -582, 0],\n resolution=(50, 50, 50),\n structural_frame=gp.data.StructuralFrame.initialize_default_structure()\n)\n\n# \u624b\u52a8\u6dfb\u52a0\u754c\u9762\u70b9\ngp.add_surface_points(\n geo_model=geo_model,\n x=[225, 460, 617],\n y=[0, 0, 0],\n z=[-95, -100, -10],\n elements_names=['Limestone', 'Limestone', 'Limestone']\n)\n\n# \u624b\u52a8\u6dfb\u52a0\u4ea7\u72b6\ngp.add_orientations(\n geo_model=geo_model,\n x=[350], y=[0], z=[-120],\n elements_names=['Limestone'],\n pole_vector=[[0, 0, 1]] # \u6c34\u5e73\u5c42\u7406\uff0c\u6cd5\u5411\u91cf\u671d\u4e0a\n)\n\n# \u8ba1\u7b97\u5e76\u53ef\u89c6\u5316\ngp.compute_model(geo_model)\ngpv.plot_2d(geo_model)<\/code><\/pre>\n <\/section>\n\n <section>\n <h2>\u4e94\u3001\u5b8c\u6574\u6570\u636e\u9700\u6c42\u6e05\u5355<\/h2>\n <p>\u4e3a\u4f7f\u7528 GemPy \u6210\u529f\u5efa\u7acb\u4e09\u7ef4\u5730\u8d28\u6a21\u578b\uff0c\u8bf7\u786e\u4fdd\u63d0\u4f9b\u4ee5\u4e0b\u8d44\u6599\uff1a<\/p>\n <div class=\"table-wrap\">\n <table>\n <tr>\n <th>\u6750\u6599<\/th>\n <th>\u683c\u5f0f<\/th>\n <th>\u5fc5\u8981\u6027<\/th>\n <th>\u8bf4\u660e<\/th>\n <\/tr>\n <tr>\n <td>\u94bb\u5b54\u5750\u6807\u4e0e\u5b54\u53e3\u9ad8\u7a0b<\/td>\n <td>Excel \/ CSV<\/td>\n <td><span class=\"tag required\">\u5fc5\u9700<\/span><\/td>\n <td>\u6bcf\u53e3\u5b54\u7684 X, Y, Z \u5750\u6807<\/td>\n <\/tr>\n <tr>\n <td>\u94bb\u5b54\u5206\u5c42\u6570\u636e<\/td>\n <td>Excel \/ CSV<\/td>\n <td><span class=\"tag required\">\u5fc5\u9700<\/span><\/td>\n <td>\u6bcf\u5c42\u5e95\u6df1\u3001\u5730\u5c42\u540d\u79f0\u3001\u5ca9\u6027<\/td>\n <\/tr>\n <tr>\n <td>\u5730\u5c42\u4ea7\u72b6\u6d4b\u91cf<\/td>\n <td>Excel \/ CSV<\/td>\n <td><span class=\"tag required\">\u5fc5\u9700<\/span><\/td>\n <td>\u503e\u5411\u3001\u503e\u89d2\u53ca\u6d4b\u91cf\u4f4d\u7f6e\u5750\u6807<\/td>\n <\/tr>\n <tr>\n <td>\u7814\u7a76\u533a\u8303\u56f4<\/td>\n <td>\u6587\u672c\u8bf4\u660e<\/td>\n <td><span class=\"tag required\">\u5fc5\u9700<\/span><\/td>\n <td>X\/Y\/Z \u7684\u6700\u5c0f\u6700\u5927\u503c<\/td>\n <\/tr>\n <tr>\n <td>\u5730\u5c42\u5e8f\u5217\u5173\u7cfb<\/td>\n <td>\u6587\u672c\/\u56fe\u8868<\/td>\n <td><span class=\"tag required\">\u5fc5\u9700<\/span><\/td>\n <td>\u65b0\u8001\u5173\u7cfb\u3001\u4e0d\u6574\u5408\u3001\u4fb5\u5165\u4f53<\/td>\n <\/tr>\n <tr>\n <td>\u65ad\u5c42\u6570\u636e<\/td>\n <td>Excel \/ CSV<\/td>\n <td><span class=\"tag optional\">\u53ef\u9009<\/span><\/td>\n <td>\u65ad\u5c42\u9762\u4ea7\u72b6\u3001\u65ad\u8ddd\u3001\u5ef6\u4f38\u8303\u56f4<\/td>\n <\/tr>\n <tr>\n <td>\u5730\u5f62\u6570\u636e<\/td>\n <td>DEM \/ GeoTIFF<\/td>\n <td><span class=\"tag optional\">\u53ef\u9009<\/span><\/td>\n <td>\u7814\u7a76\u533a\u6570\u5b57\u9ad8\u7a0b\u6a21\u578b<\/td>\n <\/tr>\n <tr>\n <td>\u5730\u8d28\u5256\u9762\u56fe<\/td>\n <td>\u56fe\u7247 \/ CAD<\/td>\n <td><span class=\"tag optional\">\u53ef\u9009<\/span><\/td>\n <td>\u8f85\u52a9\u9a8c\u8bc1\u6a21\u578b\u5408\u7406\u6027<\/td>\n <\/tr>\n <\/table>\n <\/div>\n <\/section>\n\n <section>\n <h2>\u516d\u3001\u5b98\u65b9\u6587\u6863\u4e0e\u5e2e\u52a9\u94fe\u63a5<\/h2>\n <div class=\"link-list\">\n <a href=\"https:\/\/www.gempy.org\/\" target=\"_blank\" rel=\"noopener\">GemPy \u5b98\u65b9\u7f51\u7ad9 — \u9879\u76ee\u4e3b\u9875\u4e0e\u529f\u80fd\u4ecb\u7ecd<\/a>\n <a href=\"https:\/\/docs.gempy.org\/\" target=\"_blank\" rel=\"noopener\">GemPy \u5b98\u65b9\u6587\u6863 — \u5b8c\u6574 API \u4e0e\u6559\u7a0b<\/a>\n <a href=\"https:\/\/docs.gempy.org\/installation.html\" target=\"_blank\" rel=\"noopener\">\u5b89\u88c5\u6307\u5357 — pip \/ conda \/ \u624b\u52a8\u5b89\u88c5<\/a>\n <a href=\"https:\/\/docs.gempy.org\/tutorials\/a_getting_started\/get_started.html\" target=\"_blank\" rel=\"noopener\">\u5165\u95e8\u6559\u7a0b — \u4ece\u96f6\u5f00\u59cb\u6784\u5efa\u7b2c\u4e00\u4e2a\u6a21\u578b<\/a>\n <a href=\"https:\/\/docs.gempy.org\/tutorials\/ch1_fundamentals\/ch1_1_basics.html\" target=\"_blank\" rel=\"noopener\">\u57fa\u7840\u5efa\u6a21\u6559\u7a0b — CSV \u6570\u636e\u5bfc\u5165\u4e0e\u6784\u9020\u5e8f\u5217<\/a>\n <a href=\"https:\/\/github.com\/cgre-aachen\/gempy\" target=\"_blank\" rel=\"noopener\">GitHub \u4ed3\u5e93 — \u6e90\u7801\u3001Issues \u4e0e\u793a\u4f8b\u6570\u636e<\/a>\n <a href=\"https:\/\/www.youtube.com\/@GemPy3D\" target=\"_blank\" rel=\"noopener\">YouTube \u89c6\u9891\u6559\u7a0b — \u5206\u6b65\u64cd\u4f5c\u6f14\u793a<\/a>\n <a href=\"https:\/\/www.geosci-model-dev.net\/12\/1\/2019\/gmd-12-1-2019.pdf\" target=\"_blank\" rel=\"noopener\">GemPy 1.0 \u5b66\u672f\u8bba\u6587 — \u65b9\u6cd5\u8bba\u4e0e\u79d1\u5b66\u80cc\u666f<\/a>\n <\/div>\n <\/section>\n <\/div>\n\n \n <div class=\"lang-section\" id=\"lang-en\">\n <div class=\"hero\">\n <span class=\"badge\">GemPy v3 (2024)<\/span>\n <h1>GemPy Geological Modeling Data Preparation Guide<\/h1>\n <p>A comprehensive guide to preparing borehole log data for 3D structural geological modeling using the open-source implicit modeling library GemPy.<\/p>\n <\/div>\n\n <section>\n <h2>1. Introduction to GemPy<\/h2>\n <p><a href=\"https:\/\/www.gempy.org\/\" target=\"_blank\">GemPy<\/a> is an open-source Python library for 3D structural geological modeling. It uses an <strong>implicit modeling<\/strong> approach based on universal cokriging to construct 3D scalar fields, enabling fast creation of complex models including folds, faults, and unconformities. It also natively supports probabilistic modeling and uncertainty analysis.<\/p>\n <div class=\"callout info\">\n <strong>Key Features<\/strong>\n GPU-accelerated computation | PyTorch\/NumPy backend | 2D\/3D visualization via PyVista | Export to Blender \/ LiquidEarth\n <\/div>\n <\/section>\n\n <section>\n <h2>2. Minimum Required Datasets<\/h2>\n <p>To build a geological model in GemPy, two core input data types are required. These can be extracted from borehole logs, geological cross-sections, and outcrop measurements.<\/p>\n\n <h3>2.1 Surface Points<\/h3>\n <p>Surface points are 3D coordinates marking geological interfaces (e.g., layer boundaries, fault planes, unconformities). <strong>In GemPy, each surface represents the bottom of the unit it is assigned to.<\/strong><\/p>\n <div class=\"table-wrap\">\n <table>\n <tr>\n <th>Field<\/th>\n <th>Type<\/th>\n <th>Description<\/th>\n <th>Required<\/th>\n <\/tr>\n <tr>\n <td><code>X<\/code><\/td>\n <td>Numeric<\/td>\n <td>Easting coordinate<\/td>\n <td><span class=\"tag required\">Yes<\/span><\/td>\n <\/tr>\n <tr>\n <td><code>Y<\/code><\/td>\n <td>Numeric<\/td>\n <td>Northing coordinate<\/td>\n <td><span class=\"tag required\">Yes<\/span><\/td>\n <\/tr>\n <tr>\n <td><code>Z<\/code><\/td>\n <td>Numeric<\/td>\n <td>Elevation \/ depth (negative for subsurface)<\/td>\n <td><span class=\"tag required\">Yes<\/span><\/td>\n <\/tr>\n <tr>\n <td><code>formation<\/code><\/td>\n <td>String<\/td>\n <td>Formation \/ structural element name<\/td>\n <td><span class=\"tag required\">Yes<\/span><\/td>\n <\/tr>\n <\/table>\n <\/div>\n\n <h3>2.2 Orientation Measurements<\/h3>\n <p>Orientation data describes the spatial dip of geological interfaces, constraining the geometry of interpolation. At least one orientation per structural group is required.<\/p>\n <div class=\"table-wrap\">\n <table>\n <tr>\n <th>Field<\/th>\n <th>Type<\/th>\n <th>Description<\/th>\n <th>Required<\/th>\n <\/tr>\n <tr>\n <td><code>X<\/code>, <code>Y<\/code>, <code>Z<\/code><\/td>\n <td>Numeric<\/td>\n <td>Measurement location coordinates<\/td>\n <td><span class=\"tag required\">Yes<\/span><\/td>\n <\/tr>\n <tr>\n <td><code>azimuth<\/code><\/td>\n <td>Numeric<\/td>\n <td>Dip direction (0\u2013360 degrees)<\/td>\n <td><span class=\"tag required\">Yes<\/span><\/td>\n <\/tr>\n <tr>\n <td><code>dip<\/code><\/td>\n <td>Numeric<\/td>\n <td>Dip angle (0\u201390 degrees)<\/td>\n <td><span class=\"tag required\">Yes<\/span><\/td>\n <\/tr>\n <tr>\n <td><code>polarity<\/code><\/td>\n <td>Numeric<\/td>\n <td>Polarity (usually 1; -1 for inverted structures)<\/td>\n <td><span class=\"tag required\">Yes<\/span><\/td>\n <\/tr>\n <tr>\n <td><code>formation<\/code><\/td>\n <td>String<\/td>\n <td>Associated formation \/ element name<\/td>\n <td><span class=\"tag required\">Yes<\/span><\/td>\n <\/tr>\n <\/table>\n <\/div>\n\n <div class=\"callout warn\">\n <strong>Minimum Data Rule<\/strong>\n Each structural group (Series\/Stack) requires at least <strong>two surface points<\/strong> for one surface + <strong>one orientation<\/strong>. Additional surfaces in the same group only need one surface point each, as they share orientation information.\n <\/div>\n <\/section>\n\n <section>\n <h2>3. Extracting Data from Borehole Logs<\/h2>\n <p>Borehole logs are the most common source of data for geological modeling. Below is the complete workflow to convert borehole data into GemPy inputs:<\/p>\n\n <h3>Step 1: Collect Raw Borehole Materials<\/h3>\n <p>The user should provide the following:<\/p>\n <ol>\n <li><strong>Borehole location table<\/strong>: Borehole ID, collar coordinates (X, Y), collar elevation (Z), total depth, inclination and azimuth (for deviated holes)<\/li>\n <li><strong>Stratigraphic interval table<\/strong>: Bottom depth of each layer, bottom elevation, formation code, formation name, lithology description<\/li>\n <li><strong>Structural orientation data<\/strong>: Dip direction and dip angle measured within or near the borehole<\/li>\n <li><strong>Study area extent<\/strong>: X\/Y\/Z minimum and maximum values defining the modeling domain<\/li>\n <li><strong>Stratigraphic sequence<\/strong>: Depositional order, presence of unconformities or intrusions, fault information<\/li>\n <\/ol>\n\n <h3>Step 2: Compute Surface Point Coordinates<\/h3>\n <p>For each layer bottom in each borehole, calculate the 3D coordinates:<\/p>\n <pre><code># Vertical borehole example: collar (X0, Y0, Z0), layer bottom depth D\nX = X0\nY = Y0\nZ = Z0 - D # if Z is negative downward\n\n# Deviated borehole conversion using zenith angle theta and azimuth alpha:\n# X = X0 + D * sin(theta) * sin(alpha)\n# Y = Y0 + D * sin(theta) * cos(alpha)\n# Z = Z0 - D * cos(theta)<\/code><\/pre>\n\n <h3>Step 3: Prepare CSV Files<\/h3>\n <p>Save extracted surface points as <code>surface_points.csv<\/code>:<\/p>\n <pre><code>X,Y,Z,formation\n0,0,800,rock3\n0,0,450,rock2\n0,0,250,rock1\n500,0,800,rock3\n...<\/code><\/pre>\n\n <p>Save orientation data as <code>orientations.csv<\/code>:<\/p>\n <pre><code>X,Y,Z,azimuth,dip,polarity,formation\n100,500,650,90,315,1,rock2\n100,500,450,90,315,1,rock1\n...<\/code><\/pre>\n\n <div class=\"callout info\">\n <strong>Coordinate System Note<\/strong>\n GemPy uses a right-handed coordinate system. The Z axis can be positive upward (elevation) or negative downward (depth); consistency across the model is what matters. Align with your study area’s standard coordinate system.\n <\/div>\n <\/section>\n\n <section>\n <h2>4. Basic GemPy Modeling Workflow (Demo)<\/h2>\n\n <h3>4.1 Install GemPy<\/h3>\n <pre><code># Core library\npip install gempy\n\n# Visualization module\npip install gempy-viewer\n\n# 3D visualization dependency (recommended)\npip install pyvista\n\n# For GPU acceleration (PyTorch backend)\npip install torch<\/code><\/pre>\n\n <h3>4.2 Import Data from CSV and Create Model<\/h3>\n <pre><code>import gempy as gp\nimport gempy_viewer as gpv\n\n# Data paths\ndata_path = \"https:\/\/raw.githubusercontent.com\/cgre-aachen\/gempy_data\/master\/\"\npath_to_data = data_path + \"\/data\/input_data\/video_tutorials_v3\/\"\n\n# Create geological model\ngeo_model = gp.create_geomodel(\n project_name='tutorial_model',\n extent=[0, 2500, 0, 1000, 0, 1000],\n resolution=[100, 40, 40],\n importer_helper=gp.data.ImporterHelper(\n path_to_orientations=path_to_data + \"tutorial_model_orientations.csv\",\n path_to_surface_points=path_to_data + \"tutorial_model_surface_points.csv\"\n )\n)<\/code><\/pre>\n\n <h3>4.3 Define Structural Sequence<\/h3>\n <pre><code># Map formations to structural groups (youngest to oldest)\ngp.map_stack_to_surfaces(\n gempy_model=geo_model,\n mapping_object={\n \"Strat_Series2\": (\"rock3\"),\n \"Strat_Series1\": (\"rock2\", \"rock1\")\n }\n)<\/code><\/pre>\n\n <h3>4.4 Compute the Model<\/h3>\n <pre><code>import os\nos.environ[\"DEFAULT_BACKEND\"] = \"PYTORCH\" # or \"NUMPY\"\n\ngp.compute_model(geo_model)<\/code><\/pre>\n\n <h3>4.5 Visualize Results<\/h3>\n <pre><code># 2D cross-section\ngpv.plot_2d(geo_model, cell_number=20, show_data=True)\n\n# 3D interactive visualization\ngpv.plot_3d(geo_model, show_lith=True, show_boundaries=True)<\/code><\/pre>\n\n <h3>4.6 Create from Scratch (Without CSV)<\/h3>\n <pre><code># Create empty model\ngeo_model = gp.create_geomodel(\n project_name='Model1',\n extent=[0, 780, -200, 200, -582, 0],\n resolution=(50, 50, 50),\n structural_frame=gp.data.StructuralFrame.initialize_default_structure()\n)\n\n# Add surface points manually\ngp.add_surface_points(\n geo_model=geo_model,\n x=[225, 460, 617],\n y=[0, 0, 0],\n z=[-95, -100, -10],\n elements_names=['Limestone', 'Limestone', 'Limestone']\n)\n\n# Add orientation manually\ngp.add_orientations(\n geo_model=geo_model,\n x=[350], y=[0], z=[-120],\n elements_names=['Limestone'],\n pole_vector=[[0, 0, 1]] # Horizontal bedding, normal vector up\n)\n\n# Compute and visualize\ngp.compute_model(geo_model)\ngpv.plot_2d(geo_model)<\/code><\/pre>\n <\/section>\n\n <section>\n <h2>5. Complete Data Requirements Checklist<\/h2>\n <p>To successfully build a 3D geological model with GemPy, please ensure the following materials are provided:<\/p>\n <div class=\"table-wrap\">\n <table>\n <tr>\n <th>Material<\/th>\n <th>Format<\/th>\n <th>Required<\/th>\n <th>Notes<\/th>\n <\/tr>\n <tr>\n <td>Borehole coordinates & collar elevation<\/td>\n <td>Excel \/ CSV<\/td>\n <td><span class=\"tag required\">Required<\/span><\/td>\n <td>X, Y, Z for each borehole<\/td>\n <\/tr>\n <tr>\n <td>Stratigraphic interval data<\/td>\n <td>Excel \/ CSV<\/td>\n <td><span class=\"tag required\">Required<\/span><\/td>\n <td>Bottom depth, formation name, lithology<\/td>\n <\/tr>\n <tr>\n <td>Structural orientation measurements<\/td>\n <td>Excel \/ CSV<\/td>\n <td><span class=\"tag required\">Required<\/span><\/td>\n <td>Dip direction, dip angle, measurement location<\/td>\n <\/tr>\n <tr>\n <td>Study area extent<\/td>\n <td>Text description<\/td>\n <td><span class=\"tag required\">Required<\/span><\/td>\n <td>X\/Y\/Z min and max values<\/td>\n <\/tr>\n <tr>\n <td>Stratigraphic sequence<\/td>\n <td>Text \/ Diagram<\/td>\n <td><span class=\"tag required\">Required<\/span><\/td>\n <td>Age relations, unconformities, intrusions<\/td>\n <\/tr>\n <tr>\n <td>Fault data<\/td>\n <td>Excel \/ CSV<\/td>\n <td><span class=\"tag optional\">Optional<\/span><\/td>\n <td>Fault plane orientation, displacement, extent<\/td>\n <\/tr>\n <tr>\n <td>Topography data<\/td>\n <td>DEM \/ GeoTIFF<\/td>\n <td><span class=\"tag optional\">Optional<\/span><\/td>\n <td>Digital elevation model of the study area<\/td>\n <\/tr>\n <tr>\n <td>Geological cross-sections<\/td>\n <td>Image \/ CAD<\/td>\n <td><span class=\"tag optional\">Optional<\/span><\/td>\n <td>For model validation<\/td>\n <\/tr>\n <\/table>\n <\/div>\n <\/section>\n\n <section>\n <h2>6. Official Documentation & Help Links<\/h2>\n <div class=\"link-list\">\n <a href=\"https:\/\/www.gempy.org\/\" target=\"_blank\" rel=\"noopener\">GemPy Official Website — Project homepage and features<\/a>\n <a href=\"https:\/\/docs.gempy.org\/\" target=\"_blank\" rel=\"noopener\">GemPy Documentation — Full API reference and tutorials<\/a>\n <a href=\"https:\/\/docs.gempy.org\/installation.html\" target=\"_blank\" rel=\"noopener\">Installation Guide — pip \/ conda \/ manual installation<\/a>\n <a href=\"https:\/\/docs.gempy.org\/tutorials\/a_getting_started\/get_started.html\" target=\"_blank\" rel=\"noopener\">Getting Started — Build your first model from scratch<\/a>\n <a href=\"https:\/\/docs.gempy.org\/tutorials\/ch1_fundamentals\/ch1_1_basics.html\" target=\"_blank\" rel=\"noopener\">Fundamentals Tutorial — CSV import and structural sequences<\/a>\n <a href=\"https:\/\/github.com\/cgre-aachen\/gempy\" target=\"_blank\" rel=\"noopener\">GitHub Repository — Source code, Issues, and example data<\/a>\n <a href=\"https:\/\/www.youtube.com\/@GemPy3D\" target=\"_blank\" rel=\"noopener\">YouTube Video Tutorials — Step-by-step demonstrations<\/a>\n <a href=\"https:\/\/www.geosci-model-dev.net\/12\/1\/2019\/gmd-12-1-2019.pdf\" target=\"_blank\" rel=\"noopener\">GemPy 1.0 Scientific Paper — Methodology and background<\/a>\n <\/div>\n <\/section>\n <\/div>\n\n <footer>\n <p>Guide prepared based on GemPy v3 official documentation (2024). For the latest updates, visit <a href=\"https:\/\/www.gempy.org\/\" target=\"_blank\">gempy.org<\/a>.<\/p>\n <\/footer>\n <\/div>\n\n <script>\n function setLang(lang) {\n document.querySelectorAll('.lang-section').forEach(el => el.classList.remove('active'));\n document.getElementById('lang-' + lang).classList.add('active');\n document.querySelectorAll('.lang-switch button').forEach(btn => btn.classList.remove('active'));\n document.getElementById('btn-' + lang).classList.add('active');\n document.documentElement.lang = lang === 'zh' ? 'zh-CN' : 'en';\n }\n\n function toggleTheme() {\n const html = document.documentElement;\n if (html.hasAttribute('data-theme')) {\n html.removeAttribute('data-theme');\n localStorage.removeItem('theme');\n } else {\n html.setAttribute('data-theme', 'dark');\n localStorage.setItem('theme', 'dark');\n }\n }\n\n \/\/ Restore theme preference\n if (localStorage.getItem('theme') === 'dark') {\n document.documentElement.setAttribute('data-theme', 'dark');\n }\n <\/script>\n<\/body>\n<\/html>\n\n","protected":false},"excerpt":{"rendered":"<p>GemPy \u5730\u8d28\u5efa\u6a21\u6307\u5357 | GemPy Geological Modeling Guide GemPy Gu […]<\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"open","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[11,9],"tags":[],"class_list":["post-297","post","type-post","status-publish","format-standard","hentry","category-cad","category-ai"],"blocksy_meta":[],"_links":{"self":[{"href":"https:\/\/numsimlab.com\/index.php?rest_route=\/wp\/v2\/posts\/297","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/numsimlab.com\/index.php?rest_route=\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/numsimlab.com\/index.php?rest_route=\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/numsimlab.com\/index.php?rest_route=\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/numsimlab.com\/index.php?rest_route=%2Fwp%2Fv2%2Fcomments&post=297"}],"version-history":[{"count":1,"href":"https:\/\/numsimlab.com\/index.php?rest_route=\/wp\/v2\/posts\/297\/revisions"}],"predecessor-version":[{"id":298,"href":"https:\/\/numsimlab.com\/index.php?rest_route=\/wp\/v2\/posts\/297\/revisions\/298"}],"wp:attachment":[{"href":"https:\/\/numsimlab.com\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=297"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/numsimlab.com\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=297"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/numsimlab.com\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=297"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}