Giới thiệu dự án

Sản xuất rau quả tại Việt Nam đang ghi nhận tốc độ tăng trưởng vượt bậc với kim ngạch xuất khẩu đạt 3,45 tỷ USD năm 2017 (tăng 40,5% so với cùng kỳ năm 2016), vượt qua nhiều ngành hàng truyền thống như lúa gạo (2,2 tỷ USD). Tuy nhiên, ngành trồng trọt nội địa vẫn đối mặt với bài toán chất lượng và tính chuẩn hóa quy trình canh tác, đặc biệt là khi tiếp cận các tiêu chuẩn VietGAP và hàng rào kỹ thuật quốc tế. Trong nhóm cây thực phẩm giá trị cao, giống bí ngồi xanh lai F1 (Cucurbita pepo L. hybrid LS 188) có tiềm năng năng suất lớn (40–50 tấn/ha) và giàu giá trị sinh học (protein hạt đạt 40–48%, chứa 8 axit amin thiết yếu như tryptophan 1,6%, phenylalanine 5,7%, leucine 4,7%, valine 4,2%).

Mặc dù vậy, kỹ thuật thâm canh bí ngồi xanh tại vùng đồng bằng và trung du Bắc Bộ còn mang tính tự phát. Nông dân thường áp dụng chế độ bón phân mất cân đối, lạm dụng hoặc thiếu hụt phân đạm (Nitrogen), dẫn đến tình trạng cây sinh trưởng sinh dưỡng quá mức (lốp đổ), kéo dài thời gian ra hoa đậu quả, tăng tỷ lệ nhiễm nấm bệnh (Phytophthora, Erysiphe cichoracearum) hoặc làm giảm nghiêm trọng năng suất và chất lượng thương phẩm.

Mục tiêu nghiên cứu

  1. Đánh giá các chỉ tiêu sinh trưởng và phát triển: Theo dõi thời gian qua các giai đoạn (nảy mầm, ra hoa, thu hoạch), động thái tăng trưởng chiều cao thân và tốc độ phân hóa số lá của giống bí ngồi xanh LS 188 dưới các mức đạm khác nhau.
  2. Xác định các yếu tố cấu thành năng suất và năng suất thực thu (NSTT): Đo đạc số quả trung bình/cây, khối lượng trung bình quả (g), chiều dài (cm), đường kính quả (cm), năng suất cá thể (NSCT) và năng suất lý thuyết (NSLT).
  3. Tối ưu hóa hiệu quả kinh tế: Tính toán chi phí đầu vào, tổng thu và lãi thuần ($Q = A - (B + C)$) để đưa ra khuyến nghị liều lượng đạm Urea $\text{CO(NH}_2)_2$ kinh tế và bền vững nhất cho vùng sinh thái Xuân Mai, Chương Mỹ, Hà Nội.

Phạm vi và giới hạn nghiên cứu

  • Địa điểm: Vườn thực nghiệm sinh thái, Viện Quản lý đất đai và Phát triển nông thôn – Trường Đại học Lâm nghiệp (Xuân Mai, Chương Mỹ, Hà Nội).
  • Thời gian thực hiện: Vụ Đông Xuân 2017 – 2018 (tháng 12/2017 đến tháng 03/2018).
  • Đối tượng khảo nghiệm: Giống bí ngồi xanh F1 LS 188; nguồn phân đạm khảo nghiệm là đạm Urê tinh thể chứa 46% N nguyên chất.

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

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

Thực trạng canh tác bí ngồi tại miền Bắc cho thấy sự chênh lệch lớn giữa các phương thức bón phân. Việc so sánh giữa tập quán canh tác truyền thống, quy trình khuyến nông đại trà và phương pháp thâm canh định lượng chính xác được thể hiện dưới bảng sau:

Tiêu chí phân tích Tập quán nông hộ truyền thống Quy trình khuyến nông chuẩn (Đối chứng) Giải pháp tối ưu hóa dinh dưỡng (CT3)
Liều lượng Đạm (N) Tự phát (80–120 kg N/ha hoặc >250 kg N/ha) 139 kg N nguyên chất/ha (5 kg N/sào BB) 194 kg N nguyên chất/ha (7 kg N/sào BB)
Phân bổ giai đoạn 1–2 lần bón ồ ạt, dễ bay hơi/rửa trôi Bón lót 20%, thúc 3 lần (20% - 30% - 30%) Bón lót 20%, thúc 3 lần kết hợp vun xới, tưới ẩm
Nền phân chuồng & P-K Thiếu cân đối, ít ủ hoai mục 20 tấn phân chuồng + 60 kg $\text{P}_2\text{O}_5$ + 150 kg $\text{K}_2\text{O}$ 20 tấn chuồng hoai + 60 kg $\text{P}_2\text{O}_5$ + 150 kg $\text{K}_2\text{O}$
Khả năng kiểm soát bệnh Thấp, dễ cháy lá và thối cổ rễ Trung bình, cây phát triển ổn định Cao, tỷ lệ nhiễm bệnh thấp, lá dày khỏe
Năng suất thực thu 1.600 – 1.900 tạ/ha 2.297,00 tạ/ha 2.669,33 tạ/ha (+16,21% so với ĐC)

Phân tích yêu cầu nông học theo mô hình MoSCoW

  • Must have (Bắt buộc): Nền phân hữu cơ hoai mục 20.000 kg/ha; duy trì độ ẩm đất 70–75% thời kỳ ra hoa đậu quả; kỹ thuật thụ phấn bổ sung lúc 7:00 – 10:00 sáng.
  • Should have (Nên có): Chia tỷ lệ bón đạm theo 4 mốc sinh trưởng: Lót 20% + Thúc 1 (hồi xanh 20%) + Thúc 2 (4–5 lá thật 30%) + Thúc 3 (đậu quả non 30%).
  • Could have (Có thể): Bổ sung phân bón lá vi lượng (Bo, Zn, Mg) giai đoạn nở hoa rộ để tăng tỷ lệ đậu quả hữu hiệu.
  • Won't have (Loại bỏ): Tuyệt đối không bón đạm Urê đơn độc vượt quá 250 kg N/ha; không bón sát gốc làm xót rễ non.

Thiết kế hệ thống thí nghiệm

Thí nghiệm được thiết kế theo phương pháp khối ngẫu nhiên hoàn toàn (Completely Randomized Design - CRD) với 4 công thức liều lượng đạm, lặp lại 3 lần. Mật độ gieo trồng đạt 6.944 cây/ha (khoảng cách hàng $\times$ cây: $1{,}0\text{ m} \times 1{,}0\text{ m}$; tương đương 250 cây/sào Bắc Bộ $360\text{ m}^2$). Diện tích mỗi ô thí nghiệm là $5\text{ m}^2$ với 5 cây theo dõi/lần lặp ($N = 60$ cây mẫu).

Ma trận Công thức Thí nghiệm

Nền phân bón cố định trên toàn bộ diện tích: $20.000\text{ kg Phân chuồng hoai mục/ha} + 60\text{ kg }\text{P}_2\text{O}_5\text{ (Supe Lân)} + 150\text{ kg }\text{K}_2\text{O (Kali Clorua)}$.

  1. Công thức 1 (CT1): Bón $83\text{ kg N/ha}$ (tương đương $3\text{ kg N/sào BB}$) $\rightarrow$ Lượng phân Urê $\text{CO(NH}_2)_2$ thương phẩm: $180{,}11\text{ kg/ha}$.
  2. Công thức 2 (CT2 - Đối chứng): Bón $139\text{ kg N/ha}$ (tương đương $5\text{ kg N/sào BB}$) $\rightarrow$ Lượng phân Urê $\text{CO(NH}_2)_2$ thương phẩm: $301{,}63\text{ kg/ha}$.
  3. Công thức 3 (CT3): Bón $194\text{ kg N/ha}$ (tương đương $7\text{ kg N/sào BB}$) $\rightarrow$ Lượng phân Urê $\text{CO(NH}_2)_2$ thương phẩm: $421{,}00\text{ kg/ha}$.
  4. Công thức 4 (CT4): Bón $250\text{ kg N/ha}$ (tương đương $9\text{ kg N/sào BB}$) $\rightarrow$ Lượng phân Urê $\text{CO(NH}_2)_2$ thương phẩm: $542{,}50\text{ kg/ha}$.

Phương pháp nghiên cứu và Xử lý số liệu

Số liệu sinh học nông học được thu thập theo chu kỳ 7 ngày/lần đối với chiều cao cây và số lá; 3 ngày/lần đối với số hoa. Toàn bộ số liệu thực nghiệm được xử lý phân tích phương sai (Balanced ANOVA) trên phần mềm chuyên dụng sinh trắc nông nghiệp IRRISTAT phiên bản 5.0Microsoft Excel 2016.

Mô hình tuyến tính của thí nghiệm phân tích phương sai: $$Y_{ij} = \mu + \alpha_i + \epsilon_{ij}$$ Trong đó: $Y_{ij}$ là giá trị quan trắc của công thức thứ $i$ tại lần lặp thứ $j$; $\mu$ là giá trị trung bình tổng thể; $\alpha_i$ là ảnh hưởng của liều lượng đạm thứ $i$ ($i = 1..4$); $\epsilon_{ij}$ là sai số ngẫu nhiên thỏa mãn phân phối chuẩn $N(0, \sigma^2)$. Độ tin cậy được kiểm định qua giá trị $F\text{-ratio}$, kiểm định sai khác nhỏ nhất có ý nghĩa $\text{LSD}_{0{,}05}$ và hệ số biến động $\text{CV}%$.


Implementation và kết quả

Quy trình kỹ thuật và Thuật toán tính toán dinh dưỡng

Quy trình triển khai gồm 4 giai đoạn sinh trưởng chính: Xử lý hạt giống/Ươm bầu $\rightarrow$ Bén rễ hồi xanh $\rightarrow$ Phân hóa thân lá và ra hoa $\rightarrow$ Phát triển quả và thu hoạch.

def calculate_fertilizer_schedule(area_ha, n_rate_kg_ha, p2o5_rate_ha=60, k2o_rate_ha=150):
    """
    Thuật toán quy đổi dinh dưỡng nguyên chất sang khối lượng phân bón thương phẩm
    và phân bổ theo từng giai đoạn sinh trưởng của bí ngồi LS 188.
    """
    UREA_N_CONTENT = 0.46
    SUPER_PHOSPHATE_P2O5_CONTENT = 0.165
    KCL_K2O_CONTENT = 0.60
    ORGANIC_MANURE_KG_PER_HA = 20000

    # Tổng lượng thương phẩm (kg)
    urea_total = (n_rate_kg_ha / UREA_N_CONTENT) * area_ha
    super_p_total = (p2o5_rate_ha / SUPER_PHOSPHATE_P2O5_CONTENT) * area_ha
    kcl_total = (k2o_rate_ha / KCL_K2O_CONTENT) * area_ha
    manure_total = ORGANIC_MANURE_KG_PER_HA * area_ha

    schedule = {
        "Bon_lot": {
            "Phan_chuong_kg": manure_total,
            "Supe_lan_kg": super_p_total,
            "Ure_kg": urea_total * 0.20,
            "KCl_kg": kcl_total * 0.20
        },
        "Thuc_1_Hoi_xanh": {
            "Ure_kg": urea_total * 0.20,
            "KCl_kg": kcl_total * 0.20
        },
        "Thuc_2_4_5_la": {
            "Ure_kg": urea_total * 0.30,
            "KCl_kg": kcl_total * 0.30
        },
        "Thuc_3_Nuoi_qua": {
            "Ure_kg": urea_total * 0.30,
            "KCl_kg": kcl_total * 0.30
        }
    }
    return schedule

Diễn biến sinh trưởng và Động thái tăng trưởng chiều cao

Điều kiện khí hậu vụ Đông Xuân 2017–2018 tại Xuân Mai có nhiệt độ trung bình từ $12{,}76^\circ\text{C}$ (tháng 1/2018) đến $22{,}34^\circ\text{C}$ (tháng 3/2018), độ ẩm không khí trung bình $60{,}1% - 87{,}84%$. Thời gian từ gieo đến nảy mầm của cả 4 công thức đạt mức đồng đều tuyệt đối 8 ngày với tỷ lệ nảy mầm đạt $91{,}67%$.

Thời gian từ gieo đến ra hoa có xu hướng kéo dài khi tăng liều lượng đạm: CT1 (40 ngày) $<$ CT2 (42 ngày) $<$ CT3 (44 ngày) $<$ CT4 (45 ngày). Thời gian thu đợt quả đầu tiên lần lượt là CT1 (48 ngày), CT2 (50 ngày), CT3 (53 ngày) và CT4 (55 ngày).

Bảng theo dõi động thái tăng trưởng chiều cao cây (cm)

Công thức 7 ngày 14 ngày 21 ngày 28 ngày 35 ngày 42 ngày 49 ngày 56 ngày (Cuối)
CT1 (83 kg N) 1,42 2,30 4,30 7,10 8,88 11,97 14,67 16,68
CT2 (139 kg N - ĐC) 1,63 2,77 4,60 7,30 9,20 13,28 17,40 19,30
CT3 (194 kg N) 1,67 2,93 4,47 7,40 9,47 13,33 16,39 18,63
CT4 (250 kg N) 1,80 2,87 5,00 7,37 10,47 14,13 17,32 19,47
$\text{LSD}_{0{,}05}$ 1,30 3,50 4,80 5,20 1,50 4,10 5,50 6,00
$\text{CV}%$ 14,4 10,6 7,6 5,9 7,9 15,7 17,0 16,5

Động thái phân hóa số lá trên cây (lá/cây)

Tốc độ ra lá tăng nhanh nhất ở giai đoạn 35–49 ngày sau trồng. Tại thời điểm 56 ngày sau trồng, số lá thật trung bình trên cây đạt cực đại ở CT4 ($24{,}13\text{ lá}$), kế tiếp là CT3 ($22{,}60\text{ lá}$), CT2 ($20{,}87\text{ lá}$) và thấp nhất ở CT1 ($19{,}06\text{ lá}$). Sai khác giữa CT3, CT4 so với CT1 đạt ý nghĩa thống kê ở mức $\text{LSD}_{0{,}05} = 1{,}45$ ($\text{CV}% = 3{,}9%$).

Đánh giá Năng suất và Các yếu tố cấu thành năng suất

Năng suất là chỉ tiêu quyết định giá trị kinh tế của quy trình bón đạm. Kết quả tổng hợp các chỉ tiêu năng suất cá thể, năng suất lý thuyết và năng suất thực thu được kiểm định nghiêm ngặt qua mô hình ANOVA:

                  SO SÁNH NĂNG SUẤT THỰC THU (NSTT - tạ/ha)

Bảng tổng hợp các yếu tố cấu thành năng suất

Chỉ tiêu theo dõi CT1 (83 kg N) CT2 (139 kg N - ĐC) CT3 (194 kg N) CT4 (250 kg N) $\text{LSD}_{0{,}05}$ $\text{CV}%$
Số quả trung bình/cây (quả) 6,40 6,87 7,00 7,13 1,01 6,1
Trọng lượng quả TB (g/quả) 566,00 674,00 764,67 719,33 72,02 8,6
Chiều dài quả (cm) 21,33 22,14 24,40 22,43 2,06 5,3
Đường kính quả (cm) 21,48 22,60 25,30 22,70 2,07 5,2
Năng suất cá thể (g/cây) 362,29 461,02 535,34 511,78 57,99 7,2
Năng suất lý thuyết (tạ/ha) 2.515,91 3.201,55 3.717,61 3.554,00 402,75 7,2
Năng suất thực thu (tạ/ha) 1.826,67 2.297,00 2.669,33 2.550,16 297,59 7,4

Ghi chú thống kê: Tại mức đạm $194\text{ kg N/ha}$ (CT3), trọng lượng quả trung bình đạt $764{,}67\text{ g}$, vượt trội hơn hẳn so với CT2 ($674{,}00\text{ g}$) và CT4 ($719{,}33\text{ g}$). Khi tăng lượng đạm lên $250\text{ kg N/ha}$ (CT4), số quả có tăng nhẹ ($7{,}13\text{ quả/cây}$) nhưng trọng lượng quả và kích thước quả bị sụt giảm do hiện tượng cạnh tranh dinh dưỡng giữa thân lá và cơ quan sinh sản, dẫn đến NSTT của CT4 ($2.550{,}16\text{ tạ/ha}$) thấp hơn CT3 ($2.669{,}33\text{ tạ/ha}$).


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

Đột phá khoa học và Đóng góp học thuật

  1. Xác định ngưỡng giới hạn bão hòa đạm: Nghiên cứu đã chứng minh bằng thực nghiệm rằng mối quan hệ giữa liều lượng đạm và năng suất thực thu của giống bí ngồi xanh LS 188 tuân theo hàm bậc hai (đường cong năng suất cực cận). Mức $194\text{ kg N/ha}$ là điểm uốn tối ưu sinh học; vượt qua ngưỡng này ($250\text{ kg N/ha}$), hiệu suất sử dụng nitơ (Agronomic Nitrogen Use Efficiency - NUE) suy giảm rõ rệt.
  2. Khắc phục hiện tượng thoái hóa giống: Cung cấp thông số hình thái học chuẩn (chiều dài quả $24{,}40\text{ cm}$, đường kính $25{,}30\text{ cm}$, quả trụ xanh đậm đặc ruột), giúp phân biệt và hạn chế tình trạng thụ phấn tạp giao gây xơ hóa, biến dạng quả ở nông hộ tự để giống.
  3. So sánh với các nghiên cứu tiền nhiệm:
                    SO SÁNH VỚI CÁC CÔNG TRÌNH NGHIÊN CỨU TIỀN NHIỆM

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

Quy trình thâm canh chuẩn mở rộng sản xuất

Hiệu quả kinh tế và Phân tích Lợi nhuận (ROI)

Hiệu quả kinh tế được tính toán chi tiết dựa trên tổng chi phí sản xuất (giống, phân hữu cơ, phân khoáng, thuốc BVTV, nhân công) và doanh thu thực tế từ nông sản tại giá bán bình quân $21.500\text{ đ/kg}$.

Mô hình tính toán: $$Q = A - (B + C)$$ Trong đó:

  • $Q$: Lợi nhuận thuần (Lãi thuần).
  • $A$: Tổng giá trị doanh thu từ sản lượng thực thu.
  • $(B + C)$: Toàn bộ chi phí vật tư và nhân công đầu tư trên 1 ha.

Bảng hạch toán kinh tế trên 1 ha bí ngồi xanh LS 188 (Đơn vị: 1.000 VNĐ)

Khoản mục chi phí & Thu nhập CT1 (83 kg N) CT2 (139 kg N - ĐC) CT3 (194 kg N - Tối ưu) CT4 (250 kg N)
Chi phí hạt giống 15.000 15.000 15.000 15.000
Phân chuồng hoai mục (20 tấn) 10.000 10.000 10.000 10.000
Supe Lân (364 kg) 1.150 1.150 1.150 1.150
Kali Clorua (250 kg) 2.500 2.500 2.500 2.500
Phân Đạm Urê 1.801 3.016 4.210 5.425
Thuốc BVTV & Nhân công 18.000 18.000 18.500 19.500
TỔNG CHI PHÍ ($B+C$) 48.451 49.666 52.360 53.575
Tổng doanh thu ($A$) 392.734 493.855 573.906 548.284
LÃI THUẦN ($Q$) 344.283 444.189 521.546 494.709

Hiệu quả gia tăng: Áp dụng công thức 3 ($194\text{ kg N/ha}$) mang lại lãi thuần cao hơn công thức đối chứng CT2 là $77.357.000\text{ đ/ha}$ (tăng $17{,}41%$) và vượt trội hoàn toàn so với mức bón thiếu đạm CT1 ($+51{,}48%$). Ngược lại, khi bón dư đạm ở CT4, chi phí tăng thêm nhưng lợi nhuận ròng giảm do tổn thất năng suất thương phẩm.


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

Hạn chế kỹ thuật

  • Thí nghiệm mới thực hiện trong phạm vi vụ Đông Xuân tại tiểu vùng khí hậu Xuân Mai, chưa đánh giá phản ứng của giống trong vụ Thu Đông hoặc trên các loại đất khác như đất cát pha, đất phù sa ven sông.
  • Chưa đo đạc định lượng hàm lượng nitrat ($\text{NO}_3^-$) dư thừa trong thịt quả bằng phương pháp sắc ký ion để đánh giá toàn diện tiêu chí an toàn vệ sinh thực phẩm theo chuẩn VietGAP.

Hướng nghiên cứu tiếp theo

  • Khảo nghiệm đa điểm tại các vùng chuyên canh rau lớn ở miền Bắc (Hải Dương, Bắc Ninh, Mộc Châu - Sơn La).
  • Ứng dụng hệ thống tưới nhỏ giọt kết hợp châm phân tự động (Fertigation) với công nghệ IoT nhằm tối ưu hóa lượng phân bón từng ngày theo nhu cầu bốc thoát hơi nước ($ET_c$).

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

                              CÁC NHÓM HƯỞNG LỢI TRỰC TIẾP

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

1. Điều kiện kỹ thuật quan trọng nhất để triển khai công thức đạm CT3 là gì?

Cần chuẩn bị đất tơi xốp, lên luống cao 25–30 cm để chống úng cục bộ, bón lót toàn bộ phân chuồng hoai mục ($20\text{ tấn/ha}$) cùng phân lân. Chia lượng đạm Urê ($421\text{ kg/ha}$) đúng 4 đợt theo tỷ lệ $20% - 20% - 30% - 30%$, không được bón dồn một lần.

2. Dấu hiệu nhận biết bí ngồi bị thừa đạm và cách khắc phục?

Cây thừa đạm có thân lá xanh đen đậm, vươn dài mềm yếu, xuất hiện nhiều hoa đực nhưng ít đậu quả cái, dễ bị rệp và bọ trĩ tấn công. Cách xử lý: Ngừng bón đạm ngay lập tức, tưới rửa bớt phân thừa, bổ sung phân Kali Clorua kết hợp phun Canxi-Bo qua lá để kích thích phân hóa mầm hoa.

3. Có thể thay thế phân chuồng hoai mục bằng phân hữu cơ vi sinh không?

Hoàn toàn được. Trường hợp thiếu nguồn phân chuồng truyền thống, có thể thay thế bằng $800 - 1.000\text{ kg phân hữu cơ vi sinh/ha}$ bón lót cùng với Supe lân trước khi lên màng phủ nông nghiệp.

4. Tại sao phải thụ phấn bổ sung nhân tạo vào buổi sáng?

Bí ngồi là cây đơn tính cùng gốc (hoa đực và hoa cái riêng biệt). Trong điều kiện vụ đông xuân lạnh và ít côn trùng thụ phấn, việc ngắt hoa đực quét hạt phấn lên đầu nhụy hoa cái lúc 7:00 – 10:00 sáng sẽ nâng tỷ lệ đậu quả từ dưới $50%$ lên trên $90%$, hạn chế tình trạng teo quả non.

5. Thời gian thu hồi vốn và khả năng nhân rộng của mô hình?

Cây bí ngồi xanh LS 188 có chu kỳ sinh trưởng cực ngắn (chỉ 50–55 ngày từ gieo hạt đến khi thu hoạch quả đợt đầu). Thời gian thu hoạch kéo dài 20–25 ngày. Nông dân có thể thu hồi toàn bộ vốn đầu tư và đạt lợi nhuận ngay trong vòng chưa đầy 3 tháng canh tác.


Kết luận

Đề tài "Nghiên cứu ảnh hưởng của liều lượng đạm đến sinh trưởng, phát triển và năng suất của cây bí ngồi xanh LS 188 vụ đông xuân 2017-2018 tại vườn ươm Trường Đại học Lâm nghiệp" đã giải quyết trọn vẹn bài toán xác lập định mức phân bón đạm cho cây bí ngồi tại khu vực phía Bắc.

Kết quả khẳng định Công thức 3 ($194\text{ kg N/ha}$, tương đương $7\text{ kg N/sào Bắc Bộ}$ + $60\text{ kg }\text{P}_2\text{O}_5 + 150\text{ kg }\text{K}_2\text{O} + 20.000\text{ kg phân chuồng}$) là công thức tối ưu toàn diện: đạt số quả trung bình $7{,}00\text{ quả/cây}$, khối lượng quả $764{,}67\text{ g}$, năng suất thực thu đạt kỷ lục $2.669{,}33\text{ tạ/ha}$ và mang lại hiệu quả kinh tế cao nhất. Đề nghị các cơ quan khuyến nông, hợp tác xã và bà con nông dân ứng dụng quy trình kỹ thuật này vào sản xuất đại trà nhằm nâng cao giá trị gia tăng và thu nhập bền vững trên một đơn vị diện tích canh tác.