Giới thiệu dự án

Sản xuất và chế biến chè là ngành kinh tế kỹ thuật quan trọng trong cơ cấu nông nghiệp Việt Nam, quốc gia hiện đứng thứ 7 toàn cầu về sản lượng và thứ 5 về xuất khẩu chè với tổng diện tích canh tác đạt 115.942 ha và sản lượng trên 1,09 triệu tấn búp tươi (Bộ NN&PTNT; FAOSTAT). Tại vùng kinh tế trọng điểm phía Nam, xã Phú Hội (huyện Nhơn Trạch, tỉnh Đồng Nai) là địa bàn bản địa duy nhất bảo tồn và phát triển giống chè truyền thống hàng trăm năm tuổi kết hợp với nguồn nước khoáng tự nhiên "Mạch Bà".

+-------------------------------------------------------------------------+
|                  BỐI CẢNH NÔNG HỘ TRỒNG CHÈ PHÚ HỘI                     |
|                                                                         |
|  [Đô thị hóa & KCN] -----> [Thu hẹp quỹ đất: 50 ha -> 8.5 ha]           |
|  [Chuyển dịch lao động] -> [Thiếu nhân lực trẻ & Kỹ thuật canh tác]    |
|  [Áp lực thị trường] ----> [Hàng giả/pha trộn làm xói mòn uy tín OCOP] |
+-------------------------------------------------------------------------+

1. Vấn đề nghiên cứu và các điểm nghẽn thực tế

Quá trình công nghiệp hóa và đô thị hóa nhanh chóng tại Nhơn Trạch (với 9 khu công nghiệp tập trung) đã tạo ra sự suy giảm nghiêm trọng về diện tích đất canh tác chè từ 50 ha (trước năm 1975) xuống còn 7 ha vào năm 2016 và phục hồi chậm lên 8,5 ha năm 2022. Các nông hộ đối mặt với 3 điểm nghẽn cốt lõi:

  • Xói mòn nguồn gen bản địa: Suy giảm diện tích đe dọa trực tiếp đến quần thể cây chè cổ thụ có tuổi đời hàng thế kỷ.
  • Rủi ro suy thoái nguồn tài nguyên sinh thái: Khai thác quá mức nguồn nước ngầm Mạch Bà làm thay đổi phẩm cấp lý hóa đặc trưng của sản phẩm trà hoàn thiện.
  • Bất đối xứng thông tin và thị trường: Tình trạng pha trộn chè kém chất lượng từ các vùng khác làm suy giảm giá trị kinh tế của thương hiệu chè đặc sản Phú Hội.

2. Mục tiêu dự án

  1. Định lượng hóa hiện trạng: Khảo sát, đo lường năng suất, diện tích và hiệu quả tài chính của các nông hộ sản xuất chè giai đoạn 2016–2022.
  2. Xác định các yếu tố quyết định: Xây dựng mô hình kinh tế lượng xác định các yếu tố ảnh hưởng có ý nghĩa thống kê đến quyết định bảo tồn ($P < 0,05$).
  3. Thiết kế khung chính sách: Đề xuất hệ thống giải pháp kỹ thuật canh tác hữu cơ, liên kết chuỗi giá trị và bảo hộ nhãn hiệu tập thể.

3. Phương pháp tiếp cận và kết quả kỳ vọng

Dự án áp dụng mô hình xác suất hồi quy Binary Logistic trên tập dữ liệu khảo sát thực chứng từ 62 nông hộ (37 hộ bảo tồn, 25 hộ không bảo tồn).

  • Mục tiêu kỹ thuật: Đạt độ chính xác dự báo (Overall Percentage Correct) $\ge 85%$, chỉ số Nagelkerke $R^2 \ge 0,55$, kiểm định độ phù hợp Hosmer-Lemeshow ($p > 0,05$).
  • Mục tiêu kinh tế: Cung cấp cơ sở định lượng để nâng năng suất chè búp khô từ $3,0\text{ tạ/ha}$ lên $4,5\text{ tạ/ha}$ và tối ưu hóa doanh thu mô hình xen canh đạt trên 150 triệu VNĐ/ha/năm.
  • Phạm vi & Giới hạn: Nghiên cứu tập trung tại xã Phú Hội, huyện Nhơn Trạch, tỉnh Đồng Nai; dữ liệu khảo sát trực tiếp vụ mùa 2020–2022.

Phân tích và thiết kế giải pháp

Phân tích hiện trạng

Phương thức canh tác Ưu điểm cốt lõi Nhược điểm & Rủi ro Chi phí đầu tư (VNĐ/1000m²/năm) Hiệu quả kinh tế tương đối
Độc canh chè truyền thống Dễ quản lý kỹ thuật đơn lẻ, thu hoạch tập trung. Dễ phát sinh dịch sâu bệnh hại (rầy xanh, bọ cánh tơ), thu nhập bấp bênh khi biến động giá. $8.500.000 - 12.000.000$ Thấp (chỉ đạt $18 - 25$ triệu VNĐ lãi thuần).
Xen canh cây ăn trái (Khuyến nghị) Tận dụng bóng râm điều hòa vi khí hậu, giữ ẩm đất feralit, đa dạng hóa nguồn thu. Cần kỹ thuật tỉa cành tạo tán phức tạp, đòi hỏi kiểm soát phân bón hữu cơ gián tiếp. $14.000.000 - 18.000.000$ Cao (lãi thuần chè kết hợp bưởi/sầu riêng $> 65$ triệu VNĐ).
Chuyển đổi hoàn toàn sang cây ăn trái Thu nhập tức thời cao từ cây đặc sản. Mất vĩnh viễn nguồn gen bản địa quý hiếm, triệt tiêu làng nghề truyền thống. $20.000.000 - 30.000.000$ Phụ thuộc hoàn toàn vào chu kỳ sâu bệnh cây ăn trái.

Ma trận yêu cầu giải pháp theo khung MoSCoW

  • Must-have (Bắt buộc): Mô hình xác suất nhị phân kiểm định các biến nhân khẩu học, kinh tế, khuyến nông; Kiểm soát hiện tượng đa cộng tuyến với $\text{VIF} < 10$.
  • Should-have (Cần có): Thiết kế quy trình canh tác hữu cơ theo tiêu chuẩn VietGAP kết hợp mô hình kinh tế tập thể (HTX/Tổ hợp tác).
  • Could-have (Có thể có): Xây dựng hệ thống truy xuất nguồn gốc QR code tích hợp chỉ dẫn địa lý OCOP Đồng Nai.
  • Won't-have (Không thực hiện): Khảo sát các vùng trồng chè ngoài địa bàn Nhơn Trạch trong khuôn khổ đề tài này.

Thiết kế hệ thống

Technology Stack & Chuẩn kỹ thuật

  • Nền tảng tính toán thống kê: IBM SPSS Statistics v26.0 (Phân tích hồi quy tham số nhị phân) và Python v3.10 (Thư viện statsmodels 0.14.0, scikit-learn 1.2.2).
  • Nền tảng quản trị dữ liệu: Microsoft Excel 2021 với Schema dữ liệu được chuẩn hóa kiểu biến số thực nghiệm.
  • Tiêu chuẩn an toàn & Đạo đức dữ liệu: Mã hóa danh tính chủ hộ (Anonymization), tuân thủ nguyên tắc chấp thuận tham gia nghiên cứu (Informed Consent).
+-----------------------------------------------------------------------------------+
|               DATA DICTIONARY: BỘ BIẾN SỐ ĐẦU VÀO MÔ HÌNH HỒI QUY                 |
+------+----------------------+------------+--------------+-------------------------+
| Biến | Ý nghĩa              | Phân loại  | Thang đo     | Kỳ vọng tương quan      |
+------+----------------------+------------+--------------+-------------------------+
| Y    | Quyết định bảo tồn   | Phụ thuộc  | Binary (0/1) | Biến đích               |
| X1   | Tuổi của chủ hộ      | Độc lập    | Tỉ lệ (Năm)  | (+) Dương               |
| X2   | Tổng diện tích đất   | Độc lập    | Tỉ lệ (1000m²| (+) Dương               |
| X3   | Kinh nghiệm canh tác | Độc lập    | Tỉ lệ (Năm)  | (+) Dương               |
| X4   | Tham gia khuyến nông | Độc lập    | Dummy (0/1)  | (+) Dương               |
| X5   | Trình độ học vấn     | Độc lập    | Thứ bậc (Năm)| (+) Dương               |
| X6   | Kỹ thuật canh tác    | Độc lập    | Likert (1-3) | (+) Dương               |
| X7   | Mức thu nhập từ chè  | Độc lập    | Likert (1-5) | (+) Dương               |
+------+----------------------+------------+--------------+-------------------------+

Methodology

Nghiên cứu áp dụng quy trình thực nghiệm định lượng theo phương pháp luận kinh tế lượng chuẩn mực:

  • Cỡ mẫu nghiên cứu ($N$): Xác định theo công thức ước lượng của Harris (1985): $n = 50 + p = 50 + 7 = 57 \le 62$ quan sát thực tế, đảm bảo bậc tự do và công suất kiểm định (Statistical Power $> 0,8$).
  • Phân tầng mẫu: 37 hộ duy trì bảo tồn chè ($59,68%$) và 25 hộ dừng canh tác ($40,32%$).
  • Đánh giá rủi ro mô hình: Phòng tránh lỗi phân loại sai loại I ($\alpha = 0,05$) và loại II ($\beta = 0,2$) thông qua kiểm định ma trận phân loại ngẫu nhiên.

Implementation và kết quả

Development Process & Thuật toán

Mô hình xác suất bảo tồn chè Phú Hội tuân theo hàm phân phối Logistic lũy tích:

$$P_i = P(Y_i = 1 | X_i) = \frac{1}{1 + e^{-(\beta_0 + \sum_{k=1}^{7} \beta_k X_{ki})}}$$

Biến đổi Logit tương đương:

$$\ln\left(\frac{P_i}{1 - P_i}\right) = \beta_0 + \beta_1 X_1 + \beta_2 X_2 + \beta_3 X_3 + \beta_4 X_4 + \beta_5 X_5 + \beta_6 X_6 + \beta_7 X_7 + \varepsilon_i$$

Trong đó, tỷ số Odds Ratio $\text{Exp}(\beta_k)$ định lượng sự thay đổi cơ hội bảo tồn chè khi biến $X_k$ tăng 1 đơn vị.

import numpy as np
import pandas as pd
import statsmodels.api as sm
from sklearn.metrics import classification_report, confusion_matrix

# 1. Khởi tạo Dataset từ dữ liệu thực chứng 62 nông hộ Phú Hội
def fit_logistic_regression(df):
    features = ['X1_Age', 'X2_LandArea', 'X3_Exp', 'X4_Training', 'X5_Edu', 'X6_Tech', 'X7_Income']
    X = df[features]
    X = sm.add_constant(X)
    y = df['Y_Preserve']
    
    # 2. Ước lượng mô hình Logit bằng phương pháp Maximum Likelihood (MLE)
    logit_model = sm.Logit(y, X)
    result = logit_model.fit(method='newton', maxiter=100, disp=False)
    
    # 3. Tính toán Odds Ratio Exp(B) và 95% Confidence Interval
    params = result.params
    conf = result.conf_int()
    conf['Odds_Ratio'] = params
    conf.columns = ['2.5%', '97.5%', 'Odds_Ratio']
    odds_ratio = np.exp(conf)
    
    return result, odds_ratio

# Hiển thị cấu trúc tóm tắt ước lượng thực nghiệm
print("Logistic Regression Pipeline initialized successfully.")

Testing và validation

Kết quả kiểm định thống kê trên phần mềm SPSS v26 cho thấy mô hình đạt độ vững chắc cao:

  • Kiểm định Omnibus: Chi-square $\chi^2 = 48,621$ ($p = 0,000 < 0,001$), chứng minh tập hợp các biến độc lập cải thiện vượt bậc khả năng giải thích so với mô hình gốc chỉ có hằng số.
  • Mức độ giải thích: Chỉ số Cox & Snell $R^2 = 0,544$; Nagelkerke $R^2 = 0,736$ (mô hình giải thích được $73,6%$ sự biến thiên của quyết định bảo tồn).
  • Kiểm định Hosmer & Lemeshow: $\chi^2 = 4,812$ với $p = 0,777 > 0,05$, khẳng định mô hình có độ phù hợp xuất sắc với dữ liệu thực nghiệm.
  • Độ chính xác phân loại: Tỷ lệ dự báo chính xác tổng thể đạt $88,7%$ (Nhóm bảo tồn: $91,9%$; Nhóm không bảo tồn: $84,0%$).
+------------------------------------------------------------------------------------+
|               BẢNG KẾT QUẢ HỒI QUY BINARY LOGISTIC THỰC NGHIỆM                    |
+---------------------+---------+--------+----------+----+-------+---------+---------+
| Biến số             | B (Hệ số| S.E.   | Wald     | df | Sig.  | Exp(B)  | Ý nghĩa |
+---------------------+---------+--------+----------+----+-------+---------+---------+
| X1 (Tuổi chủ hộ)    | 0.114   | 0.048  | 5.642    | 1  | 0.018 | 1.121   | P < 0.05|
| X2 (Diện tích đất)  | 0.892   | 0.364  | 6.004    | 1  | 0.014 | 2.440   | P < 0.05|
| X3 (Kinh nghiệm)    | -0.042  | 0.051  | 0.678    | 1  | 0.410 | 0.959   | Không ý |
| X4 (Tập huấn NN)    | 1.845   | 0.762  | 5.861    | 1  | 0.015 | 6.328   | P < 0.05|
| X5 (Trình độ học vấn| 0.128   | 0.155  | 0.681    | 1  | 0.409 | 1.137   | Không ý |
| X6 (Kỹ thuật canh tá| 1.412   | 0.598  | 5.575    | 1  | 0.018 | 4.104   | P < 0.05|
| X7 (Thu nhập từ chè)| 2.105   | 0.685  | 9.444    | 1  | 0.002 | 8.207   | P < 0.01|
| Constant (Hằng số)  | -12.435 | 4.120  | 9.112    | 1  | 0.003 | 0.000   | P < 0.01|
+---------------------+---------+--------+----------+----+-------+---------+---------+

Kết quả đạt được

  1. 5 biến số có tác động tích cực mang ý nghĩa thống kê ($p < 0,05$):
    • Mức thu nhập từ chè ($X_7$): Có tác động mạnh nhất ($B = 2,105, \text{Exp}(B) = 8,207$). Nông hộ có mức đánh giá thu nhập tăng 1 bậc thì cơ hội bảo tồn chè tăng gấp $8,207$ lần.
    • Tham gia tập huấn nông nghiệp ($X_4$): Hệ số $B = 1,845, \text{Exp}(B) = 6,328$. Hộ được tập huấn khuyến nông có xác suất quyết định bảo tồn cao gấp $6,328$ lần hộ không tham gia.
    • Kỹ thuật trồng ($X_6$): $\text{Exp}(B) = 4,104$, nắm vững kỹ thuật giúp tăng hơn 4 lần xác suất duy trì vườn chè.
    • Tổng diện tích đất ($X_2$): $\text{Exp}(B) = 2,440$, diện tích đất càng lớn tạo điều kiện quy hoạch xen canh hiệu quả.
    • Độ tuổi chủ hộ ($X_1$): $\text{Exp}(B) = 1,121$, phản ánh sự gắn bó văn hóa truyền thống sâu sắc của các nông hộ cao tuổi.
  2. 2 biến số không có ý nghĩa thống kê ($p > 0,05$): Kinh nghiệm ($X_3, p = 0,410$) và Trình độ học vấn ($X_5, p = 0,409$). Điều này chứng minh rằng thói quen canh tác theo lối mòn không quyết định hành vi bảo tồn bằng việc tiếp cận kỹ thuật mới và lợi ích tài chính trực tiếp.

Đổi mới và đóng góp

So sánh đối chiếu với các công trình nghiên cứu tiền nhiệm

Tiêu chí phân tích Khóa luận này (Chè Phú Hội, 2023) Nguyễn Minh Hà & CS (2020) Nguyễn Văn Lạc & CS (2021)
Đối tượng & Địa bàn Chè đặc sản bản địa (Nhơn Trạch, Đồng Nai) Công nghệ tưới nhỏ giọt trên nho (Ninh Thuận) Canh tác rau an toàn (Thừa Thiên Huế)
Phương pháp phân tích Binary Logistic kết hợp đánh giá hiệu quả tài chính mô hình xen canh Binary Logistic (21 biến giải thích độc lập) Logit nhị phân kết hợp mẫu phân tầng (N=350)
Yếu tố chi phối lớn nhất Mức thu nhập ($\text{Exp}(B)=8,207$) & Khuyến nông ($\text{Exp}(B)=6,328$) Kỳ vọng năng suất & Hỗ trợ chính sách Nhà nước Mức độ hiểu biết quy trình VietGAP & Diện tích
Độ chính xác dự báo $88,7%$ (Nagelkerke $R^2 = 0,736$) $82,4%$ (Nagelkerke $R^2 = 0,612$) $84,9%$ (Nagelkerke $R^2 = 0,658$)
Đóng góp mới Chứng minh vai trò bảo trợ vi khí hậu của mô hình xen canh $100%$ đối với chè vùng nhiệt đới Đánh giá vai trò của công nghệ tự động hóa tưới Xác thực nhận thức tiêu chuẩn an toàn sinh học

Đóng góp học thuật và thực tiễn

  • Về mặt phương pháp: Xây dựng mô hình định lượng xác thực mối liên hệ giữa các biến số phi tài chính (tập huấn, nhận thức kỹ thuật) với hành vi bảo tồn nguồn gen nông nghiệp có nguy cơ mai một.
  • Về mặt chính sách: Cung cấp bằng chứng thực nghiệm để chính quyền huyện Nhơn Trạch triển khai chương trình liên kết bao tiêu qua Hộ kinh doanh Phúc Bảo và dự án OCOP tỉnh Đồng Nai.

Ứng dụng thực tế và triển khai

Phân tích hiệu quả kinh tế mô hình canh tác (Tính trên $1.000\text{ m}^2$/năm)

+-------------------------------------------------------------------------+
|                  CƠ CẤU CHI PHÍ VÀ DOANH THU THỰC TẾ                    |
|                                                                         |
|  [Tổng doanh thu: 54.000.000 VNĐ]                                       |
|    |-- Doanh thu búp khô sơ chế: 36.000.000 VNĐ (60 kg x 600.000 đ)     |
|    +-- Doanh thu cây ăn trái xen canh: 18.000.000 VNĐ                   |
|                                                                         |
|  [Tổng chi phí sản xuất: 16.500.000 VNĐ]                                |
|    |-- Phân bón hữu cơ & vi sinh: 5.200.000 VNĐ                         |
|    |-- Công lao động & Chế biến thủ công: 8.800.000 VNĐ                 |
|    +-- Khấu hao công cụ, nước tưới Mạch Bà: 2.500.000 VNĐ               |
|                                                                         |
|  ===> LỢI NHUẬN THUẦN: 37.500.000 VNĐ / 1.000 m² (Tương đương 375 tr/ha)|
|       Tỷ suất lợi nhuận trên chi phí (BCR): 3,27                        |
+-------------------------------------------------------------------------+

Lộ trình triển khai thực địa (Roadmap)

2023 Q1-Q2: [Chuẩn hóa vườn giống F1 & Lập bản đồ gen cây chè cổ thụ]
2023 Q3-Q4: [Mở 10 lớp tập huấn kỹ thuật ủ phân hữu cơ & Cắt tỉa tạo tán]
2024 Q1-Q2: [Thành lập Hợp tác xã Nông nghiệp Chè Phú Hội & Cấp mã số vùng trồng]
2024 Q3-Q4: [Đăng ký bảo hộ Chỉ dẫn địa lý (GI) & Kết nối chuỗi du lịch sinh thái miệt vườn]

Hạn chế và hướng phát triển

1. Hạn chế kỹ thuật và dữ liệu

  • Quy mô mẫu: Mẫu nghiên cứu ($N=62$) bị giới hạn bởi tổng số hộ còn duy trì hoạt động trồng chè thực tế tại xã Phú Hội do tốc độ đô thị hóa nhanh.
  • Tính chu kỳ: Dữ liệu tài chính chưa bóc tách chi tiết tác động của các biến động thời tiết cực đoan (nắng nóng kéo dài làm giảm lưu lượng nước Mạch Bà).

2. Định hướng nghiên cứu mở rộng

  • Tích hợp mô hình Structural Equation Modeling (SEM): Đánh giá sâu hơn các biến tâm lý học hành vi, ý thức gìn giữ di sản văn hóa bản địa.
  • Nghiên cứu hóa sinh chuyên sâu: Phân tích định lượng hàm lượng Tanin, Catechin, EGCG và các hợp chất bay hơi trong búp chè Phú Hội tưới bằng nước Mạch Bà so với các vùng thổ nhưỡng khác.

Đối tượng hưởng lợi

  • Sinh viên & Học viên cao học: Tiếp cận bộ khung phân tích Binary Logistic hoàn chỉnh, ứng dụng giải quyết bài toán kinh tế nông nghiệp thực tế.
  • Kỹ sư nông nghiệp & Nông dân: Nắm bắt cơ sở khoa học để áp dụng mô hình xen canh hữu cơ, loại bỏ thuốc BVTV hóa học nhằm tối ưu hóa năng suất búp khô.
  • Cơ quan quản lý địa phương (UBND Huyện Nhơn Trạch): Có luận cứ định lượng để ban hành cơ chế trợ giá, quy hoạch khu vực bảo tồn nguồn gen chè đặc sản 8,5 ha.
  • Doanh nghiệp & HTX chế biến: Hoạch định chuỗi liên kết giá trị bền vững, nâng giá trị thương phẩm trà ướp hương lên $500.000 - 600.000\text{ VNĐ/kg}$.

Câu hỏi thường gặp

1. Yêu cầu kỹ thuật tối thiểu để nông hộ duy trì mô hình bảo tồn chè Phú Hội là gì?

Nông hộ cần diện tích tối thiểu $500\text{ m}^2$, áp dụng tỷ lệ che bóng $30-40%$ thông qua trồng xen canh cây ăn quả, tuyệt đối không lạm dụng phân bón hóa học nhằm bảo vệ chất lượng nước ngầm Mạch Bà.

2. Tại sao biến kinh nghiệm canh tác lại không có ý nghĩa thống kê trong mô hình?

Kết quả kiểm định Wald ($p = 0,410 > 0,05$) chỉ ra rằng kinh nghiệm lâu năm theo phương pháp cũ không tạo ra sự khác biệt. Yếu tố then chốt thúc đẩy bảo tồn hiện nay là khả năng tiếp thu kỹ thuật mới qua các lớp khuyến nông ($p = 0,015$).

3. Làm thế nào để giải quyết xung đột giữa quy hoạch khu công nghiệp và diện tích trồng chè?

Cần xây dựng vùng đệm sinh thái bảo tồn nông nghiệp đô thị tại Phú Hội, kết hợp mô hình nông nghiệp trải nghiệm du lịch làng nghề để gia tăng giá trị kinh tế trên một đơn vị diện tích đất.

4. Chi phí đầu tư ban đầu cho 1 ha chè xen canh là bao nhiêu và thời gian hoàn vốn?

Chi phí kiến thiết cơ bản ban đầu khoảng $80 - 100$ triệu VNĐ/ha (cây giống, hệ thống tưới tiết kiệm). Với tỷ suất sinh lời hiện tại, thời gian hoàn vốn đầu tư đạt từ 2,5 – 3 năm.

5. Giải pháp nào ngăn chặn triệt để tình trạng chè giả thương hiệu Phú Hội?

Bắt buộc áp dụng tem chống hàng giả điện tử chứa mã QR truy xuất nguồn gốc do Hợp tác xã/Hộ kinh doanh Phúc Bảo quản lý, gắn liền với chứng nhận OCOP 3 sao - 4 sao của tỉnh Đồng Nai.


Kết luận

Khóa luận tốt nghiệp của tác giả Võ Thị Hoàng Anh đã giải quyết thấu đáo bài toán thực tiễn về bảo tồn cây chè đặc sản Phú Hội tại huyện Nhơn Trạch, tỉnh Đồng Nai. Bằng việc ứng dụng mô hình kinh tế lượng Binary Logistic với độ chính xác dự báo $88,7%$, nghiên cứu đã định lượng hóa rõ nét 5 nhân tố cốt lõi chi phối quyết định của nông hộ: Thu nhập từ chè ($\text{Exp}(B)=8,207$), Khuyến nông ($\text{Exp}(B)=6,328$), Kỹ thuật canh tác ($\text{Exp}(B)=4,104$), Diện tích đất ($\text{Exp}(B)=2,440$) và Tuổi chủ hộ ($\text{Exp}(B)=1,121$). Đây là nền tảng thực chứng vững chắc để địa phương phục hồi thương hiệu "Nước Mạch Bà, Trà Phú Hội", kết hợp hài hòa giữa phát triển kinh tế sinh thái và bảo tồn di sản nông nghiệp bền vững.